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PIERRE CURIE 
(On the Fiftieth Anniversary of His Death) 


April 19, 1956, is the 50th anniversary of the 
tragic death of the outstanding French scientist 
mmerre iCurie: 

Pierre Curie was born in 1859 in Paris into the 
family of a doctor. His early education for middle 
school was obtained at home. His striking ability 
at mathematics, especially solid geometry, en- 
abled him to obtain his baccalaureate in natural sci- 
ences at the age of sixteen; he became a student in 
the physics department at the Sorbonne. 

At the age of 19 he became an assistant in the 
department of experimental physics at the Univer- 
sity of Paris; this marked the start of his scienti- 
fic career. His first paper was on the determina- 
tion of the wavelengths of thermal radiation. Short- 
ly afterwards he took up crystallography, with em- 
phasis on crystal symmetry. Crystallography was 
his sole interest from 1879 to 1891, especially crys- 
tal physics; later, he worked on radioactivity, but 
remained interested in crystals. 

Pierre and his elder brother, Jacques (a miner- 
alogist),concluded that deformation should produce 
electrical polarization in certain crystals. This 
was strikingly confirmed shortly afterwards; in 
1880 they discovered piezoelectricity, which they 
studied in great detail on quartz crystals. They es- 


tablished Curie's first and second laws for the 
longitudinal and transverse piezoelectric effects in 
quartz, which remain accepted even today. They 
also observed the inverse effect predicted by Lipp- 
man, namely deformation of a crystal linearly de- 
pendent on the applied electric field. 

The piezoelectric effect proved to be of con- 
siderable scientific importance; the brothers Curie 
used a quartz plate to measure small charges and 
currents, and the device they developed for this 
purpose proved of value in the work on radioactiv- 
TA Piezoelectricity is now widely used in radio, 
electronics, hydroacoustics, and instruments gen- 
erally. 

Pierre Curie's papers on symmetry form a 
major part of his contribution to crystallography; 
these he wrote mainly from 1883 onwards, after he 
had been appointed professor in the school of physics 
and chemistry. He derived the limiting point-sym- 
metry groups for crystals, which correspond to the 
symmetries of homogeneous anisotropic media with 
axes of infinite order. The derivation of these 
groups enabled him to give a rational classification 
of crystals on the basis of the possible presence of 
various physical properties. He linked his work on 
symmetry closely to these physical properties, as 
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when he stated that "I consider that it would be of 
interest in the study of physical phenomena to bear 
in mind the considerations of symmetry so familiar 
to the crystallographer." 

His work over many years led to the formula- 
tion of the general principle of symmetry, which he 
gave in its final form in 1893-95, as follows: "When 
certain causes give rise to known consequences, the 
symmetry elements of the causes must be repro- 
duced in the consequences. When known consequen- 
ces contain a known dissymmetry, this latter must 
lie in the causes giving rise to the consequences. 
Propositions contrary to these two are incorrect, 
at least in practice,i.e., that consequences can be 
more symmetrical than the causes giving rise to 
them." 

Pierre Curie also made a major contribution 
to the theory of crystal growth, which had been 
founded by Gibbs previously, though this was un- 
known to him at the time. The basic concept was 
that a crystal in equilibrium with a solution must 
have a shape such as to correspond to minimum 
surface energy subject to constant volume. 

In 1891, Pierre Curie began a long series of 
studies on the magnetic properties of bodies as 
functions of temperature; these formed the basis of 
a dissertation submitted in 1895. His results on 
paramagnetic, ferromagnetic, and diamagnetic 
bodies were of considerable importance; the results 
for paramagnetic bodies are known under the name 
Curie's law (susceptibility inversely proportional 
to the absolute temperature). He also observed that 
iron has a special point above which it ceases to be 
ferromagnetic (the Curie point). Similar relation- 
ships have since been deduced for ferroelectrics. 

He also participated in researches by his wife, 
Maria Curie-Sklodowska, on radioactivity; this led 
to their joint discovery of polonium and radium in 
1898, and in 1899 to the first observation of colors 
produced by irradiation (in glass and porcelain 
acted on by radium rays). They also showed that 
these rays affect the cells of organisms. In 1903, 
Pierre Curie observed the spontaneous production 
of heat by radium salts, which proved that radio- 
activity is accompanied by the release of much 
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energy; this pointed the way to the idea that radio- 
activity is a result of atomic transformations. The 
Curies made a detailed study of this radiation, and 
showed that the decay was not detectably dependent 
on the external conditions. This led Pierre Curie 
to state that decay constants could serve as time 
standards and could be used to determine the ages 
of rocks. In 1903, the Curies and Becquerel were 
awarded the Nobel prize for their work on radioac- 
tivity. 

In 1904, a new department was set up for 
Pierre Curie in the University of Paris (a resolu- 
tion of the French parliament was required to do 
this); his post in the school of physics and chemistry 
was taken over by one of his students, Paul Langevin, 
who later also became a famous physicist. In 1905, 
Pierre Curie was made a member of the Paris 
Academy of Sciences. 

The work on radioactivity was continued after 
his death (April 19, 1906) by his wife, who re- 
placed him as professor at the Sorbonne. Later, 
this department was headed by their daughter,Irene 
Joliot-Curie, who died in March of this year. The 
work of Frederic and Irene Joliot-Curie together 
led to a most important discovery, artificial radio- 
activity, which further increased the scientific 
fame associated with the name of Curie. 

The name of Pierre Curie is eternally recorded 
in the annals of science. In addition to being an 
energetic and exceptionally capable scientist, Pierre 
Curie was a warm and humane man. A good indica- 
tion of this is given by his words on receiving the 
Nobel prize: "It is not difficult to foresee that radi- 
um in the wrong hands could be extremely danger- 
ous; this raises the question whether it is best for 
mankind to know the secrets of nature, whether 
man is sufficiently adult to use them properly, or 
whether such knowledge will only do him harm. 
Nobel's discovery is here characteristic. Powerful 
explosives have enabled men to perform great 
works, but they are also a fearful means of destruc- 
tion in the hands of great evildoers who drive men 
to war. I am one of those who, like Nobel, consider 
that all new discoveries will, in the end, bring man- 
kind greater happiness than harm." 
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Formulas are deduced for the intensities in texture and polycrystalline patterns 
from the distribution function of the angles taken up by the crystallites. This 
function has two components in the case of a texture: adensity of n/27 con- 
stant at all angles (n is the number of crystallites in the specimen) and a dis- 
tribution f (~) of the axes of the crystallites around the texture axis. An ex- 
pression is derived for the integral intensity of arcs, this being independent 
of the form of f(%, and also one for the intensity at the center of the arcs, 


which is dependent on f(¢). 


The distribution-function concept ig used to dis- 


cuss the intensities of the rings given by apolycrystalline aggregate. 


Angular Distribution Function for 
Crystallites 


Part 1 [1] dealt with the intensities of spot pat- 
terns from ideal or mosaic single crystals. The 
reflection intensities were deduced by reference to 
an angular distribution function for the crystallites 
(and hence for the reflecting planes in them). The 
same concept is used here. Let f(a) be the angu- 
lar density of the normals to 2 given plane; the 
range @ to ~+do contains f (z)do crystallites, and 
the integral over all possible positions equals the 
number of crystallites in the specimen: 


\ f(2)d2=n=V’ V, (L) 


where V is the volume of 2 crystallite and V' is the 
volume of the specimen.* 

It is convenient in calculating the intensity to 
express the distribution over solid angles in terms 
of components over plane angles; condition (1) ap- 
plies to each such component. 

Texture patterns are widely used in electron 
diffraction [2,3], specimens with platy textures be- 
ing most common, in which the crystallites have 
some major crystallographic face parallel to the 
substrate but are disordered in azimuth. The direc- 
tion of the normal to this face is the texture axis. 


Needle textures are relatively rare [4]; in these the 
texture axis is a direction common to the axes of all 
the needles, which are randomly rotated around that 
direction. The axis of a platy texture is normal to 
the substrate and coincides with one of the axes of 
the reciprocal lattice; the axis of a needle texture 
lies in the plane of the substrate and coincides with 
one of the axes of the atomic lattice [3]. The re- 
ciproca] lattice for 4 texture of either type is a set 
of concentric rings and has cylindrical symmetry; 
the texture axis is an axis of infinite order. 

We have seen in Part 1 that the distribution 
function can be expressed in terms of components 
for plane angles that correspond to rotation around 
certain specified axes. 

The first such component is f (4), the density 
around the texture axis (range 277); this is constant 
and is given by (1) as 


an 


f(a) \ day =n, {(a)) = nf 2m, @) 


4G 
A. veal specimen is imperfect, so the special 
directions in the crystallites (all parallel in the 
ideal case) are spread over a certain solid angle 
(Fig. 1a), whose magnitude defines the deviation 


ys is, 
tron beam, and n is the number of crystallites in that volume, 


of course, the volume of the specimen traversed by the elec- 
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b Texture axis 


Fig. 1. a) Structure of a textured specimen, and b) com- 
ponent distribution functions f (a,) and f @») for the angles 
of the normals with respect to different planes. 


from perfection. This spread is responsible for the 
conversion of the reciprocal lattice rings (infinite- 
ly thin) into spherical belts (Fig. 1b), so the re- 
flections corresponding to the latter take the form 
of short ares [2,3]. The spread is not uniform, 


i.e., is characterized by a second distribution f (w). 


It is convenient to replace this spread in the solid 
angle a by f(a), the projection of this distribution 
on some plane angle a, (the cylindrical symmetry 
makes the choice arbitrary), which (Fig. 1b) de- 
fines the form of the nodes on the texture axis and 
of the cross sections of the other ring nodes. Let 
a'be the angular range in which f (@) is finite; then 
formula (21) of [1] shows that the density at the peak 
is about twice the mean of f (a), i.e. , that 


; Aaa (2, 
1 (4s)max © 257r\f(%2)d—q=f=t, ©) 


where @ = a'/2 is the effective angular width of 
SF (@). 

It is a matter of no importance how the inten- 
sity in an arc is distributed if we are interested 
only in the integral intensity; the angular spread 
around the texture axis is then unimportant, and 
J (@y) should not appear in the formula for the in- 
tegral intensity. But another approach is to meas- 
ure the intensity at the center of the are (i.e.,over 
a small area); this local intensity will be depend- 
ent on f (9). 
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Fig. 2. Scheme for calculating the 
intensity of a reflection for a straight 
texture. 


The distribution function for a polycrystalline 
specimen is as follows. The crystallites take upall 
spatial orientations at random, and the normals to 
any reflecting plane are uniformly distributed 
throughout a solid angle of 47. The distribution 
function with respect to solid angle is then constant 
at n/47. 


Integral Intensities in Texture 
Patterns? 


Straight Textures. _ By this is meant 
the pattern from a needle texture with the specimen 
at right angles to the beam (g = 0), and also the 
pattern from a platy texture in reflection. The tex- 
ture axis is perpendicular to the beam in bothcases. 
Such texture patterns would also occur with platy 
textures in transmission for g = 90°, but this can- 
not be realized. 

The integral intensity is calculated precisely 
as in [1] for mosaic specimens, by taking the inte- 
gral of the interference function over two directions 
in the plane of the screen, which in this case coin- 
cide with directions h, and hy in the reciprocal lat- 
tice, which is represented on a scale LA on the 
screen (Fig.2). Then formula (11) of [1] applies: 


Tayeetd oe 
hkl 0 Q | 2 a (4) 
Formula (138) of [1] gives 
sin27A sh. 
G=\7(@) “Graghgye 1% (5) 


which is defined by the integral of the interference 
function in the third direction, as modified by the 
distribution function. Only the first component in 
(2) is important for the integral intensity of an arc, 


See [5] for some preliminary evidence on this. 
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so by da in (5) we should understand an element of 
rotation around the texture axis. This elementary 
angle is dependent solely on the distance of the re- 
flection from the texture axis in the pattern, which 
corresponds in reciprocal space to the horizontal 
component of the reciprocal-lattice vector Hpi 
(Fig.2), if, as is usual, we denote the texture axis 
by c*. This horizontal component is Hpyko for an 
orthogonal lattice. The formulas may be applied to 
oblique lattices if by Hyko (or dpc) is understood 
the corresponding horizontal component. Then, 


dhs- Li 
= Bigg = nase Soot (6) 


where R is the distance of the reflection from the 
vertical axis of the pattern (a quantity that is meas- 
ured). From (2), since the integral of the interfer- 
ence function is, from (7) of [1], equal to A;/a%, 
etc., 


CS A3n any 0 A3n Lr 
ae on a eR 


(7) 


It is of interest to compare (7) with the analo- 
gous expression (16) of [1], where dnko/ 27 is re- 
placed by dy j/@, the other first factors being the 
same. The replacement of dpi by dpyyg reflects 
the fact that the angular spread in a mosaic is the 
same for all directions,i.e.,occurs around the 000 
point in the reciprocal lattice, whereas, in a tex- 


ture it occurs around a special direction (the tex- 
ture axis). The small quantity @ (the effective 
angular spread in a mosaic, about We radian) is re- 
placed by a large one, 27, so the reflection arises 
from far fewer crystals (about '49). 

Only a very small part is effective out of the 
complete circle produced from the node; the inten- 
sities in texture patterns are therefore very small 
fractions of the initial intensity, so kinematic scat- 
tering applies very closely provided, of course, 
that the blocks do not exceed the critical size. Dy- 
namic effects may occur in rare cases, but these 
are far more common in mosaic materials. The 
reason is that a textured or polycrystalline speci- 
men must contain much smaller blocks than doesa 
mosaic one of the same volume if the distribution 
over all orientations is to be uniform. 

Substituting (7) into (4), and using (2), we have 


a Lrro ip 


2 
Inna = Jon? | | Festal (8a) 
and 
Tne 2 |> ? ako 
Fase wale [| Sete? (8b) 


where V', S, and t are the volume, area, and thick- 
ness of the film, p being the repeat factor (see be- 
low). R may be measured directly on the pattern, 


in which case the dpxg of (8a) and (8b) should be re- 
placed by LA/R. 


Fig. 3. Calculation of the intensity of a reflection in an oblique-texture 
pattern. a) Inclined section of ring; b) form of this section on direction]; 


c) form in ring on arrow II. 
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The 00/ reflections (on the texture axis) do not 
obey these formulas, because they have R equal to 
0 and dpkg = dgo9 = *: their integral intensities 
may be calculated by noting that the angular spread 
of these planes over 27 is unimportant (the planes 
reflect without relation to their azimuthal position). 
But they appear in the reflecting position on ac- 
count of f (@), @, here being reckoned perpendicu- 
lar to the screen; but this means that formula (18) 
of [1] or a mosaic single crystal) is applicable to 
the 001 reflections, so the dpKo/27 of (8a) and (8b) 
should be replaced by dggj/@ , where a is the ef- 
fective angular spread of the texture axis; p = 1. 

Oblique Texture. An oblique-texture 
pattern is obtained when a platy texture is inclined 
at an oblique angle ¢ (usually about 60°); the reflec- 
tions lie on ellipses [2,3]. Figure 3a illustrates the 
intersection of the rings in the reciprocal lattice by 
a plane (the reflection sphere). The following 
points are important here. The integration for (4) 
was taken over two directions in the plane of the 
screen coincident with the h, and hy directions of 
the reciprocal lattice. An inclined specimen gives 
an element of length in the vertical direction on the 
pattern (Fig. 3) of dhy as dh,/sin g, so the inte- 
gral of the interference function in this direction 
will not be A,/aj fone of the factors in (4)], but 
A,/ai?sin ¢@. The horizontal direction also does not 
coincide with h,; Fig. 3c shows that, if the hori- 
zontal coordinate R' of the reflection is known, to- 
gether with the minor semiaxis R of the ellipse 
(i.e., the radius of the ring in the reciprocal lat- 
tice), then the element of length dhp: for this direc- 
tion is given by dhp:/dh, = R/R', and so the inte- 
gral of the interference function in this direction 
will take the value (A,/a$) + R/R'. Allowance may 
be made for the spread around the texture axis by 
integrating with respect to da,; (7) remains applic- 
able. Then we must insert factors of R/ R' and 
1/ sing in (8a) and (8b) to give 


@M |? Lx- 
Tru = Fo? | | Var ace ee) 
and 
Tk ® |? Li: 
See con || ton abl Rs © 
RSs lies) 2nR’ sin 9° |(9b) 


But sing = 1 for g = 90°; R' becomes R and 
(9a) and (9b) become (8a) and (8b), which is a par- 
ticular case of (9a) and (9b). To use (9a) and (9b), 
we need to know ¢ and R', which are not constant 
for the reflections on a given ellipse, but decrease 
as l increases. We can reduce (9a) and (9b) to a 
form analogous to (8a) and (8b) by introducing the 
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fictitious interplanar distance 


hye 
TE Resinie!@ 


ft (10) 
which is greater than dpxo for an oblique texture, 
but becomes equal to it at p = 90°. 

Reflections on the vertical axis of the pattern 
do not obey (9a) and (9b), because R' is 0; in prin- 
ciple, these intensities could be calculated from 
J (@,), as for the 001 reflections of straight tex- 
tures, but there is no need for this complicated cal- 
culation, because rotation of the specimen will re- 
move the reflection from the region with R'= 0 and 
split it into two arcs with R' ~ 0. 

We may convert from relative integral intensi- 
ties for arcs on texture patterns to relative moduli 
|@| via 

ahet fue VIET) pe 


hkl | = hkl* 


(11) 


where R! varies along an ellipse for an oblique tex- 
ture but is constant for the vertical direction for a 
straight texture, in which case it can be replaced by 
R ~ Hnko = 1/ dhko.- 


Local Intensity in a Texture Pattern 


The local intensity is the best quantity to meas- 
ure if the arcs are very broad, namely, I}; 7 atthe 
center of an arc. Consider first the case gy = 90°, 
where the angular width is governed by the form of 
J(@). The effective angular width of an arc (which 
is the effective angular width of the spherical belts 
of Figs. 1 and 2 for a straight texture) is a © a'/2, 
in which a'is the real angular width; see (3). 

Let A be the length of the small part of the arc 
used (the slit length in microphotometry). The in- 
tegral intensity is distributed, in effect, over a 
length ra, in which r is the radius of the arc onthe 
pattern. The intensity localized in A is then 


A 


rs 
Tet = Txt — 
Akl hkl > 


(12) 
with lhkl = LA/dp 7 » sO 


- Tnx! 


i ee 
hkl ia pores 5 aa 


(13) 


The local intensity at the center of an arc is 
then 


aN Linx Ink 
Tani = Jo? gee (14) 


@D [2 
3] 


This is so not only for ¢ = 90°, but also for any 
other y, and so is general, as we now show. 

The effective angular width a" of an are on an 
oblique-texture pattern is not equal to a (Fig. 4), 
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Fig. 4. Determination of a". 


but is given by 
(15) 


where x is a complicated function of ¢ [5]. The 
spherical right-angled triangle ABC of Fig. 4 im- 
plies that cos Y= cos x sin? ; cos? = coswsing; 
these transform to sin ysin? = singsinw. But 
sinw = Hp:/HpK] and sin’ = Ayko /Hpk], so (15) 
and (10) give 


Ce = Ch) SING, 


Any 


: Heo e anno (16) 
Then I}, is, as in (12), 
Ls ee Ae Adio 
Tit = Dnt =a = Tix a (17) 


We substitute (9a) and use (10) to get (14) again, be- 
cause d' and sing cancel, and the dependence on the 
angle of rotation is lost. This result is entirely in 
order: a measurement of the local intensity is es- 
sentially one of the density of the interference func- 
tion for a reciprocal-lattice ring, and this cannot 
be dependent on the angle of rotation of the speci- 
men. The local intensity can also be measured for 
hkO reflections without inclining the specimen (~= 0) 
when the reflections take the formof rings;dpkKodhk 
in (14) then becomes dt i : 

The relative local intensity is then converted to 
| 6] via 


Ores VG eiden din, p. 


See [5] on the question of p for the texture me- 
thod; a general technique for finding this factor is 
examination of the point reciprocal lattice for the 
crystal with allowance for the | ®),;| for each node 
(i.e., consideration of the structure factor) to- 
gether with calculation of the number of point nodes 
that fuse into a ring node for the given orientation 
of the texture axis. 


(18) 
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Fig. 5. Scheme for calculation 
of the local intensity from a 
polycrystalline aggregate. 


Note should be made of the cases characteristic 
of the hemihedral classes in which a ring node con- 
tains crystallographically distinct point nodes, for 
the electron-diffraction pattern can give only the 
mean value (|®p, 7I? +|®p, 7|?)/ 2. 


Intensities from Polycrystalline 
Films 


The distribution over solid angles is n/47 for 
a polycrystalline film; (4) is found by integration 
‘over two directions in the plane of the pattern. The 
point of interest is usually the local intensity in a 
ring over a short lengthA. The angular width of 
this is a, =A/r=A/Hpy LA (Fig. 5). Integration 
with respect to da, and multiplication by the plane 
angle @») gives us G, which relates the solid-angle 
distribution to the product of the two plane ones: 


A. dhs 


Ang? in2x Ash: 
é — Ie hk SIVTAgs le. ae 
i ee Ge ee ee 
Thus, 
IN Gees Ak 
hkl % 
Se 20) 


Incorporating p (a different value is needed 
here), we have from (19) and (20) that the local in- 
tensity in a length A is 


, ier 
Lk =a Lge Q V GnmLX (21) 
Then, for the relative values, 
[ore =V Tin / dip. (22) 


Formula (21) is the only intensity formula for 
electron diffraction that can be derived from an x- 
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ray powder formula with allowance for the small 
scattering angles for electrons. In this case (and 
this one only), the two situations are identical: a 
scattering volume with randomly distributed crys- 
tallites. 


Conclusions 


The kinematic approximation is used to derive in- 
tensity formulas for the reflections in electron diffrac- 
tion from textures and polycrystalline aggregates.A 
texture pattern containing shortarcs is best measured 
in terms of integral intensities, the|@| being found 
via (9a), (9b), and (11). More extended arcs make 
it preferable to measure the local intensity at the 
center of each arc, the || here being found via 
(14) and (18). The local intensity of a powder ring 
is given by (21). Textured and polycrystalline spe- 
cimens, as used in electron diffraction, have small 
crystallites, so they are much less prone to give 
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dynamic scattering than are mosaicsingle crystals; 
the scattering is usually kinematic. 
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It has been usual to deduce the Laue class andorientation of anunfacetted crystal 
from three patterns recorded on flat film. A new method is proposed for crystals giv- 
ing highly detailed patterns that employs only a single cylindrical film. The 
spherical coordinates of reflections and of points of intersection of zonal curves 


are deduced from a set of curves for p = const and @ = const. 


The results are 


marked in directly on the stereographic projection, which simplifies and short- 


ens the operations. 


Any x-ray study of a structure by oscillation, 
rotation, etc., is preceded by the setting of the 
crystal in the correct position in the goniometer. 
This may be a simple matter if the crystal is eu- 
hedral, for all that is required is optical adjust- 
ment sometimes followed by final adjustment from 
Laue patterns [1]. 

The operation becomes much more difficult if 
the crystal has no faces, or only very poor ones; 
the only means of identifying the crystallographic 
directions may then be to record patterns with poly- 
chromatic radiation. 

The main line of development has so far been 
based on the use of flat films with the crystal at 
rest; but the flat-film method has a major but not 
always obvious, drawback, for it records only dif- 
fraction directions lying at angles between 3-4 and 
60-70° to the primary beam. 

The stereographic projection given by these re- 
flections is far from complete. Umanskii and 
Kvitka [2] have used three Laue patterns with crys- 
tal settings differing by 60°, the three being used 
together in the stereographic projection and in the 
deduction of all reciprocal-lattice nodes. However, 
rather a large volume of work is involved in pro- 
cessing three Laue patterns. 
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Komkov and Frank-Kamenetskii [8] have used 
a single pattern in the zone-occurrence method to 
determine the orientation. We consider that the 
conclusions drawn from a few reflections on a single 
flat film may be ill-defined or wrong; our own ex- 
perience confirms this. 

A crystal giving a complex pattern (minerals, 
inorganic, and intermetallic compounds) can read- 
ily be made to yield more reflections by using a 
cylindrical film with a tube operated at about 40 kV 
and having an anode of Mo, Ag, or W together with 
a beryllium or other such window. 

This combination (fixed crystal, cylindrical 
film) has so far not been used in this country, be- 
cause it has been considered more difficult to con- 
struct a stereographic projection from a cylindric- 
al film; but this is very greatly simplified by a net 
of the curves for p = const and (90 — ~) = const 
drawn up by Rosch [4] as long ago as 1926, and de- 
scribed by Schiebold and Schneider [5]. 

Henry et al. [6] have given a netofthese curves 
on a scale suitable for practical use; Fig. 1 repro- 
duces this in a more convenient form. 

The curves for p = const may be represented 
as the result of projecting onto a cylinder the lines 
of intersection of the Ewald sphere with cones of 


124 


A. I. ZASLAVSKII AND D. L. ROGACHEV 


55 TT 50 a) 35 35.2 hd ge ebin 150 W 95 
LILO PO 
{ aS 
és a i Ws Ax 65 
Wy 
70 Pe leodore i \ \ 70 
i A] \ WEE It +4- i 
85 Wee ae ; 
v Ss eee : 
80 iy 80 


60 


FMA Soe 
ane 


| 


50 


/\ 


0 35 


Ai V0 


My 


ij) \ 
i 


BS 
=== 


pay 
WES KT ah, 

l\esesectagel 
ec 


SE 


Fig. 1. Net of curves for p = const and (90 — gy) = const for a cylindrical film of radius 30 mm. 


semivertex angle p. These lines are somewhat 
similar to the digit 8 in shape, and they retain 
much the same form on projection onto a cylinder. 

The curves p = const may be represented ana- 
lytically as 


a \2 
Pence ere ane 


ace | 


COs 


The second type of curve, (90 — ¢) = const, re- 
sults from projection onto the cylinder of sections 
of the Ewald sphere by a series of vertical planes 
at various angles to the primary beam. 

The (90 — ~) curves may be represented as 


12 


cos es =o) = cos? /1 + (ir) : 


The two sets of curves together give the spher- 
ical coordinates of reflections recorded on the cy- 
lindrical film, and these coordinates enable us to 
construct a stereographic projection, which gives 
a complete picture of the diffraction symmetry, be- 
cause the points fill a large part of the projection. 


If necessary, p and ¢@ can be calculated from 
the coordinates x and y on the film; 


RB 
(1+ B2)—V1i+ Bcos A’ 


ellie 


cosp= - 


V1+ B?—cosA 


tan o == ; 
3 sin A 


where A = x/R, B= y/R, and R is the radius of the 
film. 

The array of spots in the stereographic projec- 
tion gives a reliable identification of the Laue class 
as well as the angular adjustments needed to bring 
the axis of rotation into coincidence with any given 
crystallographic direction. 

We consider it much more convenient to use a 
projection constructed directly on a stereographic 
sphere for this purpose. 

The method may be illustrated by our work on 
the orientation and Laue class of a new mineral, 
which was found and passed to us by M.D. Dorfman 
of the Kola Branch of the Academy of Sciences, 
which so far has been called pink platy Fe silicate. 
Figure 2 shows the Laue pattern recorded on a cy- 
lindrical film with the crystal in a random position; 
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Fig. 2. Pattern from a crystal of “pink platy Fe silicate" recorded on a cylindrical 
film (radius 30 mm) with the crystal in a random position (Mo anode), 


Fig. 3, Scheme for the Laue pattern from a crystal of “pink platy Fe silicate" 
recorded on a cylindrical film; the zonal curves are shown, together with the 
numbering of the spots and the quartering. 


The array on the sphere revealed at once a sym- 
metry plane, with a two-fold axis perpendicular to it; 
A tracing from this scheme was superimposed the crystal therefore belongs to Laue class Cyp,-2/m. 
on the net to give the coordinates of the main points, It is also easy to establish the angular movements 
which are listed in the table. of the sphere needed to bring any crystallographic 
These results were used to mark in the prin- direction to the vertical position, e.g. , the two-foldb 
cipal crystallographic directions on the stereograph- axis (Figs. 5 and 6b). The usual x-ray methods, 
ic projection (Fig.4) and also on a stereographic e.g., that of [1], may be used to perform the final 


sphere (Fig. 6a). adjustment. 


Fig. 3 shows a scheme for this, with the zonal 
curves and numbering of the main points. 
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Fig. 4, Initial stereographic projection for the crystal. 


Fig. 6a. Projection of the crystal on a 
stereographic sphere. 


Conclusions 

A convenient method of determining the orien- 
tation and Laue class is as follows. 

1. The Laue pattern is recorded on a cylindri- 
cal film with the crystal in a random position; the 
zonal curves and points of intersection of these are 
identified. 

2. A net of curves for p and (90 — ¢) is used 
to determine the spherical coordinates of these 
points of intersection, which correspond to major 
directions in the reciprocal lattice. ! 
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Fig. 5. Result of analysis of the initial stereographic projection 
from a crystal of "pink platy Fe silicate" recorded on a cylin- 
drical film. 


Fig. 6b. Projection of the crystal on a 
stereographic sphere with identification 
of a symmetry plane. 


3. These directions and zonal curves are 
entered directly on the stereographic sphere or on 
the stereographic projection. ? 

4. The symmetry elements relating these 
points are deduced, as well as the angular displace- 
ments needed to bring any given crystallographic 
direction into coincidence with the axis of rotation. 


1 analytic computation is best for points with p and g less than 30°, 


* The stereographic sphere can also be used to advantage in the 
three-patterns method, 
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TABLE 1. Data for Stereographic Pro- 
jection Derived from X-ray Pattern of 
Pink Platy Fe Silicate Crystal 


Spot | ° | (90—@) |Quarter 
| 
1 (& 54 I 
2 40 4 ] 
17 47 73 I 
i, 66 The) I 
8 22 10 II 
16 31 56 Il 
9 49 56 Il 
10 84 56 II 
15 q5 69 Il 
18 90 79 II 
4 88 72 1\ 
3) 62 2 IV 
3 59 54 IV 
6 35 3D IV 
19 18 5 Iv 
20 29 5 IT! 
13 BD 54 Il 
2 4A 38 ia 
11 56 55 TW 
14 80 70 Ii! 


5. The adjustment is refined by the method of 
Kvitka et al. 

The work is simplified if the crystal has some 
faces, because the crystal can be given a reason- 
able initial orientation. 

This method is a rapid, convenient, and reli- 
able way of identifying the Laue class and orienta- 
tion for a crystal lacking faces. The method is es- 
pecially convenient if the camera does not allow a 
flat film to be used, especially in Weissenberg 
goniometers. 


We are indebted to M. D. Dorfman for providing 


the crystal for examination. 
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The possible structure types are deduced on the assumption that close packing 
allows only structures in which an A atom is surrounded only by B atoms, 


while each B atom has at least one A atom in its coordination sphere. 


Two 


different models for the structure type are proposed for € -Cu3Sn, which has 
hexagonal close packing and a superlattice, for which the space group and 


repeat distances are known. 


Here we consider the structure types for binary A B atom in the layer may be surrounded by B 


compounds! in which the A and B atoms together 
form a close-packed array; minor distortions are 
neglected. Many intermetallic compounds are of 
this type: all? except Ag3Sb [1] are such that an A 
atom lies at the center of a cube —octahedron of B 
atoms, while each B atom has A atoms in its coor- 
dination sphere. On this basis we assume that 
close packing allows only structures in which an A 
atom is surrounded only by B atoms, while each B 
atom has at least one A atom in its coordination 
sphere. Determination of the possible structure 
types then amounts to determination of the array of 
A and B atoms. 


Let Pp be the reciprocal of the number of A 
atoms in the coordination sphere of a B atom, and 
let PA be the sum of the Pp for all B atoms in the 
coordination sphere of an A atom; Pa is then the 
number of B atoms per A atom. 


Figure 1 shows the possible coordination of the 
A and B atoms in a layer; Pp can be Taw le. or 1. 


atoms alone. An atom with Pp = in a layer also 
has Pp = Je in space, because only B atoms occur 
at the six points 1,2,3,4,5, and 6 in cubic close 
packing and at the three points 1,2, and 3 in the 
layers above and below in hexagonal close packing 
(Fig.2). A B atom with Pp = 1 in a layer, and also 
surrounded by B atoms alone, will have a Pp in 
space not less than 1, because, at most two of the 
six atoms above and below it can be A atoms. AB 
atom with Pp = uy, in a layer will have in space a Pp 
either of '4 (Fig. 3) or ¥4; Pp in space cannot be 
less than vie so Pa in space cannot be less than vA 
* 12 = 3, i.e., there must be not less than three B 
atoms per A atom. Each A atom cannot be com- 
pletely surrounded by B atoms for A:B > 1;3.3 


1T consider structure types, not particular compounds, for the latter 
may be of variable composition, which may not even include the 
composition corresponding to the given structure type. 

>The range of compositions with B: A < 3 is not considered here. 
>The range of compositions with B: A < 3 is not considered here. 


® ie) ie) ® 


Fig. 1. Possible coordination of A atoms (@) and B atoms (©) in a layer, 
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Fig. 2. a) x , layers above and below; O, middle layer; b) 
x, lower layer; O, middle layer; 0, upper layer; @, A atoms. 


Fig. 3, Possible coordination of A atoms (@, x ,™) around a 
B atom for Pg = +4 for (a) hexagonal close packing, (b) cubic 
close packing. a) O, middle layer; x , lower and upper 
layers. b)x , lower layer; O, middle layer; 0, upper layer. 


OOS 


VY2+ e+ Ye + YerYer+Y2-3 


Fig. 4. Possible linking of hexagons of B atoms for Pg = 1) in a layer. 


@, links via vertices; —, links via sides. 


Here we consider the case Pp = ‘4 for all B 
atoms, which are equivalent (the array of A atoms 
in the coordination sphere of B is not considered). 
All the Pa, are then 3, the compound having the com- 
position AB;. We have seen that Pp in space can- 
not be less than 1 for atoms with Pp = +4 orl ina 
layer, so all the B atoms in each layer must have 
P= 1. This corresponds to linking of hexagons 
of B atoms, each of which has an A atom at the cen- 
ter; the links are via vertices and via edges nothav- 
ing a vertex in common (Fig.4). Linking as in b 
and d, or combination of b or d with the other types, 
is impossible because there would be B atoms with 
Pp Seer ig, 5b, d). For the same reason, the hex- 
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Fig. 5. Pp can take only the value 1 at points marked (1), The 
vertices of the hexagons are B atoms, an A atom being at the 
center of each. 


Fig. 6. Possible linking of hexagons in a layer when Pg in space 
is!4. The vertices of the hexagons are B atoms, an A atom be- 
ing at the center of each, 


agons cannot link via vertices along a line through 
an edge joining them (Fig. 5a). This figure shows 
that linking along such a line occurs with a gap equal 
to the length of an edge. Such sequences do not oc- 
cur if the lines are not parallel (Fig. 5c), so link- 
ing along lines that are not parallel is impossible. 
Figure 6 shows three possible types of link. Layers 
of a different type may be obtained by taking differ- 
ent combinations of links of types m and n. 


A space lattice can consist only of a layer of 
the same type placed on a given layer, and this can 
be done in only one way. Figure 7 shows possible 
arrays (three types only) in hexagonal close pack- 
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Fig. 7. Possible structure types for 
hexagonal close packing. The ver- 
tices of the hexagons are B atoms. 
O, A atom in middle layer;x, A 
atom in layer above or below. 


ing; the other types may be derived from layers 
with a different type of link. 

Mg,Cd crystallizes in type a [1], as do some 
other compounds, while type b occurs in the y 
phase of Cu—Al [2]. No representatives of the 
other types have been found in the literature; they 
should have orthorhombic unit cells with a = 2apex; 
C=Chex, b= Ginex’ 33 3a Hones 4ahexV 3, and so 
on. For €-Cu3sSn [3,4], we have a = 2ahex,c¢ 
= Chex, and b = 8ahexV3, etc.; the space group is 
Cmcm, and two models may be proposed for the 
structure (Fig. 9).4 

Figure 8 shows the layer stacking for cubic 
packing; Cu3Au is of type a [1], as are some other 
compounds. Types b and c have the A atoms inthe 


Fig. 8. Possible structure types for 
cubic close packing. x, A at 0;O, 
Avat 1410 Ac at "4p EFA ate 2A 
at heswAat ok: 


fourth layer above the points 1, those in the fifth 
above points 2, and those in the sixth above points 
3. Only three types are shown; the others may be 
derived from layers with a different type of link. 
All structures, apart from those of type a, consist 
of six layers along the 3 axis. TiAl, is a represent- 
ative of type b [1]; ZrAl; is also a representative of 
this series of structure types [1]. All of these 


*The structure of ¢ -Cu3Sn has been determined [Z. Metallkunde, 
46, 9, 692 (1955)]; a and ¢ correspond to [3,4], while b = 

= 10ahexv3, The principle is that of Fig. 8a. The same workers 
have discovered several new compounds that fit the above treat- 


ment, which have cubic close packing and superlattices. (Note 
added in proof.) 
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Fig. 9. Possible structure types 
for e-Cu3Sn. O, Sn in middle 
layers; x , Sn in layer above or 
below. The vertices of the hex- 
agons are Cu atoms, 


types arise from the possibility of placing the cube — 
octahedra of the next layer on the square tops of 
cube —octahedra in the previous layer, or onsquare 
sections of octahedra [5]. 


The arrays in multilayer structures may be 
found by superimposing layers of types a,b,c, etc. 
(Fig.6) as ind andi. For instance, the TiNis struc- 
ture type [1] occurs by superposition of type a 
layers in the order ABAC. 

We are indebted to Academician N. V. Belov 
and to Professor G.B. Bokii for valuable comments 
and advice. 
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The crystal structure of zoisite Ca,Al3[SiyO,][SiO,]O(OH) is completely ana- 
lyzed through a direct determination of the structure-amplitude signs by 


statistical equations. 


Zoisite, Ca,AlsSizOy.(OH), is the third mineral 
of the epidote group [1,2] which has been analyzed 
in the Institute of Crystallography. It is related to 
the first of these, ilvaite [CaFe’** Fe,° Si,O,(OH)]J, 
by the common Fedorov group Déf, = Pnma, and to 
the monoclinic epidote Ca,Al, FeSi30;. (OH) by a 
similar formula, Al replacing Fe. Characteristic 
of the three minerals is one short period (5.5 to 
5.9 kX), along which passes a screw axis 2,, and 
perpendicular to whichisa mirror plane. Along this 
direction in both ilvaite and epidote stretch columns 
of Fe or Al octahedra, linked by common edges, the 
atoms being disposed in only four levels; a repeti- 
tion of this may be expected in zoisite, also. In 
ilvaite and epidote, despite the orthosilicate formu- 
las, one finds the radicals [Si,O,;], predicted for 
minerals of the epidote group by Sobolev [38]. Where- 
as, however, in ilvaite all the Si are in diortho- 
groups [Si,O7], in epidote one orthosilicate group 
[SiO,] corresponds to each [Si,O,]. 

Considering zoisite as a microtwinned epidote, 
Ito [4] calculated the coordinates of the atoms in 
his first epidote structure, obtained by trial and 
error. The failings of the latter were critically 
analyzed in [2]. In order to justify his coordinates 
of zoisite, Ito confined himself to a table giving the 
similarity of the intensities of 50 arbitrarily 
selected hkl reflections, among which 16 were not 
observed; the corresponding coordinates led to un- 
acceptable distances. Some Si—O distances proved 
to be 1.45, 1.90 to 1.97, and even 2.10 kX: the 
lengths of the O—O edges in the Si tetrahedra were 
2.39 and 3.26 kX, etc. (cf. [2]). 


A method of separating out the key group of signs for 
this determination is developed in detail. 


Experimental X-ray Data from Zoisite, 
and Its Preliminary Analysis 


The specimens for examination consisted of 
crystals from the Khorutanii region, obtained from 
the Mineralogical Museum of the Academy of Sci- 
ences of the USSR. The crystals had a well-de- 
veloped zone of elongation, along which their di- 
mensions were 1.5 to 2 mm (0.2-0.3 mm across). 
The observation of transparent, colorless crystals 
made it possible to neglect the Fe content and re- 
gard the formula of zoisite as Ca,A13Si3;04. (OH). 


Unit Cell. Fedorov Group. The pa- 
rameters of the orthorhombic cell (from rotation x- 
ray photographs and from pinacoids far out on zero 
layers of oscillation patterns) were a = 16.20 kX, b 
= 5.50 kX, c = 10.14 kX, in agreement with earlier 
data [5]. From the specific gravity of the colorless 
(i.e., iron-free) crystals (3.25) and the "molecular 
weight" (454.3), we obtain the number of "molecules" 
(formula units) in the cell as 4 (3.90). 

We had available oscillation patterns of the 
zero and three-layer lines from rotation around b 
using Cu and Mo radiation, and the zero lines from 
oscillation patterns around a and c in Cu radiation 
only. 

From the Laue class and the extinctions ob- 
served, we arrived at x-ray group No. 17 [6], 
mmm Pn-a, which includes the two Fedorov groups 
Cy = Pn2,a and Dif = Pnma. 

As indicated, zoisite is a member of a family 
of minerals characterized by a mirror plane per- 
pendicular to the shortest period. Further evid- 
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ence, given subsequently, leads us, with the earli- 
er authors, to the holohedral group Dif, = Pnma. 
The general positions are eightfold and the 
special positions fourfold, divided into the follow- 
ing groups: two groups of four at centers of sym- 
metry (at levels of 0 and b/2) without parameters, 
and positions in the mirror planes (at heights yy, 
and 3b/4). In the unit cell we must put eight atoms 
of Ca, 12 Al, 12 Si, 480, 4 OH. There is no 
doubt but that the Al and Si atoms and OH groups 


fall in groups of four in particular special positions. 


The extremely strong (040) reflection, withthe 
_ short axis b = 5.50 kX, indicates that all the atoms 
are located in only four levels, as in the case of 
epidote and ilvaite, and also in Al,SiO;, BaSQ,, 
Mg,SiO,,Sb253, and a number of others. The "sten- 
cil" nature of the arrangement of the atoms in four 
levels along the short b axis (on account of the wide 
"forbidden bands" due to the large ionic radius of 
oxygen) has been repeatedly mentioned, for example 
in [1,2], where it was shown that this is an argu- 
ment in favor of the holohedral group De 5 TMS) aie 
rangement of the atoms in only four levels also fol- 
lowed from the first "lateral" xy Patterson synth- 
esis. 

Intensities of Reflections and 
Structure Amplitudes. On the x-ray dif- 
fraction patterns, the intensities were estimated 
visually by comparison with 17 blackening marks 
[7], which differed by a factor of V2 and complete- 
ly covered the range of the film emulsion. From 
the relative intensities, taking account of polariza- 
tion and kinematical factors, values of the struc- 
ture-amplitude moduli were obtained in [8]. On di- 
viding by the sum of the atomic factors for corre- 
sponding sin¥/A, all the |Fy;,7| were converted, 
apart from an unknown multiplier, into unitary 
structure amplitudes lFnx 7! ,; required for the sta- 
tistical method of structural analysis [9]. 

Comparison of the experimental material ob- 
tained from the Cu and Mo radiation patterns indi- 
cates the considerable quantitative and qualitative 
advantage of the latter (Table 1). 

The number of h0/ and hill reflections in the 
Mo patterns was almost three times greater than 
for Cu; the Mo lacked only nine of the reflections 
available in the Cu patterns: (200), (201), (002), 
(102), (202), (112), (311), (812), and (411). 

The advantage of the Mo patterns also appears 
in the fact (important for the statistical method) 
that, among the reflections with large indices, those 
for which the unitary amplitudes are large occur 
relatively often. The distribution of approximately 
a third of the h0l and hil reflections with the larg- 
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TABLE 1 
so ang Nay TES 


No. of reflections 
Radiation 


on patterns 


Total differ - 
ent 


TABLE 2 


Distribution of 
reflections with 
greatest unitary 
amplitudes 
hil 


Radiation 


31 
92 


Cuky 
MoK 


est unitary amplitudes in the Cu and Mo patterns is 
given in Table 2. 

The number 137 includes two reflections, the 
(16.0.7) and (17.0.6) which, although (sin#)/A 
< Ay A Cu, ! are absent from the Cu pattern, prob- 
ably due to strong absorption. 


Statistical Method of Determining the 
Signs of the Structure Amplitudes 


Determination of the structure-amplitude signs 
by direct methods in earlier investigations [10,2, 
11] was divided into two stages. In the first, which 
(as it seemed to the authors) made subsequent steps 
more reliable, the Harker—Kasper inequalities [9] 
were used to establish the signs of a small number 
of maximum unitary amplitudes Fak . Starting 
from this "key" group of signs, the signs of other 
unitary amplitudes were found successively at the 
second stage by means of the statistical equation of 
Zachariasen [10] (for centrosymmetric crystals) 


Srat = S (Spyig ty “Satny, Ree 444) (1) 
or 


Sy = S (Sx, -Su+K;)- 


The keys Sx; and SHAK} found from the inequali- 
ties were supplemented as work proceeded by the 
signs of other amplitudes with large unitary values. 

The method of inequalities demands the reduc- 
tion of the experimental amplitudes to absolute 
values; although this can be effected either experi- 


1Strictly speaking, the construction of the cassette made it possible 
to record reflections with sin 9/A inside the following limits: 

0.127 <= sin3/d yo X 1.360 for Mo radiation and 0.059 = sin 9/ACy 
= 0.626 for Cu radiation. 
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mentally or by calculation [12,13], it involves a 
ereat deal of work. Moretroublesome are difficul- 
ties of principle, since a number of factors in 
structural analysis (absorption, temperature, ex- 
tinction) are hard to evaluate and, furthermore, it 
is general practice to use atomic factors calculated 
for free atoms. The main failing of the method of 
inequalities is the fact that it is only suitable for 
large unitary amplitudes,? and the number of these 
is extremely small when there are large numbers 
of atoms in the cell (no heavy atom among them). 
This justified an attempt to eliminate the stage 
with the inequalities and to establish the key group 
in another way. One of such Mare [15], discussed 
later in detail for the group Den = Pnma ,consists 
of constructing a chain of signs (for example Sy, 


= SH, = —SH3) etc.) from individual equalities be- 
tween SH, and Sy, (in the present case, SH, = SH): 
SH, = —Su;,)- The individual equalities may be ob- 


tained statistically, using, in conjunction with rela- 
tion (1), those relationships between the signs of 

the amplitudes Fy 7» Faxi> Fhki> Fhki which char- 
acterize the given Fedorov group. Each chain in- 
cludes Spx] with indices of a definite type. As shown 
later, for a correct choice of three chains, in gen- 
eral one may arbitrarily assign one sign to each, 
and so fix all the signs in the chains. The signs ob- 
tained serve as key signs for finding the remaining 
Shk] With merely the help of relation (1). 


Comparison of Structure-Amplitude 
Signs. Having selected a pair of reflections from 
among those strongest in unitary amplitude® (''com- 
parison" reflections), we match these with such 
pairs of ("determining") reflections that, as a re- 
sult of the addition or subtraction of corresponding 
indices (independently of the two others) of the de- 
termining pair, the indices of the two comparison 
reflections are ox. ained, after which we write out 
the two correspond..g Zachariasen equalities. By 
means of relations characteristic of each Fedorov 
group, it is easy to decide whether the signs of the 
comparison SO OS are the same or opposite. 
For the group D, he Pnma, such relations are [16]: 


h-+l 


a) for yanking zie POTN SAD Seti ee ihe 
b) for : i ie ma ih = = ae a yay = OP ae 
e) for ; 3 . ais Mani == Fit = — Pine - > bgt == — Fy. 
ae 
d) for ‘ 3 2. y Pina == Pky o> Pint = — Pais == — Pnxi- 
(2) 
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Let us take an example. For the comparison 
reflections 808 and 8.0.20, we find the determining 
pairs 


1) 306 5.014 5) 47.4.6 9.1.14 
2) 29.014 21.0.6 6) 18.1.14 LONE 
3) Bm oe 29.0.14 Dit. taal sI6 
4) 38.0.6 30.0.14 
Then, 
1) Ssos = Ss0e X 5.0.14 = — 3305 * 5.0.14 
S¢.0.20 | S'306 x ‘5.0.14 
2) Sgog = So1.0.6 X S29.0.14 
S'g.0.20 = So1.0.6 X S29.0.14 = — Set.0.6 X S29.0.13 
3) Seog = Saz.05 X Sig0.1a = — Sa7.0.6 X S30.0.14 
S's.0.20 = S'37.0.6 X S59.0.14 
4) Sgog = 538.06 X 950.014 = Sas.0.6 X S36.0.14 


Ss.0.20 = S38.0.6 X S30.0.14 
9) Seog = Siri X Sosa = 
Se.o.20 = Siza6 X So.1.44 

6) Seog = Sis.aaa X Sh0.1.6 
Ss.o20 = Sisaaa & St0.1.6 = Sisia X Si0.16 
7) Sgog = So5.114 X Si7.1.6 
S'g.0.20 = So5aaa X Si7a6 = — Sosasa X S17.1.6 


y — 
— Sia1.6 X So.1.74 


With a "probability" of 0.72, the signs of the 
amplitudes of 808 and 8.0.20 are different. 

Relation (1) is statistical, and in comparing the 
signs it is desirable to have the largest possible 
numbers of pairs for determining whether the signs 
of the comparison pair are like or unlike, remem- 
bering, of course, that any application of (1) deter- 
mines the signs of the H amplitudes with a greater 
probability, the greater the unitary amplitudes of 
Kj and H + Kj. 

Method of Selecting the Comparison 
and Determining Pairs. For an arbitrary 
selection of the comparison pairs with positive in- 
dices h,k,/, and h)k,l,, the indices of the determin- 
ing pairs, in the case of the Doh group, are found 
from the conditions: 


h=h—h;, h=h—k, L=i—-L, 


hg =hy-E kh, “ke = hia k;, b= bl, (3) 


where hj, hj; kj, kj, Li, Lj (the indices of the deter- 
mining pairs) may be positive and negative. It was 
shown earlier [15] that the greatest number of de- 
termining pairs is obtained on comparing ampli- 
tudes with two equal analogous indices. For one of 
such indices we chose k, taking k, = kj = ky, and 


If we use inequalities including a small number (3 or 4) of ampli- 
tudes [9, 14]. 


3we shall everywhere replace the unitary amplitudes by quantities 
proportional to them, 
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then the k indices of the determining reflections will 


be any pair of numbers corresponding to the condi- 
tions kj = kj = kj. The two triplets of indices of the 
determining pair each correspond in reciprocal 
space to one lattice point in the planes y * = kj and 
yt= kj, and the triplets of indices of the compari- 
son pair to points in the plane y* = ky. If in re- 
ciprocal space we separate out all possible pairs of 
planes (y* = kj, y* = kj) with the condition kj — kj 

= ky, then each such pair gives determining pairs 
(or their lattice points) for the comparison pair, 
the points of which lie in the plane y* = ky. The 
rules for choosing the pairs are found by consider- 
ing all possible values of (hj,/j)k; and (hj, Jj ij as 
functions of the values of (hy 2 4)i, and (hg, 2 2)ko- We 
distinguish three cases: 

1) Gihy, “who; 9t, File (see Eable.3, la, [b,21¢;, 
It, Bs Flo, ZS, Zaye 

2) eA Tipe bosel p= bo Ale, Jide le, 2¢,.20;, 2e); 

Bh, ~ hy; 1, # ly (8a, 4a, 5a, 3b,-4b, 5b, 
3g, 4g, 5g, 3h, 4h, 5h); 

3) hy = ha; 2, = l_ (degeneracy condition). 

The results of this consideration are shown in 
Table 3, where for each case we give rules for the 
geometric choice of pairs in three superposed 
planes (y* =k), y* =kj, y* = kj) of the reciprocal 
lattice. 

We see from Table 3 that the maximum num- 
ber of determining pairs is obtained on comparison 
of reflections with two equal analogous indices. 

For h, = by, 1,= 1, =1, (A), the indices of the 
determining pairs, according to condition (8), will 
be 


hy -+ hg hy, —he 
2 ae 2 


h; = 


= i? & Li a l; — LF, 
and the indices of the comparison pair are formed 
from those of the determining pair by the formulas 
(hy, ko, fy) —— (h’, hj, l;) ae (h", k;, L), 
(ho, Kg, lo) = (h', ki, di) a (h", ee lj). 
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For hy = hy = hy, ly Z lo (B), 


Ly —_l, sats i 


1 
Rdg eed =it*_y, Le 


and, clearly, 


(Ao Ko, 1) = (hi, hi, LU’) — (h;, Fes ES 
(No, Ko, le) = (hi, hi, l') — (hj, Bj, €"). 


If as determining reflections we could have 
chosen those with any unitary amplitudes, then, on 
the basis of relations (2), the same number of de- 
termining pairs would support both the identity and 
the opposition of the signs for the comparison pair 
of reflections (for example, Shoko 4 = Shokpl » for hj 
+1" = 2n, kj = 2n, and Shokol , = —Shokpl. for byiti® 
= Wyn ceil. kj = 2n). 

Our problem, however, is the selection of de- 
termining pairs with maximum unitary amplitudes, 
the statistics of which lead to one definite sign. The 
scheme for the geometric choice of such pairs is 
shown in Fig. la (scheme A) and Fig. 1b (scheme 
B). Here, in the square (h,/), networks, lattice 
points with large unitary amplitudes are separated 
out, and by simple superposition the three networks 
(y* =k), y* = kj, y* = kj) are placed together in 
one figure. The H, (hj, ky,/9), Hy (hy, Kp,/ 9) and, cor- 
respondingly, Ly (My,ky,/4), Ly(Qo, kp,22) are the 
points of the comparison pair; Hj(h', kj,/;), as also 
L jj, kj,/"'), is one of the points of the determining 
pair; each of the points Hj", kj, Lj) and H jh", kj, 
Lk), justas Lj(hj, kj,/") and Lk kj, 1"), is the 
second in the determining pair [use of the points 
Hk(Lx) in the determining pairs corresponds to 
choice of reflections with negative indices 1, (hx)]. 
In the schemes, AB belongs to the kj network, and 
A,B, to the kj network; we may, however, take AB 
for the kj and A,B, for the kj, which corresponds to 


the condition kp = kj — kj. 


“In Table 3 and in the figures, h" and 1" are indicated by 0 /2, 
and h' and 1' by h, + (a/ 2)[1, + (o/ 2)]. 


Fig. 1. Schemes for selecting determining pairs. 


a) For condition hy > hg, 1, 


= 15. b) For condition hy = hy, 2, > Lo. 


ORYSTALS TRUGCLURE OF ZOIsTLE 


Any of the so-chosen combinations HjH;, HjHx, 
and LiLj, LiL, gives a determining pair for the 
comparison pairs H;H, or L,Iy. It is not hard to 
see that selection according to schemes A and B 
gives all possible determining pairs for the given 
pair of planes y* = kj, y* = kj. In order to obtain 
the maximum number of determining pairs, choice 
is made according to these schemes for all possible 
pairs of planes satisfying the condition ky = kj — Kj 

COMA ERUGWIOM Ot Cingiins emel Cmoiec 
of Initial Signs. Pairs for which the identity 
or opposition of the signs is substantiated statisti- 
cally are connected in a chain of signs through re- 
flections common to different pairs. Reflections 
with indices belonging to only one of the eight groups 
ID) 0) = Ail, IRS il, Yh Se PA) sel An, AAV BY Any. 
Pin. Zio sile 49) Bin, 2ieril, Zine Gy) Ansel, Bin, BinPile 
@)) ZapsP ik, Bia seil, Bins ) An, Ail sil, An Pils B) en Pal, 
2n+1, 2n +1, can enter into each chain. 

Out of the set available, we may assign signs 
arbitrarily to three reflections [10,17,18] and by 
this fix the origin of the cell in one of the eight 
translationally nonidentical centers of symmetry of 
the unit cell. Let us,in fact,denote the coordinates 
of the centers of symmetry in the primitive cell by 
(E14, &2, &3), where € can only take the values 0 and 
'/; then the structure factor 


7/9 


N/2 
Fry = 2 >, f; cos 2x (he; + ky; + lz;) 


j=1 


on transition to an origin in any other center of 
symmetry, transforms into 


Frat = cos 20 (he, + key + les) Fret, 


where, for cos 27 (he, +ke, + 1&3) either +1 or —1 
are possible values, i.e., the sign of F};,7 changes 
or remains the same, depending on the choice of 
coordinates for the new origin and the parity of the 
indices h, k, l. 

For reflections of the first group, the sign of 
the amplitude does not depend on the choice of ori- 
gin at one center of symmetry or another, and is 
determined exclusively by the structure. Signs may 
be assigned to the amplitudes of the remaining 
groups quite arbitrarily. In fact, since the signs of 
the amplitudes of the second group with even k and 
1 and odd h change on transferring the origin to any 
of the four centers with ¢, = 4, on assigning the 
sign of a structure amplitude of this type we fix the 
value of €;. By assigning the sign of an amplitude 
from the third or fourth groups, we fix €3 and, 
correspondingly, €,. A fixed sign for an amplitude 


137 


of the fifth group determines whether the sum &, 
+€. is integral or fractional. In order to fix the ori- 
gin in one of the eight centers of symmetry, we 
must assign the signs of three amplitudes of differ- 
ent groups in such a way that the sign of one of them 
is not a consequence of the assignment of the other 
two. 

An origin in one of the four centers of sym- 
metry of a primitive projection of the cell on, for 
example, the xz plane, is fixed by the assignmentof 
signs for two structure amplitudes Fj), of different 
groups (excluding the first group). 

Selection of the Key Group of Amp- 
litudes and Determination of Si nisi bx— 
perimental values of Fpo, were reduced to quanti- 
ties proportional to Fro , allowing for the tempera- 
ture factor exp [—-B(sind/))*] (10,19, 22]. Accord- 
ing to [13], for the comparatively refractory sili- 
cates (beginning of the endothermic process in zois- 
ite is at 1000°C, compare muscovite, 860°C), the 
factor B = 1 [kX]*. With some reserve, we took B 
= 0. TkKX( 

Initially we selected 50% of the reflections with 
the greatest Fyo,. Choice of the determining pairs 
from this group, despite their relatively large num- 
ber, led to unsatisfactory statistics, and the key 
group was reduced to 30% (137 reflections in all). 
Since this number was small for obtaining a suffi- 
cient number of determining pairs, to the key group 
were added 30% of the strongest Fy] (92 reflec- 
tions), which were used only for forming the deter- 
mining pairs. All the h0l reflections taken into the 
key group, and also the auxiliary h1ll, were setout 
in two square networks (simplified h0/ and hil 
planes of the reciprocal lattice), and with the help 
of these the selection of the comparison and deter- 
mining pairs was made by schemes A and B. The 
statistics proved satisfactory, and it was possible 
to compare approximately 60 pairs of reflections. 
For three or four determining pairs, the difference 
or similarity of the signs was established "reliably," 
and for five or more with a probability not lower 
than 0.80. Later, these pairs were united in the 
chain of signs shown in Fig. 2. 

For the assignment of signs it was natural to 
choose chains la and 2a, which included the largest 
number of reflections, obtaining the signs in terms 
of the reflections (15.0.9) and (25.0.4), which (in- 
side the chain) had the largest number of linkages. 


°As shown above, for comparing the signs of the amplitudes Fyo), 
it was sufficient to make a selection of determining pairs in the 
plane h0Z and afterwards in the plane hil. For comparing the 
signs of amplitudes Fh17, the comparison pairs were taken in the 
plane hl/ and the determining pairs one by one from h0l and hil. 
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Z \ 4 
0.20 — S901 e S 29.0.2 
\ 
— $ 20.0.4 
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e wk 04 
“eS 
— $38.08 — S30.9.10 S3z.0.16 $1004 
LE aS 
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NEY 
A) — $8.0.4 == §$48.0.16 TSet.0 1 
Ye NN ~< aa NK 
$3.0.0 S$.u.20 —— $10.0.16 bo — Sau == S410 —$6.0 
S 
——§8.0.22 $410.4 
@) S14.0-3 $§10.6.15=—S10.0.21 
Yl 
a seo10=—St6e ¢ $18.0.112—S18.0.17 
SS 
Sidr J $ 20.0.0 =— §20.0.13 


e §16.0.7=— $160.11 


Fig. 2. Chains of signs for h0l reflections. 


Putting Sj5 9.9 =—, and Ss 94 = +, we obtained the 
signs of 21 amplitudes from chains la and 2a. With 
the help of these signs and relation (1), we succes- 
sively established signs for (1.0.12), (5.0.14), 
(6.0.10), (8.0.10), (29.0.2), (88.0.2), (808), and 
(4.0.14), respectively, entering the chains 2b, 2d, 
3b, 2c, 2e, 3d, 3a, and 3c. This extended the num- 
ber of determined amplitudes from 34 to a total of 
55. In the unused chains 4a, 4b, 4c, 4d, and 4e we 
were unable to determine the sign of a single ampli- 
tude by means of these 55. 

The 55 amplitudes whose signs were fixed fair- 
ly convincingly served to provide a statistical de- 
termination of the signs of the remaining amplitudes 
selected for the key group. First we determined the 
signs of those amplitudes H for which the greatest 
number of pairs (Kj, H + Kj) could be formed from 
those already determined. Selection of the pairs 
was made (in accordance with Vainshtein [2]) by 
means of movable and stationary networks. After 
fixing the signs of the whole key group, the sign of 
each of the 137 was further checked statistically 
from pairs taken from the other 136.° Since, for a 
large number of determining pairs, the signs of 
some of the amplitudes were not confirmed, the cor- 
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responding reflections were excluded from the key 
group and not used in establishing the signs of all 
the amplitudes. In the last stage a sign was con- 
sidered established if, for 20 to 50 determining 
pairs, the probability was no lower than 0.66. 1 
Construction and Analysis of Electron- 
Density Projections 


The First xz Projection of Elec- 
tron Density. In constructing the first elec- 
tron-density projection, 278 amplitudes Fpo] were 
used, the signs being determined by the statistical 
method. Identification of the peaks in the numeri- 
cal table obtained was easy on the basis of the struc- 
tural link between zoisite and epidote. If we disre- 
gard the iron present in epidote (approximately one 
iron atom for two of aluminum), the zoisite (which 
does not usually contain iron) is a polymorphic mo- 
dification of epidote with a similar structure, since 
two dimensions of the unit cell almost coincide, 
while between the third ones there is a simple re- 
lation 2a sinB = a'. 


Epidote Zoisite 
b = 5.62 kX, b' = 5.50 kx, 
c= 10.23 kx, c'= 10.14 kx, 
a= 8.89 kx, a'= 16.26 kx 
B= 115°24', (2asin115°24'=16.06 kX), 


Zea 
Sey bye ee VSB 
Doh gee 


The orthorhombic zoisite is a microtwin of 
epidote with twinning plane (100). According to 
analysis of the extinctions (p. __), the diagonal 
glide plane n serves as this twinning plane (100). 

In the monoclinic cell of epidote (Fig.3), one 
set of four centers of symmetry coincides with the 
centers of the Al octahedra, and the other set of 
four inverts the [Si,O7] groups at levels b/4 and 
3b/4 one into the other (in Fig. 3, pairs of such 
groups are located in the middle of the a edge). De- 
spite the doubling by the twinning plane, both in 
volume and in material content, the orthorhombic 
cell of zoisite is nevertheless primitive, and there- 
fore half of the "epidote" centers of symmetry 
vanish. Consideration of the epidote cell shows that 
the Al octahedra (at 0 and c/2) already lie almost 
exactly in the plane which would correspond to a 
glide plane passing through them and through the c 
axis. In zoisite, this plane becomes the actual di- 


2 
f= Pit. 
Cth = m 


®Not strictly from 136, since not all the signs of the 137 were found 
statistically, 


"In determining the key group we had no more than 10 to 12 pairs; 


this also was responsible for the rigorous demands indicated above 
for this group. 
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Fig. 3. Structure of epidote in coordination 
polyhedra. Projection on the xz plane. 


agonal glide plane, and then the centers of sym- 
metry lying in it should vanish, since, in the group 
Df all the centers of symmetry lie midway between 
planes. 

Thus, in zoisite, the Al atoms do not lie in 
centers of symmetry, and the latter are retained 
only between pairs of diorthogroups [Si,O,] in the 
middle of the a@ edge, to which we must also refer 
the origin of the epidote cell, if we are to calculate 
the coordinates of zoisite from the coordinates of 
epidote. The translation reduces all the x coordi- 
nates of the epidote by a/2, and on taking account 
of the doubling of the a axis in zoisite, the resultis 
further divided in half. Exchanging the oblique- 


angled cell of epidote for a rectangular cell, accord- 


ing to Fig. 4, reduces all z by 


“ sin 25°24’ (« = =} = 0.38 (« — s): 

Thus, the formulas for the transformation to 

the rectangular cell (where the new coordinates of 
epidote should, to fair accuracy, be the coordinates 


of zoisite) wilk be: 


SS 

I 

nN 

| 

= 

w 

co 
SS 

8 

| 
| = 
eee? 


The coordinates obtained from these formulas 
for zoisite (Table 4) may be compared with the fig- 
ures which emerge from the statistical synthesis. 
Comparison shows that these are in fair agreement. 
Apart from this, however, the heights of the peaks 
in the synthesis did not correspond to the atomic 
numbers. Thus, the peak of the twinned Al was 
more than three times higher than that of the ordi- 


Fig. 4. Scheme of relationships between the monoclinic cell 
of epidote and the orthorhombic cell of zoisite. 


TABLE 4. Coordinates of Atoms in Zoisite 
from Second Synthesis 


No. of atoms Coordinates 
atems in cell 100 x/a | 100 2/c 
Al, 4 25.0 a= 19:05 
Alyy 4 25.0 19.05 
Alla 4 10.7 30.0 
Gar 4 BI 43.6 
(eee 4 45.4 44.2 
Si; 4 Soe Ore 
Sirs 4 44.2 28.1 
Site 4 16.2 43.5 
Oy 8 12.7 15.0 
On 8 10.1 ont 
oT 8 Bio |) 2dse 
a A R2NSM 80 
Or A 22.3 | 30.4 
On A 26,850 m tol 
Oo 4 = ("2 16.1 
Ovryt 4 0 29.6 
4 43.9 43.8 
On 4 26.85 Leh 


nary Al. Whereas Cay had height 451, the height of 
Cajj was 200. For two of the Si atoms, the heights 
were over 250, while for the third the height was 
below 200. The single Oy yy and Oyyy7 peaks could 
not be distinguished from "ghost" peaks of similar 
height. 

The reason for this inadequate revelation of the 
light atoms is clear. Owing to the rapid fall in 
atomic factor, these only introduce a substantial 
contribution into amplitudes from planes with large 
interplanar spacings and small hkl. It is the de- 
termination of the signs for just these amplitudes 
which presents greatest difficulty to the statistical 
method. 
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Fig. 5. Electron-density projection of the cell of zoisite on the xz plane 
(quarter cell). 


The Second xz Projection of Elec- 
tron Density. From the coordinates of the 
statistical synthesis we calculated the amplitudes 
whose signs were not given by the statistical me- 
thod, i.e., 407 —278 = 129 amplitudes. As indi- 
cated, a significant part of them comprised the 
amplitudes of the fundamental tones. Twenty-two 
structure amplitudes (all with large indices) were 
found to be insignificant and changed sign for small 
changes in coordinates. These were discarded in 
constructing the second synthesis, presented in 
contours in Fig. 5. The peaks of the two kinds of 
Ca became equal in height. The ratio between the 
heights of the peaks corresponding to the twinned 
Alj-]] and single Alyqj, and that between the peak of 
the twinned Alj-yj (815) with "atomic number" 2° 13 
= 26 and the peak of Ca (600) with Z = 20, were 
both well preserved. The three sorts of Siz, Siyy, 
and Siyyy have respective heights 512, 425, and 504. 
The peaks of the twinned O with heights from 362 
to 458 decompose into sharp and blurred peaks. The 
first correspond to pairs of O which are connected 
by a mirror plane and exactly overlap each other, 
and the second to pairs of O which are mutually 
superposed not from the requirements of symmetry 
but only from the principle of closest packing of 
atoms in the octahedra around the Al. The single 
OyqI and Oyqy, which evaded the first synthesis, 
are fully revealed. All the single O have maxima 
from 217 to 274. Nota single side peak ("ghosts") 
exceeds 160 (background about 100). The coordi- 
nates emerging from the second synthesis appear 
in Table 4. 

From these almost-final coordinates we calcu- 
lated the amplitudes which are compared with ex- 


perimental values in Figs. 6a to 6i, which show 
420 Fyoz, including 407 which, from experimental 
data, were nonzero and 13 zero types with h and 1 
not exceeding 10. The points on the broken line in- 
dicate that the sign of the corresponding amplitude 
was established by the statistical method. 


The coefficient of reliability for h0l reflections 


p= =I Fesch |= | Featel 
SIF 


esch | 
equals 0.22 allowing for zero amplitudes and 0.21 
without these. 


Construction and Analysis of 
Conditional Projections 


Although there was no doubt as to the y coordi- 
nates of the atoms in the zoisite structure, since 
these corresponded to the heights of the analogous 
atoms in epidote, we nevertheless verified the hori- 
zontal and vertical coordinates by the method of 
conditional or generalized projections [20]. The 

1 
corresponding \ (xyz)sin 2 ydy function for a struc- 
0 
ture similar to zoisite with atoms at only four 
levels, as we can easily see, gives a projection in 
which atoms at the levels 0 and b/2 vanish, where- 
as atoms at levels b/4 and 3b/4 are fixed as peaks 
for the first, and troughs or pits for the second. 
The function 


1 [o9) 
b \ osin 2x ydy = — 2 | > > Fy sin 26 (hx + r) | 


0 h, l=—o0 
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Fig. 6. a-i) Graphs showing the correspondence of calculated and experimental 
values of Fo] (continuous and broken lines, respectively). 


takes the following form for calculation: 


ce foo} 
4 d . ‘ ; 
eet (2 »; Fiyo sin 2n ha + 2 pa Fou sin 2nlz 
h=1 l=1 


+4 >) by Fry Sin 2x hz cos 2alz 


h, l=1 


The single primes indicate summations refer- 
ring to Fhij] with h+/ = 2n, and the double primes 
summations referring to Fphj7 withh+J/ = 2n+1. 

The Key Group of Structure Ampli- 
tudes hil. Weused the method of comparing 
signs and constructing chains which was employed 
earlier for the Fpo], and in the final stage,relation 
(1). The problem was made easier by the existence 
of around 300 amplitudes Fyg; with signs deter- 
mined earlier and confirmed by calculation. 

At the first stage (in comparing signs), weused 
172 of the greatest amplitudes Fyo] with reliable 
signs and 92 Fyi/ which took part in comparing the 
signs of the Fyo7 (in the determining pairs). The 
Fo, now entered exclusively into the determining 
pairs, one in each. The comparison and determin- 
ing pairs were chosen by the same schemes A and 
B. Of the comparison pairs, only those which had 
received the best statistical substantiation were 
chosen. The chains of signs are given in Fig. 7. 
The closing of the links enabled us to detect errone- 
ous linkages, for example in chains 4 and 5, where 
they are shown by broken equality signs. 

The x and z coordinates of the origin were fixed 
by the earlier assignment of signs for two ampli- 
tudes. In order to fix the third, the y coordinate, we 


+4 yl >: Fri, cos 2% hx sin 2rlz 


h, l=1 


arbitrarily assign the sign of one of the Fy, ,] ampli- 
tudes. Having a large number of rigidly established 
Fhol signs, and remembering the possibility of un- 
detected erroneous linkages in the Spi] chains, we 
chose as basis not the longest chain, but only the 
closed links of chain 4 with amplitudes (517), (717), 
(13.1.7), (511), and (5.1.17). Setting S,,;7 = —, we 
have Sii7= + Ssit = + S5.4.17 = + and Sy34,7 = —. With 
the aid of the five signs so furnished for the ampli- 
tudes Fp ij and also the 172 Fpo], the signs in the 
remaining chains were successively found, and then 
the signs of the key group; after verifying each de- 
termined sign by statistics, using all possible pairs 
from the same key group, we proceeded to apply the 
statistical equation (1) to the rest of the amplitudes. 
Altogether, the signs of 218 Fphij amplitudes were 
established out of a total of 277,i.e., 79%. 

Two Conditional Projections. 
conditional projection 


The 


1 
\ op sin 2x ydy 


0 


was synthesized twice: the first time from 218 
Fhij amplitudes with statistical signs, and the se- 
cond with the addition of amplitudes, the signs of 
which were calculated from the coordinates. In the 
second synthesis, we used 273 amplitudes, since, 
out of the total number of 277 experimental Fy17 

~ 0 values, six were excluded for the same reasons 
as in the case of the Fyo1, and two, (111) and (211), 
were added, these being substantial in absolute 
magnitude, but falling in the "dead zone" of the x- 
ray photograph. The improved conditional projec- 
tion is given in contours in Fig. 8. This contains 
all the peaks of the full projection, except the Aly_yj 
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Fig. 7, Chains of signs for h1l reflections. 


Fig. 8. Conditional projection fpsin2mybdy. Negative signs of the function are indicated by broken contours. 
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Fig. 9. a-h) Graphs showing the correspondence of calculated and experimental 
values of F},,7 (continuous and broken lines, respectively). 


TABLE 5. Coordinates of Basic Atoms 
of the Zoisite Structure 


Atoms 100 100 100 
ala y/b zie 
Alvi (8) 25.0 0 19.05 
Alyy (Fe) (4)| 10.7 75 30.0 
Ca; (4) | 7.2 25 43.6 
Cazy (4)| 45.4 23 4422 
Siz (4) | 8.45 25 40-2 
Siz; (4) 41.2 75 28 1 
Sizs; (A) | 16.2 | 25 43.5 
O; (8) 11.8 0} 45.0 
On; (8) 40-1 0 43.4 
On ©} 375 | 0 | %5 
o (4)} 22.3 | 25 30.4 
Oy (4)| 22.3 75 30.4 
On (4) 26.85 75 | D9. 
io: (4)} —0.3 | 25 anis6u3 
ae (4)| —0.3 Tet ecole 
Ox (4) 43.5 75 43.8 
OH (4) 26 .85 25 IoB 
| 


and three sorts of O (I, I, Ill), since, for these 
atoms, the y coordinates equal 0 or ye Within the 
quarter cell the two sorts of Ca occur at different 
heights. The Siz and Si, atoms from a single di- 
orthosilicate group and the Oyx atoms common to 
the two Si in fact occurred in a common mirror 
plane, once in the form of a triplet of peaks and 
once as a triplet of troughs. 

The coordinates of the atoms obtained from the 
conditional projection, after comparison with the 
full projection, were used for calculating the struc- 
ture amplitudes Fy1]. Figures 9a-h present graphs 
comparing Fh1] exp With Fhil calc for 295 Fhil, 
among which 277 are nonzero and 18 are zero with 
i= 40:; 


Description of the Structure of 
Zoisite 


Coordinates of the Basic Atoms. 
Table 5 shows the final coordinates of the basic 


atoms of the zoisite structure. The four groups of 
eight in general positions are each characterized by 
three parameters, and the thirteen groups of four 
in mirror planes by two. The whole structure is de- 
termined by 38 parameters (as against 35 in epidote, 
where, although there are two kinds of Al atoms, 
both are in centers of symmetry,i.e. , without pa- 
rameters). 

We used Vainshtein's method [13] to calculate 
the probable error of the coordinates associated 
with inaccuracy in the intensities. For the case of 
the projections, the corresponding formula is 


V >a 


Ar, = 0.049 at : (4) 


vs 


7 
1.35 
Zi, 


where b is the confidence coefficient of reliability, 
S the area of the projection, and Zj, Z; are atomic 
numbers. In our case (Mo radiation), the series 
was broken off at sin’/A = 1.36, and from the 
graphs given in [13], beff = 2b = 0.4. Thus, the 
probable error in the coordinates introduced by in- 
accuracy in the intensities is Arcg = 0.008kX, 
Ara = 0.012 kX, Arg; = 0.012 kX, Aro = 0.024 kx, 
i.e., is considerably smaller than inaccuracy as- 
sociated with the low resolving power of the 48? di- 
vision of the unit parallelogram. 

The sketches in Table 6 indicate the coordina- 
tion polyhedra for each crystallographically inde- 
pendent Si and Al cation; alongside are indicated the 
cation—anion distances and the length of the oxygen 
edges in the polyhedra. The Si—O distances proved 
to lie within even narrower limits than those of 
epidote (in zoisite, 1.53 to 1.70 kX; in epidote,1.51 
to 1.73 kX). The Al—O and Ca—O distances were 
approximately the same in the two structures. The 
Si-—O-Si "valence angle" of the diorthogroup in 
zoisite, ~161°, is close to the corresponding angle 
in epidote (158°) and quartz (150°). 

Description of the Zoisite Struc- 
ture. Figure 10 shows the xz projection of the 
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TABLE 6. Scheme of Environment of Cations and Interatomic Distances 
in the Structure of Zoisite, in kX 


Sij -tetrahedron 


Siyy -tetrahedron 


Siqy -tetrahedron 


v 

I lm M 
Sil 
a Dn ¥ sil 
0 
Pe Din ie 

Slo— | = 464 Si— IJI]= 1,68 Si 1) =4,70 
Si — VII = 1.53 Si — VIII = 1.44 Si— V = 1.65 
Si— IX = 1.67 Si— IX = 1.63 Si— Vi = 1.67 
I—VII =2,52 II—VIIIl =2.71 IPSS PV = 2273 
I— IX =2.73 Tit IX = 2.72 II—VI = 2.82 
VII— IX =2.54 VIII— [XX =2.77 VV es 
ee 270 WI— Ill, = 2.75 Ii— Il, = 2.75 


Aly jj -octahedron 
Ill 


7 Iv 
I 

Al— 1=2,02 

Al — III = 1.83 

AL—IV = 1.87 

Al— V =1.87 

Al— VI = 1.93 


Al — OH = 1.83 
IV — OH = 2.44 


V— VI = 2.62 

I— VI = 2,82 
J— 1p =2.94 
I—OH =2.77 


— 
III— IV 
iI— V 
Ill — OH 


yesh yy Wee 
y } = 2.58 


\ 2.88 
= 2,80 


structure, in which the oxygen tetrahedra around 
the Si and the octahedra around the two kinds of Al 
are marked out. The Ca atoms are shown by 
circles from which bonds pass to the closest O 
(OH). Figure 11 gives the axonometry of the struc- 
ture. 

The main features of the zoisite structure re- 
peat those characterizing epidote. These include 
first of all vertical "guide rods," columns of Al 
octahedra extending along the y axis and passing 
close to the middle of each quadrant in the xz pro- 
jection. In contrast to epidote, in the zoisite struc- 
ture the axis of the columns does not coincide with 
the 2, axis, and the centers of the Al octahedra do 


Al— I=2.08 
Al— I= 1.92 
Al— IV =1.87 
Al —VIII = 1.79 
IVa "Y= 2,50 
I— VIII = 2.93 
W— Iv=2.73 
I] — VIII = 2.57 
fee” Wen G08 

T— In =2.5 
II — Im = 2.75 


1. Roman figures indicate oxygen atoms 
2, Atoms I-Ip, U—I, WI—IIp are 


connected by mirror plane m. 


not coincide with the centers of symmetry. Along 
the y axis the octahedra of one column are con- 
nected only by mirror planes through a common 
edge of each two successive octahedra. The neigh- 
boring columns of Al octahedra are crystallographi- 
cally connected by diagonal and parallel glide planes 
passing through them in the z and x directions, re- 
spectively. In epidote, the single Al (Fe) octahedra 
intimately adjoin one half of the columns of Alocta- 
hedra on two sides, in accordance with the law of 
the two-fold screw axis, and are not connected with 
the other half of the Al columns. In zoisite the 
single Al octahedra are attached to both columns 
from one side only. This is the sharpest difference 


CRYSTAL STRUCTURE OF ZOISITE 149 


K( 


\ 


fy 
©, 
F( 


Fig. 11. Structure of zoisite in coordination polyhedra. Axonometry. 


between the structures of zoisite and epidote. In 
epidote, the Al atoms are in the centers of sym- 
metry, and among the vertices of each Al octahed- 
ron are either two OH groups (Al;7) or none (Al)). 
In zoisite, each Al is connected with one OH group. 


The Si atoms are of three kinds: Si; and Siy; 
enter into the diorthosilicate group and Siyqq into the 
orthosilicate. All link up with columns of Al octa- 
hedra, so that, in the xz projection, we find char- 
acteristic pentagons made up of two Al octahedra 
and three Si tetrahedra (Fig. 10). 


Despite the orthosilicate ratio Si:O = 1:4 fol- 
lowing from the formula, zoisite (like epidote) is 
nevertheless a silicate of the mixed type with two 
kinds of silicon—oxygen radicals: [SiO,] and [Si,O7]. 
Of the 12 oxygen atoms in the zoisite formula, 11 
take part in the silicon —oxygen radicals; the twelfth 
O and the OH group are not connected with Si and 


are only saturated by the Al and Ca cations. Inpro- 
jection, and still better in axonometry, it is clear 
that there are two kinds of such O atoms: Oyy and 
Ox. According to the balance of valences, the Ox 
belongs to the OH, since the smallest sum of "ca- 
tion forces" converges in this atom. 


The Ca atoms in zoisite are of two kinds. The 
coordination polyhedra for both are heptahedra, 
constituting (as in epidote and ilvaite), the combina- 
tion of a triangular prism with a hemioctahedron. 
The Cayj polyhedra "stand" with the triangular base 
of the prism on the xz plane, while the Cay hepta- 
hedra "lie" on the same plane with the side edge of 
the prism. The distances from the Ca to the sever 
O are from 2.3 to 2.8 kX. For Cay we can find two 
more distant neighbors with distances ~3 kX. The 
Cay has three neighbors at a distance rather 
sreater than 3 kX. We find a similar environment 
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Sum of 
valence forces 


*/4 "4 


2 

13/5 
4/4 2/4 
fe | 21/4 
ig Ar 7/4 1/2 
3/5 13/4 
13/4 


TABLE 7 
es See 
Aly-1I Aljyr | Cay Caqy Siz | Sirr 
O; 3/5 3/6 “18 
On 3/6 2X Js ‘l4 2 
Ont Ws 2/9 2/5 /4 
Ory 2 x 3/6 3/6 
Oy 2x */6 /g 
Ovy 2x %/¢ */s 
Ovit “Is 
Oviit “ls , 
IX ie 
OH BE<o/e */g 


for Ca atoms in almost all those cases (for example, 
in apatite, sphene, epidote, and ilvaite), where the 
Ca occupies positions not in centers of symmetry. 

The expanded formula of zoisite may be written 
as follows: 

CayCayy (Aly_y7)2(Al, Fe) O (OH) [Siz0,] [Si0,] 
= CaAl, (Al, Fe) O (OH) [SizO,] [Si0,]. 

The symbol Alj_y indicates the merging into a 
single group, in zoisite, of the two kinds of Alwhich 
are crystallographically different in epidote: 


CayCayAl,Alyy (Al, Fe) O (OH) [SiOz] [Si0,]. 


Balance of Valences. If, in accord- 
ance with Pauling's second rule we set up the sum 
of the valence forces meeting at each O, we obtain 
Table 7. 

According to the table of Si—O distances 
(Table 6), the too-small sum of forces at Oyq] is 
compensated by a reduction in the Si—Oyy distance 
in the tetrahedra. As regards Ojy, its location 
alongside the OH group at a distance of 2.44 kX sug- 
gests a hydrogen bond between Oyy and Ox(OH), 
whereby the lack of valence forces at the Oyy is al- 
so compensated. 

Certain Physical Properties of 
Zoisite. As regards cleavage, the properties 
of epidote are repeated: perfect along the (001) and 
imperfect along the (100), the apparent transposi- 
tion of the cleavage planes causing some surprise. 
From the structural plan of zoisite and epidote we 
should have expected better cleavage in the (100) 
than the (001) plane. 

In optical sign, zoisite is opposite to epidote. 
The negative character in the latter occurs because 
the Fe and Ca cations, which are energetically the 
most influential on the optical properties ,are linked 
into plane networks perpendicular to the [010] axis. 
In zoisite with not more than 2% Fe,O3, the role of 
the networks of Ca and Al (Fe) ions falls into se- 


cond place, while the first is occupied by the 
columns of Al octahedra, which determine the po- 
sitive sign of the crystal. These columns also pro- 
duce the strongly expressed striation of all the 
facets of the [010] zone. 

Structure of Zoisite According to 
Ito. A considerable inaccuracy (10 to 12%) in the 
coordinates of epidote found by Ito was revealed in 
[2]. We find the same in his coordinates for zoisite 
[4]. Even for the heaviest atoms (Ca), the discrep- 
ancy reaches 0.022 c, i.e., 0.22 kX. For Al (Fe), 
it is 0.0125 a = 0.2 kX, and for Siyyyj we have 0.025c 
= 0.25 kX. The largest are the discrepancies for 
the Oyx atom, which was also inaccurately deter- 
mined in epidote, namely 0.045a, i.e., 0.73 kX. 
For the Oyyz and Oyqy atoms the discrepancy is 
around 0.045c, i.e., 0.46 kX. 

Figure 12 shows the structure of zoisite accord- 
ing to Ito. By comparing with Fig. 10 we see its 
main defects: mutilated Al octahedra and unnatural- 
ly featured diorthogroups, caused by the incorrect 


Fig. 12, Structure of zoisite accord - 
ing to Ito, 
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coordinates of O[x, which for Ito overlaps the Ca 
atom in the xz projection. The Oyy and Oyyqyq over- 
lap in the same way. These defects have a detri- 
mental effect on the interatomic distances in the Si 


tetrahedra and Al octahedra. 


Thus, Siy—Orx 


= 1.90 kx, Sitq7—-Oyy = 1.97 kX, and Siq,—-Oyqreven 


reaches 2.10 kX. In the Siyqy tetrahedron there is al- 


he 


10. 
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30 a vastly shortened distance of Siqy—Oy = 1.45 
«X. The distances between the O atoms vary sub- 
stantially: the O~x—Oy qj] edge equals 2.39 kX, and 
the Ofx—Oyqy edge 3.26 kX. The distances in the Al 
octahedra are a little better, although even in the 
Al (Fe) octahedron the Al (Fe)- Ov distance falls to 
163° KX 
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Data on the crystal structure of indium fluoride trihydrate InF; ° 
It has been shown that the structure is of coordination type 
Two sorts of water molecule are present. 


are given. 
without isolated complexes. 


3H,O 


Their crystallochemical differences are in agreement with the chemical 


behavior of the salt. 


X-ray structural studies of indium compounds 
are of interest from the general chemical point of 
view and because the crystal chemistry of indium 
has been very little studied. The present paper is 
concerned with the determination of the structure 
type of indium fluoride, InF;- 3H,O. It is neces- 
sary to solve the basic problem of whether the 
structure is of the coordination type or whether iso- 
lated-complexes of the type InX, exist in crystal- 
line state (X= F or H,O). The latter appeared to 
be very probable, because, in most investigated 
structures of this type, indium has a coordination 
number of six and an octahedron as its coordination 
polyhedron. 

The second problem to which a structural study 
could be applied is an explanation of the role of the 
water of crystallization. 

I. Indium fluoride was synthesized by the me- 
thod of Deichman and Tananaev [1]. The crystals 
were glistening colorless tetrahedral needles; face 
development, however, was very irregular. There 
were vicinals on almost all faces, which made pre- 
cise geometric measurements very difficult (the 
precision did not exceed +10'). 

From the goniometric data within the limits of 
precision, the indium fluoride crystals belonged to 
the tetragonal class. The following simple forms 
were found: prisms {100} and pinacoids {001}. 

The optical studies by E.E. Burova gave the 
following results: ny = 1.450 and Np = 1.425, 
straight extinction. A study of the figures in trans- 
mitted light showed that the crystals were not strict- 
ly uniaxial. When the stage of the microscope was 
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rotated, the cross was always observed to break in- 
to two hyperbolic branches; this indicates that the 
crystals were in fact biaxial, but with a very small 
angle between the optical axes, i.e., they were only 
pseudo-uniaxial. 

The lattice parameters were obtained from an 
oscillating x-ray photograph taken in Cu Kg radia- 
tion: a = 7.88 + 0.05 kX, c =4.13 + 0.05 kX, den- 
sity d= 2.83 g/cm’. 

Two formula units occur in the unit cell. The 
most probable space group C#;,—P4/n is defined by 
the following laws of extinction: hk/ , all; hk0, only 
with h+k=2n, 001, all. Other possible space 
groups are C{—P4, C}-P4/m, and S,—P4. 

Il. The F’-series method was used to discover 
the basic motif of the atom distribution in the indi- 
um fluoride structure. 

Two projections of the interatomic functions 
were constructed from x-ray goniometric data on 
the zero layer lines. The photographs were taken 
in a camera which photographed the reciprocal lat- 
tice with rotation of the crystals around the x and z 
axes (Mo Ky radiation). The photographs contained 
87 and 126 independent reflections. Intensity esti- 
mation was visual. Kinematic and polarization fac- 
tors were taken into account in constructing the pro- 
jections. 

The two indium atoms can take one of the four 
twofold positions of the group P4/n: (a) 000,%% 0; 
(b) 00%, 4%; ©) 0%, hoz. 

In the first two cases, the projection of the in- 
teratomic function on the plane (100) should have a 
maximum for In—In at the point v = 4, w = 0,while 
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Fig. 1. Projection of the interatomic function on the (100) plane. 


in the third case this should be at v = 44, w = 2z. 
The In—In maximum in the projection (Fig. 1) ful- 
filled the last requirement. 

Thus, the positions of the indium atoms are 
determined unambiguously; they occupy the (c) po- 
Sitions on the four-fold axis,z ~ 0.335. 

Since fluorine and water have almost the same 
dispersin power and cannot be distinguished by x- 
ray analysis, we considered them equal as a first 
approximation. The following possibilities are 


then available for the 12 light atoms (6F and 6H,0): : 


(1) they may occupy two of the three four-fold 
positions () 44,0, 4,%0, 4% 0, 440; 
e) os MyM tas nse ty hi (f) 002, 
00z, 4%z, %42Z, and two of the three 
two-fold positions (a), (b), and (c); 

(2) all three four-fold positions (d), (e), and 
(f); 

(3) the general eight-fold positions (g): xyz; 
Ryze AK ey, 2s Voy, 2. YXEs 
YXZ5 /, saw 3 yf, +X; Z; yf, Vis yf, —X,Z, and 
any two two-fold positions; 

(4) the eight-fold position (g) and one of the 
four-fold positions [excluding (f), of which 
more below]. 


However, if account is taken of the ionic radii 
of fluorine and oxygen and the lattice parameter along 
the z axis (rp = 1.33 A, ro2 = 1.36 A, ¢ = 4,13 A), 
the first two variants are immediately rejected, 
i.e., those in which two fluorine atoms (or two 
water molecules) lie one upon another along the z 
direction. The same applies to the four-fold position 
(f) in variant (4). 

An analysis of the projections of the interatom- 
ic functions on the (100) and (001) planes showed 
that they contain only maxima fulfilling the re- 
quirements of variant (3) with the two-fold positions 
(a) and (c) occupied, while the In—F maximum ne- 
cessary for variant (4) at the point v = 4, w = % on 
the (001) projection is not present (Figs. 1 and 2). 
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Fig. 2. Projection of the interatomic function on the (001) plane, 


The general pattern of atom distribution in the 
indium fluoride structure is as follows: eightatoms 
in the general position (g) with coordinates x = 
= 0.739, y= 0.113, z = 0.835, and two atoms occu- 
pying the position (c), z= 0.835,to form an octa- 
hedron around the indium atoms (zj, = 0.335). The 
atoms which occupy the vertices (c) of the octa- 
hedra belong simultaneously to indium atoms from 
neighboring cells, so that the octahedra are linked 
in endless chains parallel to the z axis. The cen- 
ters of individual octahedra (indium atoms) of neigh- 
boring chains are shifted relative to one another by 
about one-third of the cell edge (c). The chains are 
interconnected by atoms at the vertices and thecen- 
ter of the unit cell,i.e., the position (a). 

A more detailed analysis of the structure based 
only on x-ray data is impossible because the F and 
H,O cannot be distinguished. As shown below, we 
have been able to do this by crystallochemical ana- 
lysis. However, the experimental data are suffi- 
cient to answer one of the basic questions. The 
structure is a coordination type and does not con- 
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Fig. 3. Projection of the electron density on the (100) plane. 
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tain isolated InX, complexes, but the coordination 
number of indium is six. 

To confirm the correctness of the chosen struc- 
ture variant, a projection of the electron density on 
the (100) plane was constructed (Fig. 3). 

Figure 4 shows the agreement between the ex- 
perimental and theoretical structural amplitudes 
Form — For the values of which were initially re- 
lated to an absolute scale with the introduction of 
temperature corrections. Calculation of the coeffi- 
cient 


=Fel—| Fell 
Re € t 
el 


for the h0l reflections gave a value R= 0.15. The 
temperature correction constant was B= 2.0. 

The atomic coordinates cited in Table 1 were 
found from the projections of interatomic functions 
and refined from the electron-density projection. 

The probable error in the determination of the 
coordinates of the indium atoms was 0.01 kX, and 
for the light atoms, about 0.1 kx. 

Il. It has been noted already that the x-ray 
analysis data are insufficient for a complete solu- 
tion of the indium fluoride structure. The question 
of which positions in the structure [(a), (c), and 
(g)] are occupied by which light atoms still remains. 
By turning attention to the specific crystallochemi- 
cal peculiarities of water and fluorine, an unam- 
biguous answer can be obtained. Two initial posi- 
tions arise. 

1. The connections between the individual oc- 
tahedra in a chain are by atoms at position (c) 
which simultaneously belong to two indium atoms 
and have a coordination number of two. This posi- 
tion is characteristic of halogen atoms, but is com- 
pletely excluded for water molecules. 

Thus, the (c) positions should be occupied by 
fluorine atoms (F]). The In—Fy distance is 2.07 A 
less than the sum of the ionic radii of In and F 
(2E270 A), but this is quite natural, because the 
fluorine atom in this position has a coordination 
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number of two. Similar coordination of fluorine 
with simultaneous shortening of the M—F distance 
has been found in a series of fluoride structures, 
e.g., UF;, TaF3, MoF3, AIF; [2,3]. 

2. The chains of octahedra are interconnected 
by particles (F or H,O) in special positions of the 
inversion axis 4 [position (a)]. That this position is 
taken by water molecules (H,O7) follows from the 
considerations below: First, the distance to the 
nearest indium atom (4.18 A) is much too long for 
an In—F bond; second, the symmetry of this posi- 
tion indicates a tetrahedral environment which is 
frequently characteristic for the water molecule. 

These two facts determine that the remaining 
fluorine atoms and water molecules occupy the (g) 
positions. It is known that the tetrahedral distribu- 
tion of charge on the surface of water molecules re- 
quires definite packing for them in hydrate struc- 
tures. The latter are stable only if the charge field 
of the water molecules contacts either ions of op- 
posite charge or the oppositely charged regions of 
neighboring molecules. According to this, the 
molecules H,Oy should be surrounded by two fluor- 
ine atoms and two water molecules, oriented with 
the positively charged regions toward the center of 
the molecule. 

Consequently, each indium atom is surrounded 
in the plane perpendicular to the axis of the octa- 
hedral chain by 2 F and 2 HO— Fy and H,Oy (Figs.5 
and 6). The water molecules H,Oj and H,Oj; differ 
crystallochemically, since they occupy different 
regular systems of points. It was shown above that 
a strict tetrahedral environment was required for 
the H,Oj molecules. We shall now discuss the co- 
ordination about the H,Oyjz7 molecules and how the 
compensation of charges is brought about. 

The H,Oyz molecules form infinite chains ex- 
tended along the z axis. The chains are connected 
together by H,O] molecules in the corrugated net- 
work which consists of the octagonal rings and 
which is parallel to the (110) plane — the line AB in 
Fig. 5. This is shown diagrammatically in Fig. 5. 
One positive and one negative charge of the H,O]] 


TABLE 1 
Atomic coordinates 
Atom Position in kX in fractions of cell sides 
x y z x yu Z 
In (c) 3,94 1.38 0) 0.500 0.335 
(g) 0.89 2.06 1.38 0.413 0.264 0.335 
H,O,F (c) 0 3.94 3.45 0 0.500 0.835 
(a) () 0 0) 0 0 
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QD QraQn,@m 
Fig. 5. Crystal structure of 
InF, * 3H,O. View along the 
z axis (F] atoms are notshown). 
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Fig. 6. Crystal structure of InF 
* 3H,O. View along the x axis. 


molecule compensate the opposite charges of neigh- 
boring molecules in the chain; finally, there re- 
main‘a free negative charge caused by the bond to 
indium and a free positive charge on the H,Oy; mole- 
cule. 

Thus, both the H,Oy and H,Oy; molecules have 
tetrahedral coordination. 

The presence of two types of water molecules 
in the structure of indium trifluoride trihydrate 
(ratio of H,Oy : H,Oyy = 1:2) with the H,Oy molecules 
at a considerably greater distance from the indium 
atom than the H,Oyz molecules is probably the cause 
for the 5% loss in weight when the salt is dried at 
110° [1] (the loss of one molecule of water is 7.8% — 
a value of the right order); further heating of the 
salt to 188° caused no further change in weight, i.e. 
one molecule of water is evidently less strongly 
bound than the other two. 

It is now necessary to return to the following 
structural peculiarity: Both the indium atoms and 
the water molecules H,O] are placed on the four- 
fold axes, but each is surrounded by two pairs of 
different atoms,i.e., we have, in fact, simple diad 
axes in place of the axes 4 and 4 and, consequently, 
the space group P4/n must be replaced by P2/n of 
the monoclinic class. It is true that the monoclinic 
angle is very small — less than 90°10' (since the 
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precision of the goniometric measurements was 
only +10', it could not be observed). In all proba- 
bility, this lowered symmetry explains the pseudo- 
monoaxial nature of the crystal apparent in the op- 
tical study. It is possible that this effect is only 
weakly expressed because of the incomplete order 
in the positions of the fluorine and water in the 
actual structure, which leads to an increase in the 
symmetry. The direct extinction arises because 
the length of the needles coincides with the direc- 
tion of the two-fold axis. 

The equality of the n—F and n—H,O distances 
in the octahedra (2.08 A), which was obtained on 
the assumption of a four-fold axis, seems unlikely 
from a crystallochemical point of view. The same 
should be said of the somewhat decreased H,O—H,O 
distances (2.64 and 2.56 A). Since there is some 
analogy between our structure and that of ice, dis- 
tances of the order of 2.7-2.8 A would have been 
expected. 

Bearing in mind the development of a complete 
analysis of a structure of low symmetry, itseemed 
to us expedient to take somewhat different values 
for the coordinates of the water molecules which 
would eliminate the inconsistencies noted above. 
The x-ray analysis data permit this because, first, 
we averaged the intensities of independent reflec- 
tions, assuming a four-fold axis, when calculating 
the structure factors Fh; and, second, the (100) 
electron-density projection appeared imperfect be- 
cause of strong absorption in the crystal and the 
water maxima had a diffuse form which hindered 
the precise determination of their coordinates (this 
applies particularly to the coordinates of the water 
molecules). 

The final coordinates of the atoms in the struc- 
ture of indium fluoride trihydrate are given in 
Table 2. 


Fig. 7. A layer of H,O molecules in the InF, 
* 3H,O structure (section parallel to the line 
AB of Fig. 5), Small white circles represent 
(+) charges, hatched circles (-) charges. 
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TABLE 2 
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Atomic coordinates 


Atom Position in kX | in fractions of cell edges 
{ { 
a. | | IAA the 
| 
ia (c) 0 | 394 | 1.38] o | 0.500] 0.335 
H,0, | (4) 0 | —0.25 0 | O —0.061 
HO, | (g) 5.82 | 0.89 | 1.38 | 0.742 | 9.417 | 0.290 
y | ') 0 | 3.94 | 3.45 0 f 3: 
| 
| .38 2 ‘ 


The interatomic distances (in kX) are given be- 
low: 


In— Fy = 2.07, 
In— Fy = 2.08, 
In—H,Oyy = 2.13. 


H,O] — H,077 a 2 67, 
2.73 (distance 
H,Oy — H,Oy7= ‘ 
2 2 3.19 inchains) 


Finally, it is noted that there are some analo- 
gies between the structure of indium fluoride InF, 
* 3H,O and lead oxide PbO (space group Df, — 
P4/nmm) [4]. 


The structure analyzed here can be considered 
a complicated form of the PbO type: the Pb posi- 
tions are taken by the In atoms, while the oxygen 
positions are taken by those of the H,O; molecules. 
The increased complexity of the composition has 
led to a decrease in the symmetry in indium fluor- 
ide. 


Summary 


1. The crystals of InF; - 3H,O are pseudo- 
tetragonal, space group C3,—P112/n. The lattice 
parameters are: a =b= 7.88 + 0.05 kX, c= 4.13 
+ 0.05 kX, B = 90° + 10', N=2. 


2. The structural coordination does not con- 
tain isolated InX, groups, but the coordination num- 
ber of indium is six. 

3. A complete structural analysis was accom- 
plished by means of crystallochemical analysis, re- 
vealing the different F and H,O positions and the 
character of the packing of the water molecules. 

4. In agreement with its chemical behavior, 
the structure of indium fluoride trihydrate contains 
two types of crystallochemically different water 
molecules. Both types have tetrahedral coordina- 
tion. 

LITERATURE CITED 


1. E.N. Deichman and I.V. Tananaev, "Study of 
indium fluorides," Khim. Redkikh. Elementov, 
1, 95-101 (1954). 

2. J.A.A. Ketelaar, "Die Kristallstruktur der 

Aluminiumhalogenide, " Z. Krist. , 85,119-131 

(1933). i 

G.B. Bokii, Introduction to Crystal Chemistry 

[in Russian] (1954), pp. 386-401. 

4, G.B. Bokii, Introduction to Crystal Chemistry 
fin Russian] (1954), p. 359. 


ie) 
. 


SOVIET PHYSICS ~sCRYSTALLOGRAPHY 


VOL. 15. NOe 2 MAR.-APR., 1956 


METALLOGRAPHIC AND X-RAY-DIFFRACTION STUDY OF 
ALLOYS IN THE GERMANIUM-RHODIUM SYSTEM 


N. N. Zhuravlev and G. S. Zhdanov 


M.V. Lomonosov State University, Moscow 
Translated from Kristallografiya, Vol. 1, No. 2, 


pp. 205-208, March-April, 1956 


Original article submitted October 12, 1955 


A number of alloys in the Rh—Ge system were studied metallographically 
and by x-ray diffraction, and the existence of several compounds in this 
system was established. The unit cell of Rh;Ge, was determined. 


Preparation of Alloys 


Germanium and rhodium were fused together 
in various proportions. For preparing the alloys, 
chemically pure germanium and rhodium were used. 
The melting was effected in a high-frequency fur- 
nace and a resistance furnace in quartz ampoules 
in an argon atmosphere. 

The melting point was determined approximate- 
ly for a number of alloys. The shape of a section 
of the melting diagram for the rhodium —germani- 
um system was found on the basis of the melting 
points of the alloys and from metallographic and x- 
ray data (Fig.1). 


Metallographic and X-Ray Diffraction 
Study of the Alloys 


Metallographic study of germanium-rich alloys 
(above the composition RhgGe,) causes no great dif- 
ficulty. These alloys are easily etched by nitric 
acid. The study of rhodium-rich alloys beyond the 
composition RhgGe,, however, involves some diffi- 
culty, since such alloys cannot be etched by strong 
acids and alkalis, nor even by mixtures of HNO, 
+HCl and HNO;+HF. In order to observe the phase 
composition of such alloys, sections were prepared 
by the relief-polishing method. 

Alloys from germanium to 30% rhodium con- 
sist of the eutectic Rh3Ge,—Ge and germanium crys- 
tals. Figure 2 is! a micrograph of the alloy contain- 
ing 20 wt.% rhodium; etched germanium crystals 
can be seen on a field of eutectic. Alloys containing 
30 wt.% rhodium practically constitute the Rh3Ge, — 
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—Ge eutectic. Figure 3 is a micrograph of such an 
alloy. Alloys with germanium content in the range 
between the eutectic Rh3Ge,—Ge and Rh3Ge, are two- 
phase and consist of crystals of Rh3Ge, and the eu- 
tectic Rh3Ge,—Ge. Figure 4 is a micrograph of an 


toa 


100 


Temperature, °C 


10 320% (30.60) S50 S60 70 60 90 
Rh Wt.% Ge Ge 


Fig. 1. 


lRor Figs, 2, 3, 4, and 5, the magnification is 160; for Figs. 6 and 
he It as 804; 
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Fig. 5 Fig. 6 


Fig. 7 


Figs. 2 to 7. Microphotographs of alloys in the germanium—rhodium system. 


alloy of composition RhGe, (58.52 wt.% Ge), in 
which we see light crystals of RhgGe, in a field of 
eutectic. On an x-ray powder photograph taken in 
Cu Kg radiation for this alloy, we find lines belong- 
ing to RhgsGe, and germanium. 

Annealed alloys of composition Rh3Ge, (48.47 
wt.% Ge) are single-phase. On microsections from 
these alloys we find polyhedral crystals. The x-ray 
powder photograph obtained for Rh3Ge, crystals has 
a large number of lines. Alloys of Rh3Ge, cooled in 
the ordinary way are heterogeneous. We may sup- 
pose that the compound RhgGe, is formed by a peri- 
tectic reaction. In ordinarily cooled alloys of com- 
position Rh3Ge,, the reaction does not proceed to 
completion, and the microsections (Fig. 5) show un- 
converted crystals of RhGe surrounded by a shell of 
Rh,Ge, crystals. The space between such grains is 
filled with Rh3;Ge,—Ge eutectic. Alloys of composi- 
tion RhGe (41.36 wt.% Ge), both slowly and rapidly 
cooled, are single-phase and constitute the com- 
pound RhGe. The x-ray powder photographs ob- 
tained with slowly and rapidly cooled alloys are 


identical. This excludes polymorphism of the com- 
pound RhGe. Published indications [1] of the exist- 
ence of a compound with composition Rh3Ge, were 
not confirmed. The alloys of composition Rh3Ge, 
which were prepared, both annealed and rapidly 
cooled, were heterogeneous. These alloys con- 
sisted of crystals of the compound Rh;Ge3 and a 
small amount of Rh;Ge;—RhGe eutectic. Figure 6 
gives a micrograph of an alloy with composition 
Rh3Ge,, in which we see light polyhedral crystals 
of the compound Rh;Ge3 in a field of Rh,Ge3 —RhGe 
eutectic. Annealed alloys of composition Rh;Ge3 
are homogeneous and constitute the compound 
Rh;Ge3. 


Alloys of composition Rh,Ge (26.07 wt.% Ge) 
are metallographically single-phase and constitute 
the compound Rh,Ge. Rhodium-rich alloys consist 
of crystals of rhodium and Rh—Rh,Ge eutectic. 
Figure 7 is a micrograph of an alloy of composition 
15 wt.% Ge, in which we see light crystals of rhodi- 
um in a field of Rh—Rh,Ge eutectic. 
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TABLE 1 
Ay. micro- Scatter of 
On eaG é 
Compound Poend hardness, microhardness, 
g/c kg nm? kg/mm? 
Ge oro 520 420—600 
Rh,Ge, 8.5 500 420—575 
RhGe 9.8 580 470—660 
Rh, Ges 10.7 475 420—525 
RheGe 11.0 700 500—900 
Rh 12 4 390 320—420 


Density and Microhardness of the 
Compounds 


The density o of the compounds was deter- 
mined by hydrostatic weighing [2] in carbon tetra- 
chloride (¢ = 1.585 g/cm*). The table shows densi- 
ties o determined at room temperature. The den- 
sities of germanium and rhodium are taken from 
data in [3,4]. 

The microhardness of the compounds was 
measured on the PMT-3 apparatus with a constant 
load of 10 g on the indentor. Measurements were 
made on several samples from different melts and 
the results averaged. The conditions of measure- 
ment and surface preparation were as far as pos- 
sible identical. Table 1 gives the average micro- 
hardness and its scatter. 


Solubility of Rhodium in Germanium 
and Germanium in Rhodium 


For studying the solubilities of the components 
in each other we prepared alloys with 10 wt.% of 
the component being dissolved. After prolonged an- 
nealing (50 h), the metallographic alloys showed 
heterogeneous structures. X-ray powder photo- 
graphs were obtained in the 114.4-mm diameter 
RKU camera for these alloys in Cu Kq radiation. 
On indexing and computing the x-ray photographs it 
proved that in neither the first nor the second case 
had the lattice spacings of the rhodium and ger- 
manium changed. This indicates that the compon- 


ents have only infinitesimal solubility in each other. 


Determining the Unit Cell of Rh3Ge, 


Crystals of RhsGe, were produced by growing 
in the mother liquor with subsequent separation by 
dissolution of germanium. The separated crystals 
mainly constituted concretions and were of irregu- 
lar form. The most perfect were selected for x- 
ray study out of a large number of such crystalline 
concretions. The Laue photographs showed that 
Rh3Ge, crystals belonged to the tetragonal system. 


Oscillation x-ray photographs gave the lattice spac- 
ings asa = 5.7 + 0.2 A, c=10+ 0.3 A. 

An x-ray structural investigation of the com- 
pounds RhGe, Rh;Ge3, and Rh,Ge was made in [5]. 
Independent studies of rhodium-rich alloys agreed 
closely with our own results as regards the com- 
position of the compounds and their densities. The 
results of this work and Geller's data [5] disprove 
the existence of the compound Rh3Ge,. The tempera- 
ture of the transformation into the superconducting 
state (2.12°K) determined by Matthias [1] for alloys 
of composition Rh,Ge, probably relates to crystals 
of the compound Rh;Ge3. 


Conclusions 


1. The contours of a section of the melting 
diagram of the germanium—rhodium system have 
been marked out. 

2. Four compounds are found in the germani- 
um-—rhodium system: Rh3Ge,, RhGe, Rh;Ge3, and 
Rh,Ge. It is suggested that the compounds Rh3Ge, 
and Rh;Ge3 are formed by a peritectic reaction. 
The maxima on the melting diagram relate to com- 
pounds RhGe and Rh,Ge. On fusing the components, 
these compounds are formed with a considerable 
evolution of heat. 

3. No solubility of rhodium in germanium or 
germanium in rhodium is found. 

4, The density and microhardness of the com- 
pounds in the rhodium —germanium system have 
been determined. 

5. Rh3Ge, crystallized in the tetragonal sys- 
tem with lattice spacings a = 5.7 and c = 10 A, 

In conclusion, the authors wish to thank Pro- 
fessor I.E. Alekseevskii for discussing the results 
and I.I. Lifanov, R.N. Kuz'min, and N. P. Ivanova 
for help in conducting the experiments. 
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The possibilities of the electron-diffraction method (Kikuchi-line photographs) 
in studying the degree of perfection of the surface layer in single crystals are 
examined. The method is used for studying ground and etched single crystals 


of germanium. 


Single crystals of various semiconductors (ger- 
manium, silicon,etc.) have in recent years found 
extremely wide application in electronics and radio 
technology. The properties of these single crystals 
apparently depend not only on the insignificant quan- 
tity of impurities, but also on the perfection of their 
structure, in particular that of the surface layer. 
It is well known that in order to obtain small com- 
ponents for the manufacture of various pieces of 
apparatus, massive single crystals are sawed or 
filed. This clearly causes considerable damage to 
the single crystal structure in the surface layer in 
contact with the saw or file. In order to remove 
this layer, grinding and etching are usually em- 
ployed; the nature of the damage and defects caused 
in the sawing and grinding operations, however, is 
not entirely obvious in any specific case, nor is the 
depth of their penetration into the crystal. It is al- 
so not known whether etching removes the whole 
damaged layer or only part of it. 

In view of the above, a study of the structural 
perfection of such single crystals is,in our opinion, 
of considerable interest. 

In the present investigation, we studied the de- 
gree of perfection of germanium single crystals 
subjected to various forms of treatment by electron 
diffraction. Before proceeding to a description of 
the results obtained, we shall briefly discuss the 
possibilities of the electron-diffraction method in 
investigations of this type. 

In studying massive samples by electron dif- 
fraction, the reflection method is used. The dif- 
fraction pattern formed by the reflection of anelec- 


161 


tron beam from the surface of a single crystal must 
clearly depend a great deal on the degree of perfec- 
tion of the surface. It is in fact well known that in 
studying single crystals by electron diffraction, two 
types of diffraction record differing in principle,are 
obtained: photographs formed by a regular network 
of point reflections and photographs containing 
Kikuchi lines. As regards the electron-diffraction 
photographs of the first type, their formation, as 
we know [1], is associated with the mosaic struc- 
ture of the sample,i.e., with the presence of indi- 
vidual submicroscopic blocks turned through angles 
of about 1-3° in all directions with respect to each 
other. Electron-diffraction photographs of the se- 
cond type were first obtained by Kikuchi [2] in 1928, 
when studying the scattering of electrons by sheets 
of mica. In these photographs, in addition to asmall 
number of point reflections, he found pairs of par- 
allel lines (black and white), subsequently named 
"Kikuchi lines." The formation of the Kikuchi lines 
is associated with multiple (dynamic) scattering of 
electrons, and they can only be completely ex- 
plained on the basis of dynamic theory; Kikuchi, 
however, proposed the following elementary ex- 
planation for this phenomenon. 

Whenabeamof electrons passes through sheets 
of mica about 107° cm thick, intense scattering of 
the electrons takes place. A considerable propor- 
tion of the electrons lose a large part of their ener- 
gy as a result of multiple inelastic collisions with 
atoms of the lattice. These electrons form a gen- 
eral background. A certain number of the electrons, 
however, interact with lattice atoms without seri- 
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Fig. 1. Formation of Kikuchi lines (schematic). 


ous energy loss, although their direction of motion 
may change quite severely. 

As a result of such scattering, a widely diverg- 
ing beam of electrons capable of interference forms 
inside the crystal. Meeting a system of crystallo- 
graphic planes inclined at a suitable angle in the 
course of their journey, these beams of electrons 
experience regular reflection and, as a result of 
this, increased and decreased intensities develop 
on the background of inelastically scattered elec- 
trons. 

All this is illustrated by Fig. 1. Here, one of 
the scattered beams (I) falls on a system of planes 
at the Bragg angle, undergoes regular reflection, 
and falls on point P, instead of P,. Thus, at Py, 
there is a deficit of electrons as compared with the 
general background (white point), while at P, there 
is an excess (black point). Among the scattered 
beams capable of interference there may also (gen- 
erally speaking) be one (II) which is reflected from 
the other side of the same system of planes, falling 
on point P; instead of P,. The reflection of beam 
(II) would have compensated the excess of electrons 
at P, and the deficit at P,, were it not for the fact 
that the intensity of the scattered electrons dimi- 
nishes with increasing angle of deviation from the 
original beam (Ij > Iq), so that compensation does 
not take place. 

In order to obtain a line on the photographic 
plate, there must be a continuous set of scattered 
beams over a wide range of angle; the reflection of 
these beams leads to the formation of diffraction 
cones, the axes of which are perpendicular to the 
reflecting system of planes, and the aperture angles 
close to 180°. The intersection of these cones with 
the photographic plate leads to the formation of 
_ Kikuchi lines. 

Electron-diffraction photographs with pairs of 
lines were obtained by Kikuchi from mica sheets 
about 107° cm thick. Elementary consideration of 
the formation of Kikuchi lines shows that their pres- 
ence on the electron-diffraction picture is associ- 
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ated not only with the thickness of the sample exa- 
mined, but also with its degree of structural per- 
fection. 

In fact, for Kikuchi lines to be clearly visible 
on the photographic plate (distinguished from the 
background), the number of electrons falling on P, 
instead of P,; must be greater than a certain mini- 
mum. This minimum number depends not only on 
the volume of the crystal in which the scattered 
beams (I and II) are formed and their reflection from 
the system of planes, but also on the strictly paral- 
lel nature of these planes, one to another,i.e., on 
the perfection of the crystal structure. If the sample 
studied consists of individual blocks disoriented 
with respect to one another, or if in these blocks 
the planes are not strictly parallel to each other, 
then each block will reflect beams (I) and (II) into 
points differing slightly from P, and P,, so thata 
slight disorientation (of the order of minutes) in the 
orientation of the blocks will make the Kikuchi lines 
diffuse, and for a larger disorientation (1-3°) they 
will become so wide that the contrast between lines 
and background will entirely vanish,i.e., the lines 
will not be seen; in this case, the electron-diffrac- 
tion photograph will only show the point reflections 
formed by reflection of the primary beam S, from 


various systems of planes. 
If the sample being studied has no such blocks, 


and the reflecting systems of planes are strictly 
parallel to one another (perfect single crystal), a 
primary beam §) with an extremely small angular 
divergence has very little probability of falling on 
several systems of planes at the Bragg angle, and 
only one or two chance point reflections may form 
on the diffraction photograph. As regards the scat- 
tered beams (I and II), these, according to the fore- 
going discussion, lead to the formation of patterns 
with Kikuchi lines on condition that the dimensions 
of the sample in the direction of the beam Sy (or the 
thickness of the blocks protruding from the surface 
in the case of reflection photographs) are not 
smaller than a certain critical value. This value 
may be calculated in any specific case, as was done 
in [3]. The calculation leads to the value ~107~°cm. 
If the block size is smaller than critical, there will 
be no high-intensity scattered beams formed in 
them, or if these are formed they will have little 
probability of meeting a system of planes at the 
Bragg angle in the same block and giving reflec- 
tions at Py or P,. If, however, the block sizes are 
much greater than critical, the electron beamswill 
experience strong absorption and incoherent scat- 
tering, leading to an increase in the general back- 
ground and a weakening of the Kikuchi lines. 
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Thus, the presence or absence of Kikuchi lines 
in the electron-diffraction photograph, as well as 
their sharpness, depends on the degree of perfec- 
tion of the crystal: the greater the perfection of the 
crystal, the sharper and clearer will the Kikuchi 
lines be, and vice versa. 

In one paper [4], photographs with Kikuchi lines 
were used in studying the plastic deformation of 
aluminum single crystals. 

In the present investigation, the method de- 
scribed was applied to the study of the degree of 
perfection in ground and etched single crystals of 
germanium. The samples were ground with corund- 
um powder (120 min) and etched in hydrogen per- 
oxide. 

At the initial stage of the work, we made some 
special experiments with the aim of determining 
(albeit roughly) the thickness of the layer removed 
by the etchant. This value was found by weighing 
the sample before and after prolonged etching. As 
a result, it was established that, on etching in the 
peroxide at room temperature, 3 to 4 uw were re- 
moved per hour, while on etching in boiling per- 
oxide, 3 to 5 uw disappeared per minute. 

In studying germanium single crystals, it was 
important to obtain an electron~-diffraction photo- 
graph from the "natural" face, not subjected to any 
kind of treatment, so as to be able to compare it 
with photographs from the etched samples. Such an 
electron-diffraction photograph, obtained from the 
(111) cleavage face, appears in Fig. 2. 

The presence of sharp Kikuchi lines on this 
photograph, and the absence of a network of point 
reflections, indicates that, in the sample studied, 
there was no mosaic structure, but that in its de- 
gree of structural perfection it approximated to a 
perfect crystal. In this connection, it is interesting 
to note that in one paper [5] the (111) face of ger- 
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manium was studied microscopically and interfero- 
metrically, and that a small number of hexagonal 
depressions appeared, which the author described 
as mosaic structure. In our view, the depressions 
were not the result of mosaic structure, but had a 
quite different nature. 

Let us turn to an analysis of photographs from 
germanium samples treated in various ways. 

In the electron-diffraction photographs from 
filed samples, as a rule, there were point reflec- 
tions (together with strong background) sometimes 
accompanied by weak, short Kikuchi lines. These 
photographs led us to conclude that, when germani- 
um was filed or sawed, individual, quite large 
blocks (>107° cm) were formed in the layer in con- 
tact with the tool, these being turned relative to 
one another and to the original single crystal. On 
the photographs from ground samples (Fig. 3) the 
point reflections are accompanied by Debye half- 
rings, indicating the presence of polycrystalline 
dust on the ground surface. Thus, in grinding ger- 
manium, tiny crystallites (~10~° cm) are torn out 


Fig. 3. Electron-diffraction photograph from a ground ger- 
manium single crystal (L = 700 mm). 


Fig, 2. Electron-diffraction pkotograph from the cleavage 
face of a germanium single crystal (L = 250 mm). 


Fig. 4, Electron-diffraction photograph from a germanium 
single crystal etched for 15 sec (L = 700 mm), 
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of the crystal surface, together with the formation 
of mosaic structure. After brief etching (~15 sec) 
of the ground sample, the electron-diffraction 
photograph contains only point reflections, forming 
a two-dimensional network, or else set along 
straight lines (Fig.4). 

The presence of reflections up to the sixth 
order from the (111) face on the photograph of Fig.4 
corresponds, according to the Bragg —Wulff form- 
ula, to a block disorientation of about 1-3° (A =0.05 
SE OE SN. 

In the course of later work, we made a more 
detailed study of the effect of etching time on the 
degree of perfection of the surface layer of single 
crystals. This study was made on germanium 
samples, which after grinding were etched in boil- 
ing peroxide several times for 15 sec, each etch 
being effected in fresh etchant. As a result, it was 
established that a pattern containing Kikuchi lines 
appeared after a half-minute etch; at this stage, 
however, the Kikuchi lines were greatly diffused 
and accompanied by point reflections, indicating 
mosaic structure in the surface layer. Increasing 
the etch time to 1 min led to a slight increase in 
the sharpness of the Kikuchi lines, without any 


marked change in the extent of the mosaic structure. 


An analogous picture was found for samples etched 
for 2 min. The electron-diffraction pictures of 
samples etched for 4 min (Fig. 5) showed point re- 
flections up to the fourth order, indicating a block 
disorientation of ~1°30'. Apart from the reflections 
mentioned, this photograph also showed fairly 
sharp Kikuchi lines; comparison of this with the 
photograph corresponding to the cleavage face, how- 
ever, Showed that even the 4-min etch had not en- 
tirely removed the layer damaged in grinding. 

The electron-diffraction photograph shown in 
Fig. 6 was obtained from a sample etched for 10 
min (40-yu layer removed). Just as in the case of 
the cleavage face, this gave sharp Kikuchi lines 
and bands, indicating a fairly perfect structure for 
the surface layer. 

In some cases it was necessary to increase the 
etch time to 20 min in order to obtain such pic- 
tures; this corresponded to the removal of a layer 


80 pw thick. 
Thus, for complete removal of the layer with 


mosaic structure formed on grinding a germanium 
single crystal, etching in boiling hydrogen peroxide 
for 10-20 min is required. It also follows from the 
above data that the structural damage (mosaic) pro- 
duced in the surface layer of the single crystal on 
grinding penetrates to a depth of 40 to 80 wu. We 
note that on the diffraction photographs of germani- 


Fig. 5. Electron-diffraction photograph from a germanium 
single crystal etched for 4 min (L = 250 mm), 


Fig. 6. Electron-diffraction photograph from a germanium 
single crystal etched for 10 min (L = 250 mm), 


Fig. 7. Electron-diffraction photograph from a germanium 
single crystal etched for 4 min (L = 700 mm), 


um samples etched in boiling hydrogen peroxide 
with added alkali solution (5 drops of a 5% NaOH 
solution to 10 cm?’ H,O,) for 10-20 min, haloes from 
some kind of amorphous film formed on the surface 
of the single crystal during etching were system- 
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atically observed in addition to the Kikuchi lines 
and bands. 

In conclusion, let us examine briefly the 
method of obtaining the diffraction photographs. In 
the reflection of an electron beam from the surface 
of a single crystal with perfect structure, a pattern 
containing Kikuchi lines is formed, the diffracted 
beams extending over a large angular range. Fora 
more complete recording of this picture, it is de- 
sirable to operate with a small distance between 
sample and photographic plate. This, however, in- 
volves the inconvenience associated with the small 
scale of the diffraction picture. 

In view of this, we took photographs with two 
different positions of the crystal-holder, with L 
= 250 and 700 mm, respectively. Comparison of the 
diffraction photographs shown in Figs. 2 and 7 dem- 
onstrates the advantages of this method. In the dif- 
fraction picture of Fig. 7 we notice the curvature of 
some Kikuchi lines, associated with the refraction 
of the electrons [1], and the extremely high inten- 


165 


sity of the "mirror" reflection situated on the 
Kikuchi line. 
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Unit cells (all monoclinic) are reported for 56 specimens of various hydro- 
micas; the parameters inthe basalplane are a = 5.16-5.29 kX, b = 8.90-9.20 


kX. 


Specimens with a repeat consisting of a single three-level silicate layer 


have c = 10.1-10.3 kX and 6B = 99.5-101.5°; the two-layer hydromicas have c 


= 20.0-20.6 kX and B = 93.5-96°. 


The diffraction patterns and layering allow 


the specimens to be classed in five groups. 


The hydromicas fall between the micas and the 
montmorillonites as regards structure and physico- 
chemical properties; the amounts of K and H,O are 
largely responsible for this. The K cations lie be- 
tween the silicate layers, and their number con- 
trols the adhesion between the layers, which thus 
is largest in the micas and least in the montmoril- 
lonites. 

The H,O in excess of that required by the struc- 
tural pattern is considered [1,2] to consist of H,O 
ions, whose position and role resemble those of K. 

Our ideas on the structures of the hydromicas 
and montmorillonites are rather general, because 
these minerals are highly variable, as are clay 
minerals generally. This means that any structur- 
al characterization must be based on a large group 
of specimens, which should include (apart from the 
various hydromicas) specimens of micas and seri- 
cites as being representatives of the layer silicates 
having the most perfect structures. We have exam- 
ined specimens previously studied in other ways 
and also mineralogically identified specimens, 
mainly from M.F. Vikulova's collections. 

Clay minerals lose the interlayer water during 
examination by electron diffraction, so there is no 
essential difference in structure between hydrom- 
icas and montmorillonites. The diffraction patterns 
are Similar, and the detail is governed by the struc- 
tural perfection (see Table 1). 

The specimens have been arranged in order of 
the qualitatively evaluated degree of structural per- 


fection. The first 37 have a more or less strict c 
repeat distance, as the reflections with k ~ 3k' are 
fairly readily distinguished. The other specimens 
have only approximate spatial repeats, because only 
the reflections with k = 3k' are clearly distinguish- 
able [3]. 

All the patterns have largely the same set of 
reflections: 131, 131, 201, 201 (second ellipse) and 
331, 331, 061 (fifth ellipse), which represent 
periodicity within a single type of layer with three 
levels. On the other hand, the 111, Lee , 021 reflec- 
tions (first ellipse) show several types of array as 
regards number and intensity. The reflections 
with k = 3k' indicate that a monoclinic cell should 
be adopted for all specimens; then the k ~ 3k' re- 
flections indicate that the layers may have a repeat 
of one layer or two (see the c repeat distances).The 
c and 8 then make the monoclinic cell with a one- 
layer repeat equivalent to an orthorhombic cell with 
a repeat of three layers (because cos B © a/3);and 
the monoclinic cell with a two-layer repeat is 
equivalent to an orthorhombic one with a repeat of 
six layers. These unit cells can be correlated with 
all the structural modifications of the micas esta- 
blished by x-ray analysis [4], each in accordance 
with the spatial symmetry of the atomic array, 
apart from the modification with a 24-layer repeat. 
This requires analysis of the intensities,of course 
_ A purely qualitative examination of the 021,11 
11/ reflections shows that there are several differ 
ent groups, each group having similar intensities 
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TABLE 1. Unit Cells of Specimens 
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Experiment : 

xo. Specimen No, Preliminary identification B° 

1 289 95.5 

2 3050 5.16 8.90 20.3 96.0 

3 Zh-14/27 r Sericites 5.19 8.96 20.1 96.0 

4 Zh-14 5.16 8.92 20.3 96.0 

5 Zh-12/6 5.18 8.98 20.3 95.5 

6 736 Muscovite 5.16 8.90 20.0 95.9 

7 Weathered muscovite peli 8.90 20.0 96.5 

8 Phlogopite 5.29 SHINS 10.3 100.0 

9 Biotite 5.28 9.20 LO, 100.0 
10 Vermiculite 5.30 9,20 10.4 99.5 
11 Zh-12/3 5.18 8.96 20.3 94.5 
eZ 756 Sericites 5.18 8.98 20,3 94.0 
tS 3060 els 8.92 20.3 96.0 
14 577 Glauconite Bole 9.04 10.1 HOU S 
15 342 Hydromuscovite 5.18 8.95 20e2 95.0 
16 Zh-4/2 te 5.18 9.10 IO, 3 95) 
17 Zh-1/4 eee 5.17 8.95 20.1 95.0 
18 L-11 Hydromica 5.18 8.95 10.2 100.0 
19 Zh-12/4 Sericite 5.20 9.00 20.2 95.5 
20 769 Kaolin clay 5.20 9.00 20.6 94.0 
21 420 eydromicas Sah 8.94 10.3 100.5 
22 118 5.19 8.97 10.2 101.0 
DS 750 Altered kaolinite 5.16 8.92 20.6 93.5 
24 122 Hydromica 5.19 8.96 10.3 INOULE') 
25 T-21 Hydromica or chlorite 5.20 9.00 10.3 100.0 
26 $-140 Hydromica 5.18 8.98 10.3 101.0 
PA | T-19 Chlorite 5220 9.00 OR 100.5 
28 N-76 5.20 9.00 10.3 100.5 
29 S-150 : o.20 8.99 10.2 100.0 
30 492 ey couicas 5.18 8.96 10.1 100.5 
Bi L-9 5.20 9.00 10.0 100.0 
ay 241 Quaternary hydromica 5.19 8.95 10.2 100.0 
33 OND 5.18 8.95 10.2 100.0 
34 iif 5.18 8.96 20.0 96.0 
35 124 ela 8.94 20.0 95.0 
36 453 Hydromica clays 5.18 8.96 20.0 95.0 
37 415 5.18 8.98 10.2 100.0 
38 S-151 5.18 8.98 10.1 100.0 
39 134 5.18 8.95 MOR2 100.5 
40 273 Devonian hydromica Oplll 8,94 OZ 100.0 
41 478 5.20 9.00 10.0 100.0 
42 81 TyarOtmentele? 5.16 8.91 10.2 100.5 
43 459 5.18 8.95 10.1 101.0 
44 3875 5.18 8.95 10.2 100.5 
45 488 Devoiian hydromica 5.20 9.00 HOSZ 100.0 
46 K-110 5.17 8.94 10,2 100.5 
47 B-564 5.18 8.94 O33 100.5 
48 388 5.17 8.94 iL. 100.0 
49 B-541 ‘ 5.18 8.95 10.3 100.0 
50 446 ee 5.20 9.02 10.1 100.0 
51 458 5.17 8.95 10.2 100.0 
52 569 5.20 9.02 10.0 100.0 
(513) 144 5.18 8.94 OZ 100.0 
54 940 Quaternary hydromica 5.16 8.94 LOS: 100.5 
55 443 f 5.18 8.98 10,2 100.0 
56 387 Beare cay 5.17 9,02 10.2 100.0 
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and positions for the 021,111, 11/1 reflections, 
while groups differ either in intensity or in position. 

Group I includes 289, Zh-14/27r, Zh-14, 
Zh-12/6, 736, the weathered muscovite, 342, and 
769. The arc reflections for 289, 342, and 769 are 
much shorter than those for the others of this 
group, which indicates that the flakes of these three 
are isometric, the others being elongated. 

Group II includes 3050, 3060, and (perhaps) 
Zh-12/4 and 750, whose reflections lie roughly as 
for group I, but which have a distinctive intensity 
distribution. 

Group III consists of the phlogopite, biotite, 
vermiculite, and Zh-4/2; the patterns demonstrate 
a strict one-layer repeat. The c and # indicate 
that the interlayer water is lost from vermiculite 
as it is from montmorillonite and halloysite [5,6]. 

Group IV (Zh-12/3, 756, and Zh-1/4) has a 
distinctive intensity distribution, the reflections 
being especially uniformly distributed in height 
above the zero line. The unit cells given in the 
table could be replaced by orthorhombic ones of c 
repeat extended to six layers. 

Group V is the most numerous and consists of 
specimens of relatively low structural perfection, 
which give patterns of glauconite type; 577, L-11, 
420, 118, 122, T-21, S-140, T-19, N-76, S-150, 
402. 1-9, 212. ete. 

It is difficult to establish the distinctive 
features of the specimens Nos. 38-56 because the 
patterns have either very indistinct k ~ 3k' reflec- 
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tions or ones that completely overlap others. 

The patterns indicate as hydromicas some spe- 
cimens whose preliminary identifications were dif- 
ferent (769, 750, T-19). 

The a and b parameters of 3050, Zh-14, 736, 
etc., are somewhat unusual for micas. These are 
higher or lower than the values given in the litera- 
ture for phlogopite, biotite, and vermiculite; this is 
due to the isomorphous replacement within the struc- 
tures of the particular specimens. 

I am indebted to M.V. Vikulova for making 
available the various materials and for valuable in- 
formation on the mineralogical characteristics of 
the specimens. 
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Optical characteristics are reported for the UFP series of ultraviolet 
polarizers, which operate between the limits of 240 and 750 my. 


Polaroid was invented some twenty years ago 
and has led to rapid extension of the use of pola- 
rized light in science and technology. Existing 
polaroid materials have excellent characteristics 
in the visible region and equal nicol prisms in this 
respect, but they absorb strongly in the ultraviolet 
and are largely useless in that region, although 
there is a great need for thin film polarizers of 
high aperture for the ultraviolet. 

Research has been done at this institute on 
ultraviolet polaroids based on polyvinyl alcohol 
films, which are of adequate transmission in the 
range 240-400 mu. These films produce polarized 
ultraviolet after treatment with suitable reagents. 
The products are neutral in color and transmit a 
little blue light in the crossed position. The region 
of best ultraviolet polarization is dependent on the 
treatment; some specimens also polarize visible 
light and so may be used with or instead of ordinary 
polaroid. 


Methods 


All the spectral measurements were made with 
an SF-4 photoelectric spectrophotometer. All such 
instruments polarize the radiation to some extent; 
this is unimportant with ordinary materials, but 
serious errors can arise with dichroic materials. 
The simplest method is to use polarized light inall 
measurements, by inserting a polarizing prism in 
the system. A polarizing device! was used, which 
consisted of a polarizing prism and a means of ro- 
tation. The prism was installed directly behind the 
exit slit in a support that could be rotated via a gear 
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Fig. 1. Dichroism curves for UFP specimens 
made in various ways. 


train through 180° about an axis parallel to the beam; 
the angular setting could be read to 15'. The speci- 
men was examined with its plane of polarization 
first parallel to that of the prism and then perpen- 
dicular to it. The required orientation was obtained 
by setting the prism in the position giving minimal 
transmission; the accuracy was increased by insert- 
ing a neutral filter of 1-2% transmission in the re- 
ference beam. All specimens were used uncemented 
and were 0.07-0.08 mm thick. Reflection losses 
were neglected. The spectral transmission was 
measured to 0.5%. 


Optical Characteristics 


The UFP’ series of ultraviolet polaroids was 
examined; Fig. 1 shows dichroism curves for four 


I This device was made at the P/Ya 412 plant at our request. 
Ultraviolet film polarizers. 


169 


170 Gry Le 


I 


300 400 500 600 700 800 X,mu 
Fig.2, Dichroism curves for UFP No, 4, 


specimens made in different ways. These differ in 
properties; each gives the best polarization over 
some particular range. For instance, No. 1 may 
be reckoned the best for 240-265 mu, No. 2 for 
265-290, No. 3 for 290-340, and No. 4 for 340-400 
my. The last is also suitable for the visible region 
(Fig. 2). 

In 1949, Barer [1] reported the ultraviolet 
transmission of British polarizing films; Fig. 3 
shows his results. This material gave complete 
polarization only over the range 290-300 my ,and 
the transmission was very low throughout the ul- 
traviolet. 

It is clear that the new polarizers are much 
more efficient. The set should be used together, 
each material being employed only in its best 
range. The set as a whole gives reasonably high 
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Fig. 3, Dichroism curves for British polaroid, 


transmission and a high degree of polarization 
over the range 240-750 my. 

These new materials should find application in 
optical instruments for many purposes requiring 
polarized ultraviolet light in science and technology. 


Conclusions 


1. New ultraviolet film polarizers have been 


developed. 
2. The optical characteristics of these are re- 
ported. 


I am indebted to E.V. Parvova for preparing 
the materials and for considerable assistance in 
the work. 


LUE RA TU Rae .Cle tp 
1. R. Barer, "Some experiments with polarizing 


films in the ultraviolet," J. Sci. Instr. ,26,10, 
325 (1949). sy 
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for the visible region. 
from conjugated double bonds. 


Introduction 


Russian polaroid is made by iodinating films of 
polyvinyl alcohol and has good optical character- 
istics [1], the degree of polarization and transmis- 
sion being high throughout the visible region. The 
color is almost neutral or very slightly greenish. 

Figure 1 shows the dichroism. The molecular 
compound of iodine with polyvinyl alcohol is therm- 
ally unstable; the material is damaged if heated 
above 80-90°C, which makes it unusable in vehicles, 
stereoscopic projectors, and in any optical equip- 
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ment developing much heat. 

There is clearly a need for better heat resist- 
ance in combination with good optical character- 
istics. Two reviews by Land, the inventor of pola- 
roid [2,3] refer to heat-resisting polaroid (type K) 
for use in vehicles, which is made by thermally de- 
hydrating films of polyvinyl alcohol. Polaroid K is 
of excellent optical quality; there is some fall in 
the degree of polarization in the red region, so the 
material transmits a red color in the crossed posi- 
tion. Figure 2 shows the dichroism. Land's re- 
views give only the very general statement that po- 
laroid K is made by heating stretched films of poly- 
vinyl alcohol to 160°C in the presence of a suitable 
catalyst; it is assumed that the dichroism is due to 
conjugated double bonds. 

Ushakov et al. [4], in 1949, made a fairly de- 
tailed chemical study of the dehydration of poly- 
vinyl alcohol. 

Films of this type have been made in this insti- 
tute and have been found to be of good thermal sta- 
bility and to act as polarizers in the visible region. 


Fig. 1. Dichroism ratio of iodinated polyvinyl] alcohol films. 


The structural changes during heating are very 
complicated and are dependent on the pH, degree of 
stretching, catalyst, and other conditions; they may 
be followed from the infrared absorption as meas- 
ured with polarized radiation. It has proved possible 
to make heat-resisting polarizers with good optical 
characteristics; the origin of the dichroism has al- 
so been established. 


Methods 


Alkali-saponification polyvinyl alcohol (made 
by the Polymeric Products Research Institute, 
Ministry of the Chemical Industry) was used; films 
0.1-0.15 mm thick were made by evaporating 8-10% 
aqueous solutions on glass. These films are very 
highly elastic and are readily oriented by stretch- 
ing, especially when hot. They were stretched to 
5-6 times the initial length in a frame with a worm 
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Fig. 2. Dichroism of American Polaroid K, corrected for 
reflection losses, 


drive; they were heated in a drying oven to various 
temperatures and over an electric hotplate, the 
temperature in the latter case being 100-110°. 

Deformation during heat treatment was avoided 
by clamping the films in holders in the stretched 
state. The average thickness was 0.07-0.08 mm. 

The spectrophotometry was done with an SF-4 
(quartz) spectrophotometer fitted with a polarizer. 
Infrared absorption spectra were recorded with an 
IKS-1 spectrometer (rocksalt prism), the polarizer 
being of KRS-5. The spectrum was recorded photo- 
graphically; spot measurements were also made. 
The films were measured unmounted; no corrections 
were applied for reflection. The errors of measure- 
ment were 0.5% (SF-4) and 1-2% (IKS-1). 


Results 


These films were found to be extremely stable 
on heating; no water was lost unless a special cata- 
lyst was added. The optical quality is not good un- 
less the water loss is fairly low and is produced 
rather rapidly, because numerous oxidation prod- 
ucts are readily formed, as can be seen from the 
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Fig. 3, Dichroism of specimen 115. 
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Fig. 4. Dichroism of specimen 64 before and after second 
stretching, 


C=O absorption (the heating was done in air, not 
inert gas). 

Various catalysts were tried: HCl, HBr, FeCls, 
ZnCl,, AlCl3, etc. The best results were obtained 
with HCl, HBr, and AICl;. 

HCl was used in most cases, this either being 
added to the solution before the films were poured, 
or being used as vapor with the stretched film. 
Temperatures of 110-160° were used, with treat- 
ment times of 30-60 sec. 

Figure 3 shows the dichroism of specimen 115, 
which was made at 160° (30 sec) with HCl as cata- 
lyst. The dichroism is high, but the overall trans- 
mission is poor. 

Better characteristics are obtained if the stress 
in the dehydrated material is allowed to relax, the 
film then being stretched again. The relaxation may 
be produced in hot water or in 1-5% aqueous boric 
acid. The contraction factor is 2-3. Figure 4shows 
the improvement (the total transmission shows 
about 10% increase), which arises because the 
stretching increases the degree of orientation. 
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Fig. 5. Dichroism of specimen 80 before (1) and after (2) 
heating to 130° for 250 h. 
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Fig. 6. Dichroism in the Tange 5.5-95.6 
of three pieces of specimen 52 set parallel 
(each pieces 0,03 mm thick), 


This method gives materials of very low de- 
gree of dehydration. More prolonged treatment 
produces a dark material of low transmission, so 
it is necessary to remove the catalyst and stabilize 
the material. The following stabilizers were tested: 
NaCl, Na,CO3, Na,B,O7, and H3BO3 mixed with 
Na,B,O,. Aqueous solutions of these remove the 
catalyst and enter the film, which relaxes substan- 
tially (mo special relaxation treatment is needed); 
the film is then stretched again, the transmission 
rising. 

The products were tested for thermal stability 
in a thermostat at 125-130°. Figure 5 shows the re- 
sult for specimen 80 after 250 h; the general trans- 
mission was increased, but the degree of polariza- 
tion in the red region dropped. 

The optical parameters change only during the 
first 20 h and thereafter remain virtually constant. 

Infrared spectra reveal the source of the di- 
chroism, as they are very sensitive to the struc- 
ture of macromolecular films; polarized radiation 
enables one to determine the orientations of groups 


Jo 

we II 
40 

ft 
30 
20 
10 

: {__>»— 
400 500 600 700 A,Mp 


Fig. 7. Dichroism of specimen 52. 
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Fig. 8. Dichroism at 5.5-6.5 y for a 0.03- 
mm stretched film of polyvinyl alcohol. 


and bonds. Dehydration has a complex mechanism 
very much dependent on the conditions used. ! 

Treatment in air at 110-160° always gives rise 
to double bonds and to the production of fresh carb- 
onyl groups. 

Figure 6 shows the dichroism in the range 5.5- 
7 w for specimen 52 (dehydrated for 5 min at 160°), 
which has pronounced dichroism in the visible re- 
gion (Fig.7). The infrared absorption is compara- 
tively weak, so three pieces of a specimen were 
used in the parallel position. Figure 8 gives the 
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The detailed mechanism is dealt with in a paper presently in 


preparation, 
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curves for untreated film. The band at 5.8-5.9 uw in 
both spectra indicates that C—O groups are pres- 
ent initially and after treatment, and that these 
have no obvious preferred orientation. Heat treat- 
ment gives rise to a new broad band at 6.1-6.3 yn, 
which corresponds to double bonds parallel to the 
stretching direction. Various types of double bond 
are present. 

More prolonged treatment gives a clearer pic- 
ture; specimen 70 (1 h 40 min at 150°) gives sharp 
dichroic bands at 5.9 and 6.3 yw (Fig.9), the latter 
corresponding to conjugated double bonds of the 
type -CH= CH—CH= CH—CH=CH_-, which lie 
parallel to the stretching direction. The C=O 
groups are perpendicular to this direction, which 
shows that the heat treatment produces crosslink- 
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ing. The second stretching clearly has tie functions 
of reducing the number of these links and of in- 
creasing the degree of orientation parallel to the 
stretching direction. These conjugated double bonds 
are thus the chromophore responsible for the di- 
chroism in the visible region. 


Conclusions 


1, Heat-resisting polarizers have been made 
from dehydrated polyvinyl alcohol films. 

2. A method has been devised for improving 
the optical parameters and for stabilizing the ma- 
terial up to relatively high temperatures. 

3. Infrared spectra show that conjugated 
double bonds are responsible for the dichroism in 
the visible region. 
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A description is given of the principal features of a method of growing large 
crystals of this tartrate, together with a novel method of orienting and mount- 
ing seeds; the products are virtually free from inhomogeneity caused by re- 


generation. 


The compound ethylene diamine tartrate (EDT) 
has crystals that can be cut as crystal filters of low 
temperature coefficient. 


The piezoelectric properties were first exa- 
mined by Mason in 1947 [2]; EDT crystals are now 
widely used abroad, and a plant solely devoted to 
this has been built in the USA [3]. No very exten- 
sive use has been made of EDT here, as the materi- 
al has been rejected in favor of potassium tartrate, 
which has some advantages in crystal filters. How- 
ever, EDT has some other uses, and the production 
of its crystals has not previously been described in 
the Russian literature, so some details are given 
here of techniques of value in large scale produc- 
tion. In addition, some new information is presented 
on the growth of the crystals, which is of interest 
in relation to the growth of other substances. 

The results presented here were obtained by 
the author jointly with A.A. Shternberg between 
1947 and 1952. 

The compound is not available for purchase in 
finished form, so users have to prepare it. I.F. 
Bogdanov, in 1947, devised methods of making EDT 
in the form of aqueous solutions and small crystals 
with ethylene diamine hydrochloride as starting ma- 
terial. 

Ethylene diamine tartrate CgH,4N,O, is the 
neutral salt resulting from the interaction of tart- 
aric acid with ethylene diamine. The crystals ob- 
tained from aqueous solutions usually contain no 
water of crystallization, but a monohydrate is also 
known, and this is deposited from aqueous solutions 
below 40.6°C. 
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Only anhydrous EDT crystals are used; the 
simplest method of preparing the salt is to neutral- 
ize tartaric acid with an aqueous solution of ethyl- 
ene diamine, the latter being available as aqueous 
solutions containing up to 70% of the diamine of tech- 
nical purity. The neutralization produces much 
heat, so the reaction must be carried out by slow ad- 
dition of the acid to the solution of the base. The 
temperature must not be allowed to rise above 50- 
55°C, because the base starts to decompose above 
this. Some decomposition occurs even below 50°C, 
but this is relatively slight; the solution becomes 
progressively more yellow as time passes. The de- 
composition products are removed by filtering the 
solution through a layer of activated charcoal, the 
filtrate being almost colorless. Exposure to heat 
and light rapidly causes it to turn yellow again. 

EDT solutions are highly viscous; the literature 
gives the viscosity at 50°C for a saturated solution 
as about 240 millipoises. The salt is very soluble in 
water, the solubility increasing substantially with 
temperature; a liter of water dissolves about 1.3 kg 
at room temperature and about 1.6 kg at 50°C. EDT 
crystals can thus be grown by cooling the solution, 
especially since the solubility is a linear function 
of temperature. The decomposition products tend 
to reduce the solubility. 

Practical purposes are satisfied by the above 
approximate results, which relate to solutions of 
pH 7.0-7.2, which corresponds to a slight excess 
of the base. 


*Readat the First All-Union Conference on Piezoelectricity, Insti- 
tute of Crystallography, November 26, 1952 [1]. 
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EDT will form highly supersaturated solutions, 
as for many salts that form highly viscous solu- 
tions; no crystals are formed even when the solu- 
tion is supercooled by 12-15°C. Here we disagree 
with Mason [3], who states that EDT solutions show 
a considerable tendency to crystallize spontaneous - 
ly (more so than Rochelle salt or ammonium phos- 
phate); our observations indicate precisely the 
converse. 

EDT crystals fall in the dihedral axial class of 
the monoclinic system (Fig.1); 6 is 105°30'. The 
density is 1.54, and there is a prominent cleavage 
on {001} 3 

We used small crystals in solutions at rest to 
establish the main features of the growth. There 
is a pronounced tendency to preferential growth in 
the negative sense of the b axis (a diad symmetry 
axis). The {110} faces around the other end of this 
axis scarcely grow at all, as do the pinacoid faces. 
EDT resembles ammonium phosphate in usually 
giving needle crystals, which are of no practical 
value. Crystals of large cross section require the 
use of broad seeds, which must also be cut from 
such large crystals. The best seeds are the heads 
of broad crystals, which consist of the rapidly 
growing dihedral faces. These broad crystals were 
produced by a method previously employed for am- 
monium phosphate, in which highly supersaturated 
solutions are employed to produce parallel growths 
of small crystals on the pinacoid faces, which ac- 
centuate the growth of the latter. Several growth 
cycles based on these heads give reasonably broad 
crystals with heads suitable for use as seeds. No 
more effective method of making broad crystals 
(e.g., adjustment of solution composition, as for 
ammonium phosphate) is presently known. We ini- 
tially used heads cut perpendicular to the b axis and 
attached with cement to a flat horizontal disk. The 
crystals grew along the b axis (Fig.2); the effects of 
supersaturation on the growth were examined on 
these. Highly supersaturated solutions gave only the 
{oo1} and {100} pinacoid and {110} and {110} di- 
hedron faces. The faces of the {011} dihedron grew 
at the positive end. Low supersaturations caused 
some growth on the {101} pinacoid, {111} dihedron, 
and {010} monohedron faces. Some faces of types 


P. G. POZDNYAKOYV 


Fig. 2. 


(hkl) and (h0l) were also observed, especially in 
impure solutions. 

Inhomogeneities are produced within a certain 
range of supersaturations. EDT crystals were 
formerly grown at low supersaturations, which fall 
short of that range, but the growth is then only 0.8 
mm/ day, so it requires about three months to grow 
a reasonably large crystal [4,5]. The lower end of 
the hazardous range corresponds to a supercooling 
of 3-4°C. We have found [6] that the upper end of 
that range corresponds to supercooling by 6-8°C; 
supercooling by 12-15°C is possible, and this gives 
growth rates of 5-8 mm/ day, which means that a 
crystal weighing over 1 kg can be grown in about 20 
days. Lower supercooling gives inhomogeneity be- 
cause additional faces arise, these having growth 
rates sharply dependent on the supersaturation; they 
give rise to inclusions of solution. In particular, a 
fast-growing {010} face usually produces cloudy 
layers parallel to itself; this face is favored by me- 
dium supersaturation, and its growth rate increases 
very rapidly with the latter, which depletes the so- 
lution, the residual mother liquor becoming trapped. 
These additional faces are lost at high supersatura- 
tions, so no cloudy inclusions are produced. 

EDT crystals regenerate readily, so any ac- 
cidental damage is healed and does not give rise to 
large or persistent defects. 

We then made tests on the effects of impurities 
and solution composition generally on the growth. 

Neutral solutions gave falling growth rates and 
tapering crystals, very much as occurs with am- 
monium phosphate if the solution is too acid. 

Alkaline solutions (ones with excess of the 
base) give higher rates and eliminate taper. We 
have not tested the effects of other bases, in par- 
ticular ammonia, which has been recommended by 
American workers; Walker [5] states that ammonia 
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sufficient to give pH 7.5 eliminates tapering and 
doubles the growth rate (to 1.5 mm/ day). We have 
produced rates more than three times this, so we 
did not examine the use of ammonia. 

Experience over more than four years has 
shown that EDT crystals are very insensitive to a 
great many impurities, whereas the growth rates 
of ammonium phosphate and potassium tartrate are 
very sensitive to impurities. 

These features make it especially simple to 
grow large crystals; the initial supersaturation 
need not be set very precisely, nor need the rate of 
temperature reduction be kept strictly constant. Al- 
so, the solution does not need to be highly pure or 
of closely controlled composition. 

There is no marked effect on the growth rate if 
the crystal is moved in the solution, on account of 
the high viscosity; halts in the motion, if not very 
prolonged, also cause no inhomogeneity (contrast 
ammonium phosphate and potassium tartrate). 

Defects (cloudiness, cracks) arising during 
growth are very rare in this method of growing 
EDT crystals, though cracks sometimes arise to- 
ward the end, especially while the crystals are be- 
ing removed. These cracks take a highly charac- 
teristic form, being confined entirely to the sur- 
face; they arise mainly on the {001} cleavage. 

Some defects are caused by incorporation of de- 
composition products; the crystals become yellow- 
ish if the solution is so impure as to be brownish 
and largely opaque. We have no evidence that the 
properties of the crystal are seriously altered, nor 
is it clear what products are incorporated. The de- 
composition products tend to depress the growth 
rate. It is clearly undesirable to allow excessive 
contamination by these products, but it is extreme- 
ly simple to check on this from the color of the so- 
lution. Filtration (see above) removes the impuri- 
ties. 

Cylindrical vessels (such as are used for am- 
monium phosphates) may be employed in large- 
scale production; these may be of stainless steel 
and of volume up to 150 liters, with viewing parts, 
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fittings for thermometers and temperature regula- 
tors, crystal holder, and close-fitting covers. It is 
not essential to use a reciprocating crystal holder; 
continuous rotation can be employed. The holder 
may be a disk of plastic or stainless steel, with 
several seeds mounted near the edge. Several disks 
may be mounted on the same shaft if the crystal- 
lizer is of sufficient volume. 

Heads of crystals give the best seeds for rou- 
tine production. The lower end of the crystal (that 
containing the head) then usually contains defects 
arising from regeneration, so this part must be dis- 
carded. 

Tests were done to devise a method of grow- 
ing EDT crystals free from these defects; a spe- 
cial orientation of the seed may be used together 
with a particular form of mounting in the holder. 
These tests showed that large EDT crystals canbe 
grown almost free from these defects arising from 
regeneration. In addition, the new type of seed re- 
quires a much smaller volume of crystal. 

New seeds are cut from a head formed by in- 
tersection of the faces of the {110} dihedron, which 
is sliced along a series of diagonal planes parallel 
to the length of the crystal (Fig.3). This gives sev- 
eral triangular plates, each of which can act as the 
seed for a crystal of the same width as the initial 
ohe. The plates obtained from the sides of the head 
are small and are of no real value for growing fur- 
ther crystals. 

Defects arising from regeneration are mini- 
mized by reducing the surface area of the seed; 
point seeds would be the best on this basis, but they 
do not allow one to grow broad crystals. The fol- 
lowing method of holding the plate seeds provides 
essentially growth of the pinacoid faces from a re- 
stricted area. The plate is placed in a deep slot 
closed at the bottom with the vertex upward and 
slightly projecting from the slot (Fig.4). The crys- 
tal grows outward from the vertex, and the {110} 
faces that gradually grow out from the slot con- 
stantly stimulate the growth of the pinacoid faces, 
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Fig. 5. EDT crystal grown on the vertex of a seed 
plate (early stages of growth). 


which continues until the base of the plate is in- 
scribed within the perimeter of the crystal. The 
cross section of the finished crystal is thus gov- 
erned by that of the base. 


This method may be used with other substances 
that show a pronounced tendency to grow in one di- 
rection only. Figure 5 shows an EDT crystal in the 
early stages of growth from a seed plate. 


High supersaturations should not be produced 
until the crystal growing from the seed has reached 
the required size, This implies special control of 
the early part of the temperature-reduction curve, 
if this is the method employed to cause growth, 


EDT decomposes rapidly above 50°C, while the 
monohydrate may be deposited below 40°C, so we 
made some tests on crystal growth at a fixed tem- 
perature of about 45°; the high solubility of EDT 
makes it quite feasible to produce growth by solv- 
ent evaporation. 

Solvent removal was effected by a very simple 
method (layers of filter paper inserted under the 
edges of the cover). Water condensing on the cover 
is carried away by the capillary action of the paper. 
The volume of solution in the vessel is kept constant 
by adding EDT solution saturated at a lower tem- 
perature. 


EDT crystals were also grown by reducing the 
temperature of the supersaturated solution. The 
initial supercooling was 8 to 12°C. The linear tem- 
perature dependence of the solubility makes it very 
simple to draw up a temperature program, espe- 
cially since the crystal grows in one direction only. 
The only quantities that need to be known are the 


Pha Gs CRO ZDNMAROM 


volume of solution, number, and cross section of 
the crystals, and the permissible growth rate. The 
last was taken as 5 mm/day in our calculations. A 
fixed rate of reduction is used, apart from the ini- 
tial stage, where a reduced rate is used until the 
seeds have grown out to the required cross section. 
In addition, the fixed rate is exceeded somewhat to- 
ward the end of a run in order to compensate for 
the slight tendency to lateral growth. 

We have also used other seed plates for some 
time. These were cut parallel to the {110} dihed 
ron face as slices 1.5-2.5 mm thick from the en- 
tire crystal (Fig.6). These seeds are held on flat 
inclined holders mounted on a common base. These 
seeds regenerate satisfactorily and give large uni- 
form crystals, but they have the disadvantage that 
the end of the crystal is oblique,i.e., parallel to 
the {110} face. On the other hand, they are better 
than plates cut from crystal heads, because there 
is considerably less waste from cutting. In addi- 
tion, the complete absence of defects increases the 
proportion of the crystal that can be utilized. A 
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Fig. 8. EDT crystal (74 natural size). 


slight difficulty arises over regeneration for seeds 
of very large area. Figure 7 shows a crystal grown 
from such a plate; in this case, the temperature 
may be reduced uniformly, as when complete heads 
are used. 

Although it is possible in principle for the 
monohydrate to be deposited if the temperature is 
taken below 40°C (which could bring the run to an 
end), we have encountered no case in which this oc- 
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curred, even when the temperature was reduced to 
25°C. There may be other reasons for preferring 
the isothermal method (operated above 40°C). 

EDT crystals are very easily grown by simple 
techniques to be entirely free from defects. Figure 
8 shows an EDT crystal at slightly less than natural 
size. We have produced crystals of slightly over 1 
kg, but there appear to be no obstacles to the pro- 
duction of larger ones. 

We are indebted to Academician A.V.Shubnikov 
for comments on the draft of this paper. 
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A new type of nomogram was proposed in [1] 
for calculation of monoclinic and hexagonal struc- 
ture factors, as these factors were not reducible 
to one-term derivatives of the cos 27hx cos 27ky- 
cos27Tlz type. The new method was proved to be 
very useful in a series of recent structural studies 
[2-5], carried out with the participation of the pres- 
ent authors; in particular, it has been found worth- 
while for orthorhombic structures as well, where, 
until recently, the logarithmic nomogram method 
predominated, which was itself proposed by one of 
the present authors, and "mechanized," as de- 
scribed in [6]. 

Logarithmic nomograms are inconvenient, be- 
cause it is necessary to use the same scale forboth 
sines and cosines, and because of the effort in- 
volved in remembering the signs corresponding to 
a given quadrant. 

By using nomograms of the type described in 
[1], the need no longer arises to bring the structure 
factor formulas to the "logarithmic" form, which 
takes up considerable space in the new edition of 
the International Tables. Groups with their origin 
at a center of symmetry contain no sines at all. Let 
us examine the structure factor calculation for the 
following symmorphic groups: 


D3, = Pmmm, DY, =Cmmm, D5} = Fmmm, D3, = Immm. 
A =n {cos 2n (hz + ky + Iz) + cos 2x (he + ky + Iz) 
+ cos 2x (ha + ky + Iz) + cos 2a (hx + ky + 1z)}. 


As is usual in structure analysis, we begin 
with the structure factors for the hk0 reflections, 
i.e., with the expression A = n{cos 27 (hx + ky) 

+ cos 27 (hx — ky)}. Let x = 0.35and y= 0.21. On 

millimeter graph paper (Fig.1) we chose a central 
bold line, on which we place a cosine scale (twice), 
at a scale of 1 = 20 cm. On one of the bold lines to 


the left of the central scale we plot the values of hx 
onal=40cm scale, writing h=1,2,... at the 
ends of the segments. If hx is more than one (40 
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cm), the corresponding fractional part is plotted. 
On the symmetrically opposite bold line (on the 
right), the values of ky are put but, in this case, 
twice; once counting from the bottom of the 40-cm 
scale, and once from the top (again, only fractional 
parts). In the second case, the figures are marked 
with minus signs (1;°2,°3)...). 

A straight line (drawn on tracing paper), which 
connects a given h on the left-hand scale with a 
given k on the right-hand scale, cuts the central 
line at the value of cos 27 (hx +ky). The straight 
line h to k gives the second term, cos 27 (hx — ky). 

If a group is not symmorphic, then, depending 
on whether there is an even or odd number of halves 
involved in the argument of the second term in the 
structure factor expression, the corresponding 
cosine will be either added or subtracted. Thus, 
for the group Dy = Pnma, and for the hk0 reflec- 
tions, the structure factor will be 


A = 4cos 2e (hx + ky) + 40s 2 (hx + hy) 
for even k, and 
A= 4cos 25 (hx + ky) — 4c08 20 (he + ky) 


for odd k (h is always even for Deis according to 
systematic absence rules). 

The structure factors are usually calculated 
according to the layer lines. For the /-th line, the 
corresponding expression (starting off with the 
group Dé, = Pmmm), is: 


2 {cos 2x (hx + ky + Iz) + cos dz (ha + ky + Iz) 
+ cos 2x (ha + ky + Iz) + cos 2x (ha + ky + I2)}, 
which is rewritten in the form 


2 {cos 2 (lz + hu + ky) + cos 2n (lz + haz + ky) 


+ cos 2n (lz + hae + ky) + cos 2x (/z + hx + ky). 


On the left-hand scale we fix the value of lz, 
and we begin to plot the values of hx only from this 
point, with, as before, the counting starting at the 
bottom when the 40-cm scale is run over. From 
the same Jz point, we plot the values of hx along 
the opposite direction, marking them Loe 

On the right-hand scale, in the same order as 
previously, we mark off k and k. A line from h to 
k cuts the central scale at the value of the first 
cosine: h to k, h to k, and h to k give the second, 
third, and fourth cosines. With nonsymmorphic 
groups, the second, third, and fourth terms may 
have either plus or minus signs, depending on 
whether there is an odd or even number of transla- 
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tional halves in the argument, for example, for the 


group Dp Pnma, there are four possible cases: 

hep l = 2n peal scales Oe 

oa ==) 

h-+ ; == on ale) ee Z 

ee == 

ees pcr = a 

| ir 
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For the group DJ, = Pnma, with the general 
structure factor expression 


A = 2 {cos 2e (hx + ky + Iz) + cos 2s (hr + ky + Iz) 


+ cos 2 (he thy +le+ “ =) 
| 
J 


+ cos Qn (he + ky + le aL _ ‘) 


there are, in all, two possible cases: 


h+l=2n a 


h+—l=2n+41 4+ 4+ — —. 

If the group is not centrosymmetric, it is 
necessary to separately calculate the cosine part 
(A) and the sine part (B) of the structure factor.The 
first part is calculated by the method already de- 
scribed, and the second part by raising the central 
scale a distance of 5cm. It is convenient to have 
the central scale movable (compare [1]), not only 
for changing frem cosines to sines, but also so that, 
instead of marking the value of Jz on the left-hand 
scale, the central scale can be simply moved by the 
distance /z, working to a scale of half the size. 

For the group DS = P222, 


A = 0s 2x (ha + ky + Iz) + cos Qn (he + ky + Iz) 
+ cos 2m (ha + ky + Iz) + cos 25 (hc + ky + Iz), 


B = sin 2n (ha + ky + Iz) + sin Qe (he + ky + lz) 
+ sin Qn (he ky + Iz) +4 sin 2h (he ky + lz), 


which can be rewritten as 


= (iz + hx + ky) +008 2n (Iz + he + hy) 
ees lene + ky) 4 +. cos 27 (lz + hz + ky), 

Iz + hz + ky) — sin 2x (lz — ha + ky) 
n(lz-+ he + hy) + sin 2c (iz + hz + ky), 


i.e., the cosine Sequence is++++;the sine Sequence, 
pif 
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For the group Pca2,, with the origin on the 2, 
axis, the general structure factor expression is 


A = cos an (he + ky + 1z) + cos 2e (ha + hy + lz + 


h--l 
3 


+ cos 2m (ha + hy + let 2) + cos 2x (ha + ky + le + 5), 


B= sin 2 (he + ky + lz) + sin 25 (ha: + ky + la+ 


} 


+ sin 2x (he + ky tlg+t =) + sin 2a (hic + ky +lz + 


and we have the four sign sequences: 


(PND, 1=2n “A+ 
Bw 

ss Pp l=2n- lA 
file 

Sie heap j + = 
B= 


is Pte. Pe i Seah Ay ae 
B A- 


It is easy to make up similar sequences for all 


the 230 space groups.! 


za 


“+ = 


ae ge 


+4 


ht ‘) 


l 


a 
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Some very full x-ray studies [1] of alloys from 
the Bi—Sb system have shown that these components 
form a continuous series of solid solutions, withan 
unlimited range of solubilities. It appeared of inter- 
est to investigate the thin layers obtained by simul- 
taneous condensation of the vapors of these twocom- 
ponents in a vacuum, since these layers actually 
demonstrate a continuous change in the concentra- 
tion of each of the metals. 

The films were condensed either on celluloid 
or on a fresh sliver of rocksalt crystal. The tem- 
perature of condensation or subsequent annealing 
was varied from experiment to experiment, from 
room temperature up to about 100°C for the cellu- 
loid, and up to 300-350°C for the rocksalt. As the 
condensation and annealing temperatures were in- 
creased, the degree of orientation of the film crys- 
tallites was improved. Ata 50:50 composition, 
several diffuse reflections were observed. 

Over the whole range of concentrations studied 
by the electron diffraction method, an identical 
structure was observed, corresponding to that of 
antimony or bismuth, with continuously variable 
unit cell constants; a was equal to 4.53 + 0.01 A 
(on the Bi side) and 4.29 + 0.01 A , and c was equal 
to between 11.88 + 0.02A and 11.24 + 0.02 A. 

The films prepared on celluloid had their crys- 
tallites oriented with their basal faces parallel to 
the substrate and in a random azimuthal arrange- 
ment. By annealing at 100°C it was possible to ob- 
tain good orientation with any composition. 

The films prepared on NaCl, particularly those 
subjected to annealing after condensation, did show 
azimuthal orientation. From each of the hki0 ring 
cleavages, 12 point reflections were observed, in- 


dicating two equally probable dispositions; any vec- 
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tor of the [100], [110], [010] type in the metal lat- 
tice could be parallel to the [110] and [110] direc- 
tions in NaCl. 

A series of specimens on NaCl crystals, rep- 
resenting the whole range of concentrations, was 
annealed at t = 200°C for 15 h. No changes were ob- 
served in the diffraction patterns as a result of this 
annealing. 

The formation of well-oriented films, even when 
condensing on substrates at room temperature, 
shows that the Bi and Sb atoms are quite mobile on 
these substrates under the experimental conditions 
used. We therefore assumed that the fact that 15 h 
of annealing at 200°C had no effect on the diffrac- 
tion pattern was proof that these alloys could not 
attain an ordered arrangement of their Bi and Sb 
atoms. 

Experiments were carried out on the oxidation 
of films of the pure components and of their alloys. 
An Sb film on rocksalt oxidized on heating to t 
= 180°C in air at a pressure of 10°!-107 mm Hg. 
In this case, senarmontite, Sb,Og, was invariably 
formed, this having a cubic lattice with a = 11.14 
A. In all cases, the senarmontite crystallites were 
oriented with their (111) planes parallel to the sub- 
strate. If the antimony specimen was azimuthally 
oriented, then the Sb,O, crystals repeated this 
orientation. The relative orientations of the Sb and 
its oxide are a reflection of the well-known rela- 
tionship between cubic and hexagonal lattices; 

[111 ]Joub || [001]hex and [110 ]eub ll [100 Ihex- 

When Bi films were oxidized, or Bi vapor was 
condensed under preliminary vacuum conditions (at 
200-250°), oxide phases were apparently observed; 
these had cubic lattices with a, = 4.65 A (1) and a> 
= 5.49 A (II). The crystals of these phases were 
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Fig. 1. Electron-diffraction pattern from a Bi film Fig. 2, Electron-diffraction pattern from a film of 
containing a cubic oxide phase with a = 4.65 A, azimuthally oriented Bi, containing a cubic oxide 
phase with a = 5.49 and NaCl. 


regularly oriented (Figs. 1 and 2). 

In Bi specimens without azimuthal orientation, 
the crystallites of (I) were oriented with their cube 
faces parallel to the substrate, forming a mosaic 


A, is evidently identical to a phase described in [2]. 
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The question of the manner in which the ability 
of crystals to form regular overgrowths depends on 
the relative difference in parameters of the over- 
growing flat networks has been considered by 
Yanulov [2]. He quotes some old experiments by 
Frankenheim [4] on the preparation of growths of 
various halide salts (NH,I, KI, KBr, KCl) on the 
(001) cleavage planes of biotite and muscovite. 

On the basis of experimental studies, Yanulov 
came to the conclusion that the relative number of 
oriented crystals forming on the mica was a simple 
linear function of the relative difference between 
the parameters of the flat networks growing to- 
gether (Fig.1). The substances were crystallized 
by evaporating drops of solution on freshly cleaved 
flakes of the mica. Excess solution was removed 
with a filter paper. In Yanulov's opinion, this ex- 
cluded the possibility of nonoriented growth being 
obtained through precipitation of crystals from the 
upper layers of the drop. 

The present study was to investigate the possi- 
bility of deriving quantitative relationships between 
the overgrowth behavior of KI and KBr on muscovite 
and the degree of similarity of their parameters. 

The process of crystallization of KI and KBr 
from solution onto muscovite can be represented as 
follows. When the drop of solution is placed on the 
surface of the crystal, evaporation of the liquid 
takes place, leading to the formation of two kinds 
of crystals. One of these arises in the immediate 
vicinity of the substrate, and is oriented!; the other 
is formed outside the zone of influence of the crys- 
tal substrate surface; these latter very fine crys- 
tals, when "their weight is sufficient, fall onto the 
surface, where they may lay in any position" [5]. 
This type of crystal generates nonoriented over- 
growths. 

Yanulov's conclusion that there was a linear re- 
lationship between the formation of regular over- 


growths and the degree of similarity of the crystal 
parameters does not seem to be very likely to us, 
because, in his experiments, he did not completely 
eliminate the causes of crystals precipitating from 
the upper layers of the liquid. 

In order to avoid crystals precipitating from 
the upper layers of the drop, we used a vertical, 
flat crystallization vessel, with the mica plate 
placed in it in a vertical position also. On this, 
under the action of surface forces, only oriented 
crystals arose. The observations were made using 
a microscope. 

The solution with the mica plate immersed in 
it was either left several days to stand or was 
cooled rapidly with dry ice. In the latter case, 
growths were obtained after 10-15 min. The experi- 
ments were carried out in the KI—muscovite and 
KBr—muscovite systems. During the crystalliza- 
tion of KBr on muscovite, when the relative differ- 
ence between the network parameters was 11.4%, 


3 6 9 1215 18A% 


Fig. 1. Relationship between the 
difference in parameters of inter - 
growing flat networks (A %) and 
the percentage orientation of 
crystals (%), 1) Crystallization 
on muscovite; 2) crystallization 
on biotite (after K.P. Yanulov). 


‘Bradley [3] has established that the zone of surface forces of mica 
extends to 10°* cm when NHI is crystallized onto its surface. 
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Fig. 2, Oriented crystallization of KBr on 
muscovite. Magnification 120x, 
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Fig. 3. Graph of Inw against for 


co 


the KI—muscovite system. 


Yanulov obtained only 39% of oriented crystals (Fig. 
1), while in this case we found complete orientation 
of all the 500 KBr crystals studied (Fig. 2). 

Thus, where in crystals the relative difference 
between their parameters and that of muscovite is 
increased from 11.4% (crystallization of KBr on 
muscovite with full orientation) to the values per- 
mitted by Royer's limit rule (about 15%, for example 
as in KCl), there must be an abrupt change in their 
ability to form regular overgrowths. 

In our view, therefore, it would be wrong to 
assume that there is a gradual change in percent- 
age orientation, depending on the relative differ- 
ence between the parameters of the flat networks 
growing together. Our experiments indicate, rather, 
that the proportion of oriented crystals varies in 
stepwise fashion. 

As the quantity which varies according to the 
difference in network parameters, we can examine 
the value of the "limiting" supercooling [1], i.e., 
that value of solution supersaturation at which pre- 
cipitation of crystals onto the support occurs. The 
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TABLE 1 


Time of appearance of crys- 
tallization centers (av. of 
many measurements),* sec 


Temp. at start of 
crystal separation, 
eC 


Saturation 
temp., 
[G 


10.5 150 
20 10.0 100 
20 9.0 50 
20 8.5 10 


* The time of appearance of crystallization centers is the re- 
ciprocal of the probability of their formation, T = 1/w [1]. 


supersaturation of the solution may be expressed by 
the quantity (In (c/co), where c is the equilibrium 
concentration, and c,, is the concentration atwhich 
precipitation of crystals begins. 

In this connection, a relationship of consider- 
able interest is the straight-line graph obtained by 
plotting the logarithm of the probability of forma- 
tion of crystallization centers (Inw) against a quan- 
tity inversely proportional to the square of the 
supersaturation, 


This relationship was first established by 
Ovsienko [1] for the crystallization of aqueous NaBr 
and NaCl solutions on the surface of galena. Our ex- 
periments on the KI—muscovite and KBr —muscov- 
ite systems confirm this relationship. 

Part of the experimental data on the KI— mus- 
covite system is given in Table 1. 

From the results obtained, a graph was con- 
structed (Fig.3). The results obtained by Ovsienko 
and by us allow conclusions to be drawn as to the 
change in the value of the "limiting" supercooling 
with change in the difference between the param- 
eters of the flat networks. 

Thus, in the case of the NaCl—PbS system (A 
= 5.7%), precipitation of crystals of the dissolved 
substance takes place on the PbS at a degree of 
supersaturation one quarter of that pertaining be- 
fore the introduction of the PbS crystal, while with 
NaBr solutions (A = 0.2%), the supersaturation is 
reduced by a factor of 10 [1]. 

As the probability of appearance of the crystals 
is related to the supersaturation of the solution, it 
might be considered that changes in the parameters 
of the intergrowing flat networks would have a di- 
rect effect on the probability of appearance of crys- 
tal growth nuclei. Table 2 shows how the time of 
appearance of crystallization centers alters depend- 
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TABLE 2 


Saturation | Temp. at start of Time of appearance 
System A % temp., | crystalseparation, | of crystallization centers, 
°G 1G sec 
KI— muscovite 3.81 20 10 100 
KBr—muscovite| 11.4 20 10 Crystals: did not appear 
KI— muscovite 3.81 40 20 Crystals appeared moment- 
arily 
KBr- muscovite| 11.4 40 20 30 


ing on the relative difference between the network 
parameters (A %). 


Conclusions 


1. The percentage of oriented crystals de- 
pends, above all, on the experimental conditions. 
It is possible to obtain complete orientation of one 
substance growing on another even when a con- 
siderable difference exists between the parameters 
of their flat networks. 

2. The percentage orientation does not change 
gradually, but abruptly, close to a limiting value. 

3. The drop in the value of the "limiting" 
supercooling is due to the reduction in the relative 
difference between the parameters of the two flat 
networks concerned, which in its turn affects the 
probability of the appearance of crystal growth 
nuclei. 


The author is grateful to V.A. Frank- 
Kamenetskii for his guidance in this investigation. 
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It is well known that a- and @-quartz, in the 
immediate vicinity of the phase transition point, 
have a rapid temperature variation of the index of 
refraction. Consequently, if the phase transition is 
not completed simultaneously throughout the entire 
volume of the crystal, but is gradually propagated 
from several centers of formation of the new phase, 
under certain conditions it is evidently possible to 
observe an intense and almost regular reflection of 
light at the a- and 8-phase boundaries, and a re- 
distribution of light in regions adjacent to these 
boundaries. 

This circumstance has been used for carrying 
out the following simple method of observing the 
entire process of phase transition in large single 
crystals of quartz. 

The quartz specimens, in the form of cubes 
with the linear dimensions 20-25 mm, were placed 
in a parallel beam of light, as shown in Fig.1. In 
Fig. 1, A is a light source, Ba condenser, Ca 
circular point diaphragm, Da collimator lens, Ea 
cylindrical vacuum furnace, F the specimen, and 
K a screen or photographic plate. 

With such an arrangement, the shadow image 
of the specimen, shown in Fig. 2a, was obtained on 
the screen. 

The furnace construction ensured uniformity of 
the temperature field, necessary for accomplishing 
a phase transition in large quartz specimens with- 
out damaging them. Slight temperature gradients of 


the order of 0.1-0.2°C/cm in the cross section of 
the furnace enabled phase transition fronts to be ob- 
tained parallel to the light rays. 

The apparatus was used for following the phase 
transition processes in a free crystal, in a crystal 
exposed to the action of a strong dc field, and also 
in a high-frequency field. The entire phase transi- 
tion was photographed on a motion picture film. 

Figure 2b shows the occurrence of centers of 
formation of a-quartz in a single crystal of high- 
temperature 8-quartz. When a center of the new 
phase was formed, temperature redistribution oc- 
curred in its immediate vicinity and large gradients 
of the indices of refraction were produced, result- 
ing in distortion of the outline of the shadow of the 
Specimen, a redistribution of the light flux, and the 
appearance of a dark area around the place of forma- 
tion of the new phase. The center of formation of 
the new phase itself was revealed in the photograph 
(Fig. 2b) in the form of a bright spot. The centers 
produced on the lateral boundaries were concealed 
by the mounting and are not visible. Their presence 
may be judged only by the distortions of the contour 
and the occurrence of a dark region. 

In Figs. 2c-2l are shown the photographs of dif- 
ferent successive stages of development and move- 
ment of the phase transition fronts. The contours of 
the part of the front which is almost parallel to the 
light rays are revealed in the form of a bright line 
in the photograph. 


BRIEF COMMUNICATIONS 189 


Bags 2. 


Observation of the character of movement of 
the front showed that the sites of origin of the cen- 
ters of formation of the new phase were determined 
entirely by the temperature field of the furnace. 
The time required for the complete accomplish- 
ment of the phase transition was determined by the 
rate of supply of heat to the specimen, and could 


be indefinitely increased by the choice of suitable 
heat-supply conditions. In the case of slow conduct 
of the phase transition, the specimen retained its 
single-crystal structure even in repeated transi- 
tions, despite the fact that the formation of the new 
phase always commenced simultaneously at a num- 
ber of points in the specimen. In the case of rapid 
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transitions, of the order of 2-3 min, the specimens 
sometimes twinned. On approaching the phase- 
transition temperature, the boundaries of the twins 
were clearly observed in the form of thin bright 
lines. Distortion of the front of formation of the 
new phase occurred at the boundary of twins, and 
the passage of the front across a twin boundary ap- 
peared to be difficult. The velocities of the front 
were different in different stages. A sharp in- 
crease in velocity was observed at moments pre- 
ceding the joining together of a number of fronts. 
In the case of sufficiently uniform external 
temperature field and slow heat supply, the form 
and direction of movement of the phase-transition 
fronts were determined basically by the character 
of anisotropy in thermal conductivity of the quartz. 
The patterns of development of the phase transition 
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were different when observed along different crys- 
tallographic axes. For each of the axes, the char- 

acter of these patterns remained the same in all the 
specimens. 

If no twinning of the specimen occurred in the 
process of the a— 8 transition, the 8 —a@ transition 
repeated strictly all the stages of the first transi- 
tion in the reverse order. 

It should be mentioned that a number of rapid- 
ly moving dark bands were always formed in front 
of the phase-transition front. 

The method described enables one to follow 
quite reliably not only the phase-transition process 
itself, but also the character of the variation of the 
temperature field in the specimen during phase 
transition, and it may be used in work on the exact 
determination of the phase-transition temperature. 
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In a recently published paper by V.I. Mikheev 
[L], it was asserted that we, as it were, "had ex- 
pressed the opinion that the phenomenon of the dis- 
persion of the optical indicatrix... was incom - 
patible with the homology of crystals." In fact, in 
our earlier papers on homology [2,3,4], we denied 
the possibility of a definite connection between the 
homology of the external form of crystals and that 
of- the indicatrices of their physical properties, in 
particular, their optical indicatrices. Considered 
individually, the external form and optical indi- 
catrices can and do possess properties of homology, 
but not a single operation of form homology (apart 
from those simultaneously constituting symmetry 
operations) enters into the group of operations for 
the homology of the optical indicatrices. 

As we see it, our actual opinions and those 
ascribed to us do not quite coincide. This, however, 
is not the point. 

In the paper cited, V.I. Mikheev tries to re- 
fute our assertion in the following way. He takes a 
polyhedron of the tetragonal system with an optical 
indicatrix for two wavelengths inscribed in it (the 
indicatrix is laid out quite incorrectly, but this is 
not important in the present case). He then sub- 
jects this model of a tetragonal crystal to extension 
along some direction and finds, as we should ex- 
pect, a model of a monoclinic crystal possessing 
the properties of homology which he (the author) 
wished to obtain. In this discussion, however, V.I. 
Mikheev makes the great mistake of arbitrarily as- 
suming that the particular deformation (in this case, 
extension) whereby the tetragonal form of the crys- 
tal may be imagined to be converted into monoclinic 
and the deformations whereby the optical indi- 
catrices of the tetragonal crystal are simultaneous- 
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ly transformed into indicatrices of a monoclinic 
crystal are one and the same. It is in fact quite 
well known that these deformations cannot be iden- 
tical. 


If we take V.I. Mikheev's point of view, then 
we deduce that, on subjecting crystals of the cubic 
system to uniform elastic compression along one 
of the normals to the cube and octahedral faces, the 
crystals must necessarily acquire optical proper- 
ties of a uniaxial-negative crystal, since the spher- 
ical optical indicatrix of such crystals should under 
these circumstances always transform into an 
oblate ellipsoid of revolution. In fact, in suchcases, 
as shown hy direct experiment, very different re- 
sults may be obtained. Rocksalt, for example, on 
compression along the normal to the (100), be- 
comes optically uniaxial-negative, while sylvine 
under the same conditions becomes uniaxial-posi- 
tive. On compression of fluorite along the normal 
to the (111), the crystal becomes optically uniaxial- 
positive, while under the same conditions sylvine 
becomes uniaxial-negative ([5], p. 481). 


Analogous phenomena are found for crystals of 
all other systems, not only for mechanical, but al- 
so thermal deformation. If we are to be guided by 
V.I. Mikheev's views, then for all crystals of the 
cubic system the spherical optical indicatrix should 
expand uniformly in all directions on thermal ex- 
pansion, i.e., the refractive index should increase. 
Actually, the thermal expansion of such crystals as 
rocksalt, sylvine, fluorite, and aluminum —potas- 
sium alum leads to a fall in refractive index, while 
in other crystals (diamond, zinc blende) under the 
same conditions the refractive index increases [5], 
p.451). 
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The lack of agreement between the thermal de- 
formation of a crystal and that of its optical indi- 
catrices is found not only for cubic crystals but for 
crystals of all systems. It is well known, for 
example, that, on heating, a calcite crystal ex- 
pands along the optical axis and contracts in a per- 
pendicular direction, whereas the optical indicatrix 
under the same conditions expands in both these di- 
rections ((5], p. 452). 

V.I. Mikheev asserts that on elongation of a 
crystal its optical indicatrix should also expand, 
but then one can assert with equal justification that 
the Fresnel ellipsoid should do likewise. This, how- 
ever, would contradict common sense, since the 
elongated indicatrix should correspond to an oblate 
Fresnel ellipsoid. 

It is perfectly clear from our discussion that 
no definite connection exists or can exist between 
the homology of the external form of crystals and 
that of their optical properties. 

In V.I. Mikheev's article one may also find a 
whole series of "errors" naively imposed upon us. 
For example, Mikheev asserts that we ascribe 
some "spherical indicatrix of electric charge dis- 
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tribution" to a monoclinic crystal. In fact, in our 
earlier paper [3] we were considering not the 
spherical indicatrix, but a sphere ground out of a 
crystal. Our other "errors" are of the same nature. 
There is no need to discuss them. 


LITERATURE CITED 


1. V.I. Mikheev, "Possible connection between 
the dispersion of the disposition of optical in- 
dicatrices and the homology of crystals," Zap. 
Min. Obshch., Part 84, 170 (1955). 

2. A.V. Shubnikov, "On the so-called homology 
of crystals," DAN SSSR, 838, 453 (1953). 

3. A.V. Shubnikov, "What the homology of crys- 
tals is," Tr. Inst. Kristallogr., 9, 35 (1954). 

4, A.V. Shubnikov, "On the absence of a definite 
connection between the homology of the ex- 
ternal form of crystals and that of the indi- 
catrices of their physical properties," Zap. 
Min. Obshch., Part 84, 121 (1955). 

5. F. Pockels, Lehrbuch der Kristalloptik, pp. 
451, 452, 481. 


SOVIET PHYSICS — CRYSTALLOGRAPH 


og 


VOR. tl, NOV 2 MAR.-APR., 1956 


REVIEW 


APPLICATION OF COMPUTERS TO X-RAY 


STRUCTURAL CALCULATIONS 
M. A. Porai-Koshits 


N.S. Kurnakov Institute of General and Inorganic Chemistry 
Translated from Kristallografiya, Vol. 1, No. 2, 


pp. 248-263, March-April, 1956 


Original article submitted November 9, 1955 


I. The mechanization of computing processes 
in x-ray structural analysis is one of the most 
pressing problems in x-ray instrument making. 
Naturally, the main subjects for mechanization are 
those calculations which are most laborious to exe- 
cute, most urgently required, and most standardin 
form. In the field of x-ray structural analysis, this 
category, above all, includes the solution of prob- 
lems arising in the study of atomic structure in 
crystals of complex composition. 

Three main factors govern the difficulty of 
structural investigations: (1) the so-called "initial- 
phase problem" (the absence of any direct experi- 
mental data on the phases of diffracted rays); (2) 
the duration of the experiment and the preliminary 
analysis of experimental data; (3) the difficulty of 
the computing operations involved in the determina- 
tion of atomic coordinates. 

The statistical method of determining the co- 
ordinates of atoms, developed and successfully ap- 
plied in recent years (especially in the Soviet 
Union), solves the first (and most basic) of these 
difficulties, at least as far as centrosymmetrical 
crystals are concerned. The construction of dif- 
fraction cameras with ionization recording of the 
rays and automatic adjustment of the crystal and 
receiver makes the actual experiment much faster. 
Mechanization of the main computing processes 
solves the third problem. Taking all together, we 
may well expect soon to witness a radical change 
both in the speed and volume of structural investi- 
gations, and hence also in the very position of x- 
ray structural analysis in relation to fields of sci- 
ence constituting the "clientele" of structural 
studies, chemistry, geology, physics, biology, etc. 
The hundredfold increase in the speed of calcula- 
tions achieved with the help of electronic computers 
constitutes a powerful stimulus to reinforced ef- 
forts toward the perfection of experiment and the 
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development of methods of structural investigation. 
It is thus clear that the maximum acceleration of 
computing operations in the analysis of crystals, 
complex in structure,is important in principle as 
well as technologically. 

II. The main computing problems in the field 
of atomic structural analysis are the following. 

1. Finding the distribution of electron density 
over the cell (or the interatomic function in inter- 
atomic vector space), i.e., Summing the two- and 
three-dimensional Fourier series, in general of the 


type 


0 (xyz) = va {F'4 (hkl) cos 2x (hx + ky + Iz) 
hkl 
+ Fp (hkl) sin 25 (ha+hky+1z)}. 
2. Calculating the structure amplitudes (or 


unitary structure amplitudes) of the diffracted rays 
from the coordinates of the atoms and the atomic 
amplitudes (per-unit coefficients), i.e., calculating 
a sum of type 


NN 
OU) = Dy {; cos 2x (he; + ky; + lz, 


eh 


), 


y 


Ip (h kl) —- ye Wi sin Wigs (hz; _ hy; ates ey. 


j=1 
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3. Refining the coordinates of the atoms by 
successive approximations to the electron-density 
distribution and the initial phases of the reflections. 
The most convenient form of this refinement is the 
method of differential synthesis, at each stage of 
approximation requiring the solution of the system 
of equations: 


Ann Ax + Anz Ay + Apr Az+ An = 0, 
Anz [ep a= Arr Ay + Ax Az aie Ax, =) 


Ap; Ax +- Ay, Ay + Ay, Az + A; = 0, 
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where Ax, Ay, Az are unknown corrections to the 
coordinates, andA}y, App, Anz etc., are Fourier 
expansions of the corresponding derivatives of elec- 
tron density at a point corresponding to the original 
coordinates of the atoms: 


at ae 
Pe Ve ei UyeG ah du* eeuyey 


o2 
Ani = ( - ) » etc. 


Ox Oy FU j=; 


4, Determining the signs of the structure amp- 
litudes in a centrosymmetric crystal by the sta- 
tistical method. In the simplest form, this problem 
reduces to the setting up of all possible combina- 
tions of reflections hkl, h'k'' andhth', k+k’, 

1 +1' belonging to a group of strong reflections, 
then to a wider group of strong and medium, and, 
finally, to all reflections including the weak (but of 
course existing) ones. 

In its more complex setting, this problem in- 
cludes the calculation of the probability of satisfy- 
ing the sign relationship s(hk/) +: s(h'k'l') = s(h +h} 
k+k', 1 +1'), where s is a symbol for the sign of 
the structure amplitude, i.e., +1 or —1. 

5. Calculating the quadratic form, i.e., find- 
ing 1/d? or (sin ¥/A)* from the lattice parameters 
and diffraction indices. In the general case, the 
problem reduces to solving the equation 


Ont = AW® + Bh? + CP + Dhh + Lhl + FRI 


for all necessary combinations of hkl. 

In recent years, a great deal of attention has 
been paid abroad to mechanization of the main com- 
puting operations of structural analysis. Many ma- 
chines and mechanisms have been developed for 
solving the first three of the problems just listed. 
The largest laboratories have in some way or other 


been provided with the necessary computing devices. 


In the Soviet Union, also, considerable work 
has been done on constructing computers for struc- 
tural analysis. Ata conference on the use of com- 
puting techniques called by the Commission for X- 
Ray Diffraction: at the Institute of Crystallography, 
Academy of Sciences of the USSR in Moscow, on 
April 14-15, 1955, six contributions were devoted 
to the construction of computers, either in service 
or being developed. 

III. Before passing on to a description of va- 
rious types of machines, it is useful to give some 
classification as regards their main features. 

All machines may be roughly divided into two 
groups: on the one hand, none-too-complex spe- 
cialized machines intended for use universally in 
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x-ray and crystallographic laboratories (these we 
shall subsequently call "laboratory" machines) and, 
on the other hand, high-speed (specialized or uni- 
versal) machines, automating the computing pro- 
cesses to the highest possible degree, and intended 
mainly for large computing centers. 

It must be emphasized that machines of the 
first type will answer their purpose only on satis- 
fying certain conditions as regards dimensions, 
cost, and simplicity of servicing, and machines of 
the second type are conditioned by the accessibility 
of the computing center for systematic and "inertia- 
free" fulfillment of the demands of the structural 
laboratories. 

In principle, the operation of all computing 
machines may be divided into two groups: (1) con- 
tinuous-action machines or analog computers, in 
which the computing operations are imitated bycer- 
tain physical processes, and, (2), discrete-action 
machines (i.e., finally, digital computers). 

Analog machines may be used on mechanical 
motion, electrical, electromagnetic, optical, and 
other processes. Present-day digital machines 
may be separated into four main groups: purely me- 
chanical, electromagnetic pulse, electronic, and 
magnetic machines. 

Analog machines are in all cases specialized, 
intended for the solution of definite particular prob- 
lems; digital machines may be either specialized 
or universal. 

Below are listed the main types of computers 
applicable to structural analysis. The letter L in- 
dicates laboratory types and LS denotes large-scale 
machines. 


Machines acting continuously: 
1. Mechanical (L) 
2. Electrical (L) 
3. Optical: 
fly's eye (L) 
x-ray microscope (L) 
photo-adder (L) 
4, Electronic Pepinsky (LS) 


Discrete-action machines; 

1. Mechanical (L) 

2. Electromagnetic: 
specialized (L, LS) 
universal (LS) 

3. Electronic: 
specialized (LS) 
universal (LS) 

4, Magnetic: 
specialized (LS, L) 
universal (LS) 
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Continuous -acting machines are mainly of the 
laboratory type. An exception is the electronic ma- 
chine of Pepinsky. Discrete-action machines, apart 
from mechanical, are mainly large scale. Special- 
ized electromagnetic machines may be considered 
as intermediate between laboratory and large-scale 
(see below); in the field of magnetic specialized ma- 
chines we may expect the creation of laboratory- 
only desk computers (see below). 

IV. The least perfect, until recently, were 
the purely mechanical continuous-action machines. 
Quite a number of mechanisms simple in principle 
have been proposed ("hand" machines) both for cal- 
culating structure amplitudes and for summing 
series. The first group includes the machines of 
Evans and Peiser [1], Cox,! Booth [3], Bragg [4], 
Vand [5], and Clastre and Gay [6]. The main fail- 
ing in all these machines is evidently the slowness 
of the computing process, associated either with 
the fact that calculations are made for each reflec- 
tion separately and often even for each term ina 
summation one after the other, or else with the fact 
that the initial parameters are given in too complex 
or inconvenient a form. Another failing in most of 
the machines mentioned is their low accuracy. 

Mechanical computers for summing Fourier 
series have failings of a somewhat analogous kind; 
these include those of Brown and Wheeler [7], 
Robertson [8], Shilton [9], McLachlan and 
Champaygne [10], Rose [11], Vand [12], and 
Beevers and Robertson [13]. 

In all probability, however, the failings of the 
mechanical computers (constructional rather than 
matters of principle) can be overcome. Purely 
mechanical or electromechanical machines also 
have a number of advantages over any others as re- 
gards reliability of operation, compactness, and 
possibly cost. 

Deserving of attention is the general idea of 
Booth's second machine, which solves the problem 
of calculating structure amplitudes completely and 
simultaneously with respect to all the terms of the 
summation’ (described in Booth's monograph [14]), 
and Robertson's discrete-action machine for the 
summation of series which, judging from a short 
communication [15], constitutes an arithmometer, 
an analog of the ordinary method of strips. The 
machine accelerates the work by roughly a factor 
of ten, with a considerable increase in accuracy. 
The details of its construction have not yet been de- 
scribed. 

V. Among the constructionally simplest, high- 
speed, compact, and inexpensive systems for find- 
ing two-dimensional electron-density distributions 
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in two-dimensional assemblies of structure ampli- 
tudes Fy; are those based on an optical analogy. 

The optical analogy enables either the direct 
or the inverse Fourier transformation to be made 
at will (the crystal cell transforming into a distribu- 
tion of structure amplitudes in reciprocal space or, 
conversely, the reciprocal lattice into a distribu- 
tion of density in the crystal cell). 

The following systems described in the litera- 
ture (especially in a series of monographs [16,17] 
and also Eller's review [18]) are well known: (1) an 
apparatus for finding structure amplitudes based on 
the method of modeling the projection of the struc- 
ture (fly's eye); (2) a machine of the "x-ray micro- 
scope" type solving the inverse problem, i.e., 
transforming a two-dimensional reciprocal lattice 
into a density distribution; (8) a photographic de- 
vice for summing Fourier series acting on the prin- 
ciple of the superposition of sinusoidal black—white 
two-dimensional distributions on the same photo- 
graphic plate. 

Among this group of machines, we must give 
more detailed attention to that of Eller, belonging 
to the third group and manufactured in France by 
the Baudouin company [19]. The principle of action 
is shown in Fig. 1: Sis a linear light source, Aa 
screen with vertical "columns" of sinusoidally vary- 
ing blackening, and P a photographic plate. Varia- 
tion in the orientation of the sinusoidal waves rep- 
resenting the individual terms of the Fourier series 
is achieved by rotating the photographic plate in its 
own plane, variation in the frequency of the waves 
by moving the screen A toward or away from the 
photographic plate, variation in the initial phase of 
the wave by moving the screen A in a horizontal di- 
rection in its own plane. A general view of the ma- 
chine is given in Fig. 2. A cassette with photo- 


Fig. 1. Principle of construction 
of the photographic device of 
Eller. 


1Described in Booth's monograph [2]. 
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Fig. 2. Eller's machine. 


graphic plate P is rigidly connected to a disk to 
which is attached a graphical representation of the 
reciprocal lattice network; the screen S is connected 
by a system of blocks to an indicator which can 
move along a vertical rod situated in front of the re- 
ciprocal lattice network. On rotating the disk with 
the network, the orientation of the cassette changes; 
on moving the indicator along the rod, so does the 
screen—cassette distance. The adjustment of the 
apparatus consists of two operations: (1) a recip- 
rocal lattice point is brought onto the vertical 
diameter of the disk, (2) the indicator is set oppo- 
site to the point. The first operation gives the ori- 
entation of the sinusoidal distribution, and the se- 
cond its frequency. Thus, the adjustment process 
in the apparatus is extremely simple. 

The projection is obtained on a photographic 
plate 17 X 23 cm in size, and the indices of the 
terms being summed may lie between 0 and 100. 
The time needed for a complete projection with 100 
coefficients of the series is 30 min. 

The machine has very varied applications. 
Apart from the calculation of ordinary and various 
kinds of complex projections (in particular the 
minimalized projection of Burger [20]), it can be 
used to give pictures of the Fourier transformations 
of any system of points, the values of structure 
amplitudes [21], and even to solve the problem of 
the refinement of coordinates by the method of dif- 
ferential series [22]. 

The rapidity of obtaining results, together with 
simplicity of equipment, is the main advantage of 
the photographic method. A failing is its low ac- 
curacy, difficulty of assigning numerical charac- 
teristics to the results obtained and, in all proba- 
bility, insufficient sensitivity when the structure 
contains atoms with sharply differing scattering 
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powers. The use of optical machines is apparently 
most widespread in France and Italy [23,24,25].In 
the Soviet Union, unfortunately, hardly anyone has 
been engaged in the construction of optical comput- 
ing systems. Only by V.A. Shamburov at the Insti- 
tute of Crystallography, Academy of Sciences of the 
USSR, have any developments been made. At the 
Institute of the Chemistry of Silicates, a machine 
acting on a similar principle has been made for the 
solution of rather different problems. 

VI. Certain electrical and electromagnetic 
machines acting on the analog principle are de- 
scribed in the literature. Except for the machine of 
Maxwell [26], all the rest are very similar in con- 
cept and differ mainly in constructional details 
(Rymer and Butler [27], Hagg and Laurent [28], 
Beevers and McEvan [29], Ramsay, Lipson, and 
Rogers [30], and Azaroff [31]). In all the machines 
listed, which are intended for summing one-dimen- 
sional Fourier series, the values of the trigono- 
metric function are obtained by the sinusoidal sub- 
division of resistances or transformers (voltage 
taken off cos 27x or sin 27x); multiplication by Fh is 
effected by potentiometric variation of the voltage, 
and summation by the addition of currents or volt- 
ages in corresponding circuits. 

Machines of this type are mechanized analogs 
of Robertson's method of strips, since the electric- 
al circuits only give the function F cos 27x directly 
at certain points ..x = 1/n,.2/n, 3/nssn/4neee 
Transformation from F cos 27x to Fcos 27hx is ef- 
fected by a corresponding change in the order of 
connecting the sections to the commutation panels. 
Each coil (or resistance) is joined to its own panel, 
and the order of connecting its sections corresponds 
to a definite h (for example, for h= 2, the zeroth, 
second, fourth, etc. sections are connected one 
after the other). Passage from one terminal to the 
next in turn is effected simultaneously and auto- 
matically on all the panels by means of the standard 
walking mechanism of a telephone apparatus. In this 
way, currents (or voltages) corresponding to 
Fy) cos 27hx for various h and given x are summed in 
a common circuit. 

The group of electrical analog machines for 
summing one-dimensional Fourier series includes 
one developed by engineers A. L. Dvorkin and E.G. 
Lipshtein for the X-Ray Laboratory of the INEOS, 
Academy of Sciences of the USSR. In this machine, 
the coefficients of the series are introduced by 
means of telephone dials; the maximum h index is 
24, and the period is divided into 96 parts. The in- 
troduction time for calculating a one-dimensional 
series is a minimum of 5 min, while operation and 
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calculation of the results takes 15 min. Correspond- 
ing to this, calculation of a two-dimensional series, 
in the simplest case of a projection with two lines 
of symmetry, with separation into 48 x 48 parts 
and 100 coefficients, requires approximately 10h 
which, generally speaking, is only a little less than 
the time required to make an analogous calculation 
by the hand method of strips. Undoubtedly, a change 
in the system of assigning the initial data and the 
attachment of a printing device at the output will 
enable the operating time to be reduced at least by 
a factor of two. 

VII. The second group of computers applic- 
able to x-ray structural analysis is made up of auto- 
matic systems which, according to their fundament- 
al parameters (rapidity of action, scale, complex- 
ity of servicing, and cost), are certainly intended 
only for the largest structural laboratories or com- 
puting centers. Before passing on to a description 
of individual types of machines,it is appropriate to 
discourse on some general principles of their con- 
struction [32-37]. 

All these machines are of the pulse type: each 
single-digit number is assigned in the form of a 
corresponding number of electrical impulses pass- 
ing successively through the circuit. The main 
operating units of the machines are automatic re- 
lays. Each relay can exist in two states, excited 
and quiescent, depending on whether it has re- 
ceived an impulse or not. Corresponding to this re- 
lay, an output circuit (or several of these) closes 
or opens and hence either directs pulses or pre- 
vents their direction into other parts of the machine. 
Thus, each relay fulfills a definite logical opera- 
tion, answering the question "yes" or "no" and con- 
veying this reply further. 

The relays of automatic machines are joined in 
individual groups or registers, playing the part of 
receivers or stores for the numbers. A number of 
pulses falling on a group of relays connected into a 
register bring it into a definite state "imaging" the 
corresponding digit. The single-digit registers are 
connected into groups for representing multidigit 
numbers. The connection is effected in such a way 
that, on completing a full "ten" in one register, a 
single pulse is developed and sent into the register 
of the following digit, all the units of the first 
register returning to the zero state. 

Clearly, questions of economy in the number of 
relays frequently compel one to change from the 
ordinary "scale-of-ten" to other systems. The most 
suitable (though not most economical) is the binary 
system of counting, which is widely used, especial- 
ly in electronic machines; in this, each "ten" con- 
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tains only two digits, 0 and 1, and hence requires 
only one relay (having two states). 

The successive direction of two groups of 
pulses into the same counters leads to the addition 
(or, if desired, subtraction) of the numbers applied 
to them. The pulses stored by the counters may be 
"held" in them for as long as required. In this re- 
spect, the counters play the part of memory cells. 
By feeding in a special group of pulses equal to a 
whole number of tens, hundreds, etc. , the number 
stored in the counter may be taken from it and 
passed to any other counter. 

The memory may be either temporary or perm- 
anent, i.e., the counter canbe cleared after opera- 
tion, or the stored number can be preserved inde- 
finitely. In machines intended for structural ana- 
lysis, the permanent "memory" cells are used to 
establish the values of the trigonometric functions. 

Thus, the counters may be used for logical 
operations, arithmetical operations, or as memory 
cells. Hence, the machine can carry out any in- 
structions in an assigned sequence. The problem 
consists only in decomposing a complex procedure 
into simple acts and in setting up a program of 
operations most suited to a machine of the particu- 
lar type. 

The introduction of numbers into a computer is 
usually effected by means of punched cards or tapes. 

On passing the cards (tapes) through the ma- 
chine, electrical contact is made only in places 
where there is a perforation, so that the corre- 
sponding circuit is closed and the right number of 
pulses is fed into the machine. 

As regards the physical principles of the opera- 
tion of the basic elements (relays), all impulse 
machines can be divided into three types: electro- 
magnetic, electronic, and magnetic. In the first, 
the instructions of the operator are carried out by 
electromagnetic relays opening or closing certain 
contacts; in the second, switching in the circuits is 
effected by feeding pulses to a network of electron- 
ic tubes connected in smallish groups (the so-called 
trigger circuits); in the third, the same operationis 
effected by changing the magnetic state of a certain 
magnetic element. Since the operation of electro- 
magnetic relays is associated with mechanical mo- 
tion, while that of electronic and magnetic devices 
involves merely a change of field, the computing 
speeds are very different. In electromagnetic ma- 
chines, the elementary operation takes fractions of 
a second, and in electronic machines fractions of a 
microsecond. 

The memorizing of numbers in pulse counters 
demands a very large number of relays. Hence, 
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electronic machines use other means of memoriz- 
ing; these are also based on electronic processes 
(cathode-ray tubes or mercury delay lines) and have 
a large capacity, together with a rapidity of action 
adequate to the computing speed of the electronic 
machine. In order to increase the memory of elec- 
tronic machines further, recording can be made on 
magnetic tapes or drums. 

Electromagnetic machines do not have this kind 
of memory structure. Hence, their capabilities are 
comparatively limited. Extra memory devices can 
be provided in the form of punched cards, punched 
by the machine itself at intermediate stages of the 
calculation. 

The transfer of such cards into the input de- 
vice of the machine, however, demands the inter- 
vention of an operator, so that the whole process of 
calculation is split into several stages and cannot be 
fully automated. 

VIII. Universal electromagnetic computers 
are known under the name of computer-analyzers. 
The genus of computer-analyzers includes not only 
the fundamental counting mechanisms or tabulators, 
but also a whole series of auxiliary devices: 
punchers, reproducers, controllers, sorting ma- 
chines, and multipliers. 

The first paper on the application of computer- 
analyzers to x-ray structural investigations was 
published in 1946 (Shaffer, Schomaker, and Pauling 
[88]), and further progress in computing processes 
was made, mainly in the United States and England, 
in 1946 to 1953 [39,40]. 

A general advantage of computer-analyzers is 
the possibility of using them for solving all three 
fundamental problems of structure analysis: the 
summation of Fourier series, the calculation of 
structure amplitudes, and the calculation of dif- 
ferential series for improving the coordinates. 

The simplest method of calculating Fourier 
series on computer-analyzers reduces in essence 
to the mechanization of the Beevers — Lipson strip 
summation. Each strip is replaced by a punched 
card containing the same data. Selection of the 
punched cards from the whole pack, their restack- 
ing, and the checking for correctness of choice are 
effected by hand. The tabulator only carries out 
the operation of adding the one-dimensional series. 
Two-dimensional and three-dimensional syntheses 
reduce to the successive calculation of one-dimen- 
sional series. In the Crystallographic Laboratory 
of Cambridge University, a block of punched-card 
strips has been prepared; this provides for the 
separation of a period into 120 parts and contains 
values of F}, with two significant figures, and 
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values of Fj, cos 27hx with four. The pack contains 
2460 different cards and, in all, including the 
necessary number of duplicates, 18,860 cards, 
disposed in eight boxes 60 x 19 x 10 cm in Size. 

According to the data of this laboratory [40], 
for the calculation of a two-dimensional series con- 
taining 200 terms, with division of a centrosym- 
metric cell into 60 x 60 parts, seven hours (for 
two operators) are required,i.e., roughly 7 sec 
per point; in calculating a three-dimensional 
series of 1000 terms some 50 sec per point are 
needed. 

This computing speed is quite inadequate for 
the systematic calculation of three-dimensional 
distributions. An advantage of this very simple 
method, however, is that it requires only one 
basic machine, the tabulator, which considerably 
broadens the practical possibility of setting up cal- 
culations and even the direct acquisitionby the 
structural laboratory of a ready-made machine, re- 
quiring no new development. 

Analysis shows that more than half the time 
required for calculating a series by the method de- 
scribed goes into the hand selection of the punched 
cards and their return to the pack, andalso into the 
checking for correct choice. Other, more complex 
methods of using computer-analyzers make it pos- 
sible to shorten this time substantially, and simul- 
taneously to reduce the required number of cards. 
In this case, however, practically the whole gamut 
of machines is required: puncher, controller, re- 
producer, multiplier, and sorting machine. Hence, 
the calculations can only be made in a computing 
center. 

For the calculation of structure amplitudes and 
the refinement of structures by differential synth- 
esis, the whole set of machines is required in all 
cases, the procedure for their use being quite com- 
plex and split into several stages. 

Despite the fact that the application of com- 
puter-analyzers to problems of structure analysis 
has become quite widespread in the United States, 
and especially in England, the method as a whole 
cannot be regarded as entirely satisfactory. Its 
chief disadvantages are: (1) relatively low gain in 
time (completely insufficient for three-dimensional 
calculations); (2) the need for a qualified operator 
(qualified as regards structural analysis) working 
in a general-purpose computing center; (3) insuffi- 
cient automation of the calculation for the com- 
plexity of its procedure, bringing with it the possi- 
bility of operator error. 

The method could be considerably improved by 
supplementing the computer-analyzer with elec- 
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tronic control for the sequence of operations. 

IX. Considerably more successful in the solu- 
tion of structural questions are specialized electro- 
magnetic-impulse machines, made up mainly of 
standard parts from computer-analyzers and tele- 
phone apparatus. 

In the Soviet Union, a machine of this type is 
being developed in the Physics Faculty of Moscow 
State University (under Professor A.N. Tikhonov 
and Engineers B.V. Sokolov and A.V. Luk'yanov). 
The machine is intended to sum Fourier series. 
The main unit of the machine (making it specialized) 
is a set of a special kind of electrical strips or 
permanent memory cells storing values of the trig- 
onometric functions for the first four periods. 

The second important element specializing the 
machine is a set of telephone walking mechanisms 
effecting simultaneous switching on the commuta- 
tion panels and transforming functions cos 27x into 
Cos 2Thx. 

The coefficients of the series are given in the 
form of four-figure numbers. The period is di- 
vided into 48 parts, the maximum index being 48. 
The authors estimate the machine time for calcu- 
lating a one-dimensional series as 2 min, the time 
for putting out in the form of punched cards 1.5 min, 
and the printing time, 2 min. Thus, the total time 
is of the order of 3.5 to 4 min, which means that 
a two-dimensional series (given in the form 


= 
pa F147, cos 2x hz cos 2xky) 
hk 


could be calculated in approximately 50 to 70 min. 
To this we must add the time for preparing the 
punched cards with the Fp, coefficients. Since a 
series consisting of 100 terms requires in all 10 to 
20 punched cards, this time is in general insignifi- 
cant (a trained operator can punch 250 cards per 
hour). 

The whole machine consists of puncher, read- 
ing mechanism, main computing mechanism, auto~ 
matic puncher, and printing system. 

The machine as a whole is to be considered as 
intermediate between a laboratory-only machine 
and one belonging to a computing center. In this lies 
its main defect. For a laboratory-only machine it is 
too cumbersome, and probably too expensive; in re- 
lation to machines of the second type, it is insuffi- 
ciently versatile (solving only one problem), and 
gives altogether too low a computing speed. Fur- 
ther development of the machine is to be desired, 
since the general advantage of electromagnetic over 
electronic computers is in reliability of operation, 
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lower cost, and relative simplicity of servicing. 
Improvements should be directed toward increas- 
ing the field of problems which the machine can 
solve, and also the computing speed. There are 
certainly possibilities in this direction, as can be 
seen from the example of another specialized im- 
pulse machine constructed by Beck and Cormarie 
in France, using more fully the operating prin- 
ciples of automatic telephone exchanges. The ma- 
chine simultaneously gives the values of p(xyz) at 
all (64) points x for given values of y and z. 

X. For all their advantage in reliability of 
operation, relative simplicity of servicing, and 
relative cheapness, electromagnetic computers 
cannot emulate electronic types in computing speed 
and flexibility in setting up problems. 

All the demands of the most laborious and 
variegated structural calculations can be complete- 
ly met, given the existence of a universal high- 
speed electronic computer running continuously 
and accessible for the customers of the computing 
center. 

Much has been said of the applications of uni- 
versal electronic computers [50-58], in particular 
in contributions to the International Congress of 
Crystallographers in Paris in 1954 [54-57]. The 
possibilities of electronic machines may be judged 
from the data presented by Cruikshank [57]: fora 
total machine time of 290 h, the Manchester 
University machine made calculations for 18 struc- 
tures, including 56 calculations of structure ampli- 
tudes and 50 differential syntheses. 

The application of a universal electronic com- 
puter for solving structural problems in Moscow 
was effected by Fellow of the Faculty of Computer 
Mathematics, Moscow State University, N. P. 
Trifonov (Faculty Head Academician S. L. Sobolev) 
in collaboration with others of the Crystallochemi- 
cal Laboratory, Institute of General and Inorganic 
Chemistry, Academy of Sciences of the USSR. 

The main preparatory work connected with the 
setting up of structural problems on the universal 
electronic computer is the programming of the 
problem, i.e. , finding the most convenient and 
economical sequence of operations and writing this 
in the form of coded machine orders. Programming 
must take into account not only the possible econ- 
omy in machine time, but also the loading of the 
machine "memory," the control of calculations at 
various stages, and the suitability of the setting up 
of the problem in view of its possible modification 
or complication at later stages, etc. 

N.P. Trifonov set up and solved the following 
problems: (1) calculation of structure amplitudes 
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from "trial" coordinates of the atoms and atomic 
factors; (2) calculation of electron-density distribu- 
tions; (3) linking the two problems into a single 
cycle to be carried out automatically, i.e., finding 
Fealc from Xj, Vjo js and f;, transferring the sign 
of Fealce (or in general the initial phase Dgajc) to 
the experimental Fexp, and calculating the electron- 
density distribution; (4) refining the coordinates of 
the atoms by the method of differential synthesis, 
uniting all the successive-approximation cycles in- 
to a single automatically executed process. 

With respect to the two-dimensional problem, 
programs were created for all plane symmetry 
groups; with respect to the three-dimensional prob- 
lem this was done for space groups C}), and Di 
The programming of the problem for the other space 
groups, however, does not now present any tech- 
nical difficulties. 

Let us consider a few details. 

1. The calculation of structure amplitudes is 
effected with allowance for the dependence of the 
atomic factors on the hkl indices. The atomic fac- 
tors are given to the machine in the form of a table 
of values of Sj as a function of sin’v/A with a 0.1 in- 
terval in sin’ /2X. The machine is also given the 
parameters of the quadratic form (sin’/A )? = Ah? 
+Bk? +Cl? +Dhk +Ehl + Fki , and it calculates 
sin 3/A for each hkl reflection and then finds the 
value of f j kd) by interpolation of the tabulated 
data. Thus, another problem of independent interest 
is solved in passing: the calculation of the quadratic 
form. 

2. In the calculations of electron-density dis- 
tributions made on the machine, the period was di- 
vided into 60 parts. In principle, however,the elec- 
tronic machine, not being connected with strips, 
can just as easily effect division into arbitrary 
parts, and in practice the most convenient are 
hundredths of the period. 

3. The solution of problems on the universal 
machine is so flexible that it permits any auxiliary 
calculations to be included at will. Thus, at any 
stage in the refinement of coordinates by the me- 
thod of differential synthesis one has the possi- 
bility of switching in a program for calculating the 
whole electron-density distribution in order to in- 
spect the general picture. It is possible to include 
in the process of successive refinements not all, 
but only some of the atoms, leaving the coordinates 
of the rest unchanged or disregarding them entire- 
ly, only including them in the general process at 
some intermediate stage, when the process of re- 
finement has reached the required level (this opera- 
tion does not require a halt in the operation of the 
machine). 
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The machine solves in passing the problem of 
comparing Fealc and Fexp, and reducing Fexp to 
the absolute scale. 

At the concluding stage, a calculation is made 
for the "theoretical series" in order to find the cor- 
rection for the breaking off [of the series] (Booth) 
and to calculate the functions 24,7) (Fexch — Foalo> 
Chk b? (Fexch — Fealc), etc., which are needed in 
order to estimate the accuracy of determination of 
the coordinates from the well-known computing 
formulas. 

For the characteristics of the machine time 
necessary, it is most convenient to take specific 
examples. 

In studying crystals of Ni(NH3)3(SCN)., a uni- 
versal electronic computer was used to calculate 
the ordinary and weighted projections on the XY 
plane. For 250 independent hk0 reflections and a 
plane group of electron density pgg, calculation of 
the projection by the hand method of strips would 
require approximately 30 h. Roughly the same 
time would go into calculating the signs of the struc- 
ture amplitudes from tables or nomograms. Thus, 
a single cycle of refinement by hand operation 
would require around 60 h. On this problem the 
universal electronic computer spends 50 min ma- 
chine time if the programs for calculating the struc- 
ture amplitudes and electron density are not con- 
nected together, and in all, 30 min on joining these 
into a single cycle (structure-amplitude values not 
printed). Use of the electronic machine, however, 
requires preliminary treatment of the initial data: 
transforming the data to the binary system and put- 
ting it on punched tape. Some 2 h are spent in this 
operation for a problem of the size in question. 
Furthermore, for nine independent atoms in the 
cell, not less than an hour is spent on interpolation 
calculations for finding the coordinates of the 
maxima. Hence, in all, the first cycle of refine- 
ment requires some three and a half hours. On us- 
ing the method of differential synthesis, as a result 
of the curtailment of print-out operations, the ma- 
chine time for one cycle of the present problem is 
reduced to 5 min. The elimination of the interpola- 
tion work and administrative operations in trans- 
ferring the material for analysis to the laboratory 
after each cycle of refinement reduces the whole 
working time to simply that of the machines, plus 
the time required for preliminary preparation of 
the initial data. In the present problem, eight 
cycles of refinement were used, and these took the 
machine 40 min. 

The whole calculation thus required 2 h 40 min 
instead of the 14 h needed when the F-series me- 
thod was used in the same machine, or 480 h when 
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using strips and nomograms. 

The refinement of atomic coordinates in the 
Cs,CoCl, structure was effected by the differential- 
synthesis method, applied to a three-dimensional 
density distribution. The number of independent 
hkl reflections was 950; the space group was Djf, 
= Pnam. The number of structurally independent 
atoms was six. The computing time for one cycle 
was 30 min. Three cycles of refinement were 
needed in all. In an hour anda half, this gave re- 
sults which would have called for months by the 
method of strips and at least a week on any of the 
other machines mentioned above. 

Work on broadening the field of problems 
soluble on the universal machine is continuing. At 
the present time work is proceeding with the set- 
ting up of another of the most important structural 
problems: the choice of combinations of reflec - 
tions hk/, h'k'?', h+h', k+k', 1 + Z' for determin- 
ing the signs of the structure amplitudes by the 
statistical method. The same problem has been set 
up in England by Cochran and Douglas [58], but the 
solution is very cumbersome and not entirely con- 
vincing. The setting up of the sign relationships on 
a different machine is also occupying a group of 
engineers associated with the x-ray laboratory of 
the INEOS, Academy of Sciences of the USSR. 

It follows from this discussion that the use of 
universal electronic computers has great possi- 
bilities, especially with the continually changing 
demands of x-ray structural analysis, with the de- 
velopment of new methods of solving structural 
problems, and with modifications of old ones. 

A fundamental failing in the use of universal 
electronic computers is the time occupied by the 
preparatory operations in setting the initial data 
and the final operations in printing the results, this 
being because the standard auxiliary apparatus of 
the universal machine used for these purposes can- 
not take account of the specifications of the prob- 
lem (such as the fact that Fexp is given to two or 
at most three significant figures, and that the re- 
sults, the coordinates of the atoms, require no 
more than three places of decimals). 

In this respect, as in certain others, special- 
ized discrete-action electronic computers shouldbe 
more suitable. Specialization of the machine should 
substantially reduce its size and cost, and gain time 
mainly in the preparatory and final operations. The 
capabilities of such a machine will however fall sig- 
nificantly below those of the universal type. Ma- 
chines of this type have evidently not been developed 
abroad. 
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In the Soviet Union, on the initiative of G.S. 
Zhdanov, it is proposed to construct a machine for 
solving three problems: finding the electron-den- 
sity distribution in two-dimensional cross sections, 
finding F(hkl) for a given 1, and calculating the 
quadratic form. 

XI. Another type of specialized electronic ma- 
chine is the Pepinsky X-RAC [59], an analog ma- 
chine having in principle much in common with the 
photographic systems based on the superposition of 
two-dimensional sinusoidal distributions. The idea 
of the machine arose from considering the mecha- 
nism for the transmission of television signals: 
that of converting the picture into electrical signals 
and conversely, the signals into an image. Since 
the taking of a motionless image constitutes a 
periodically repeatable process, the electrical sig- 
nal produced may be represented as a function of 
time in the form of a double Fourier series, and 
its transformation into an image on the screen as 
a translation of a function of time into a function of 
the coordinates of a point on the screen. The beam 
of an oscilloscope oscillates in two directions, ver- 
tical and horizontal, at constant frequencies Wyert 
and Whoriz, Where Wporiz >> “vert: The unmodu- 
lated beam will trace out a continuous band on the 
screen of the cathode-ray tube; this moves along 
the horizontal and is displaced along the vertical. 

If a modulated sinusoidal voltage at frequency 
hhoriz, is applied to the grid of the CRT, the 
sereen shows vertical bands (h bands over the 
width of the screen). If the frequency of the applied 
field equals kwyert, horizontal bands develop (k 
bands over the height of the screen). On applying 
frequency hwhoriz +kwvert, corresponding in- 
clined bands are obtained. The contrast of these 
will depend on the amplitude of the applied field. 
Clearly, on applying a field equal to 24;,.Fy, cos 27 
* M&poriz tk’yert), the oscilloscope screen will 
give an image of the electron-density distribution. 
Thus, the main problem lies in obtaining an elec- 
trical signal constituting the sum of a double 
Fourier series. This calls for the generation of 
the necessary set of frequencies in two ranges (low- 
frequency kwyert and high-frequency hw} oyjz), the 
assignment of their amplitudes and initial phases, 
and then their superposition. In this review, the 
technical details of constructing such a complex 
machine cannot be given; the size would be very 
great and so, undoubtedly, would be the cost. 

The combination of the basic X-RAC machine 
with the auxiliary system S-FAC [60] for calculat- 
ing structure amplitudes [61] makes the Pepinsky 
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machine certainly the most flexible of all the spe- 
cialized machines. Such a machine can solve any 
structural problems, beginning with the summation 
of series and finishing with minimization of the in- 
teratomic function and the calculation of sign rela- 
tionships. 

At the same time, one cannot fail to note the 
economic unsuitability of manufacturing such a spe- 
cialized machine. Its cost and servicing approach 
those of the universal machine; its functions are 
substantially narrower; it would be very problem- 
atical to ensure constancy of loading for the enorm- 
ous operating speeds. One may well suppose that 
the correct way toward the coordinated solution of 
structural problems lies not in the creation of elec- 
tronic machines such as the X-RAC, but in the de- 
velopment of specialized discrete-acting magnetic 
systems [63-65] which in the near future should 
compete substantially with electronic machines. 

XII. The first communication regarding the 
application of a magnetic machine of the universal 
type, MAGYC, to the solution of structural-analy- 
sis problems was given by Booth at the Paris Con- 
gress [62]. 

The replacement of electronic trigger circuits 
by magnetic devices evidently enables the scale and 
cost of machines to be sharply reduced and the re- 
liability of the work to be increased. Reducing the 
size and cost in turn makes it possible to raise the 
number of computing units and modify the general 
structure of the machine, i.e., to construct spe- 
cialized machines having almost the same capa- 
bilities, as regards modification and broadening of 
the field of problems soluble, as those of a uni- 
versal electronic computer. 

It is also not excluded that table machines for 
the solution of the simplest structural problems 
(summation of one-dimensional series) will prove 
most compact, inexpensive, and reliable on using 
magnetic counting units. So far, however, there 
are no real data on this. 

XII. By comparing all the available data on 
computers for the structural analysis of crystals, 
we may draw the following conclusions. 

1. As comparatively cheap and compact mech- 
anisms, suitable for use in laboratory conditions, 
the most appropriate are: (a) purely mechanical 
machines for summing one-dimensional series; (b) 
an optical machine of the Eller type for rapid (but 
not exact) calculation of two-dimensional density 
distributions and for computing structure ampli- 
tudes. 

2. For the solution of complex problems of 
structural analysis, it is appropriate either to use 


universal electronic machines in a general-purpose 
computing center, or to create specialized but fair- 
ly versatile machines for the structure-computing 
center, based on electronic trigger circuits or mag- 
netic counting units. 


LITERATURE CITED 


I. Mechanical Systems 


1. R.C. Evans and H.S. Peiser, "A machine for 
the computation of structure factors," Proc. 
Phys. Soc. ,54, 457 (1942). 

2. A.D. Booth, Fourier Technique in X-Ray Or- 
ganic Structure Analysis (1948), p. 67. 

3. A.D. Booth, "Two calculating machines for x- 
ray crystal structure analysis," J. Appl. Phys. 
18, 664 (1947). 

4. W.L. Bragg, "A device for calculating struc- 

ture factors," Acta Cryst.,5, 4, 474 (1952). 

Vv. Vand, "A simple mechanical structure-fac- 

tor computing aid, " Acta Cryst.,5, 3, 390 

(1952). a 

6. L. Clastre and R. Gay, "Machine a calculer 
les facteurs de structure," Acta Cryst. ,7, 10, 
623 (1954). rp 

7. S.L. Brown and L.L. Wheeler, "The use of a 
mechanical synthesizer to solve trigonometric 
equations," J. Appl. Phys.,14, 30 (1943). 

8. M.J. Robertson, "A simple harmonic continu- 
ous calculating machine," Phil. Mag. ,13, 413 
(1932). <j 

9. A. Shilton, "A machine for harmonic synthesis," 
Proc. Phys. Soc. ,56, 130 (1944). 

10. Dan McLachlan, Jr., and E. F. Champaygne, 
"A machine for the application of sand in mak- 
ing Fourier projections of crystal structures," 
J. Appl. Phys.,17, 1006 (1946). 

11. A.J. Rose, "Machine a calculer permettant la 
determination de fonctions periodiques et leur 
introduction dans des calculs, " Centre Nat. 
Recherche Sci. Lab. de Bellevue (Paris, 1950) 
Walls We 


12. V. Vand, "A Fourier electron-density bal- 
ances UJ acoCl wINSth. 20 elon Looe 

13. C.A. Beevers and J.M. Robertson, Collection: 
Computing Methods and the Phase Problem in 
X-Ray Crystal Analysis (X-Ray Crystal 
Analysis Lab., Pennsylvania State College, 
1952), p. 119. 

14. A.D. Booth, Fourier Technique in X-Ray 
Organic Structure Analysis (1948), p. 76. 

15. J.M. Robertson, "A desk computer for Fourier 
analysis, " Acta Cryst.,7, 624 (1954). 


o 


> 


Il. 


REVIEW 203 


Machines Based on an Optical 


Analogy 


16. 


Gs 


18. 


its) 


20. 


21. 


22. 


23. 


24, 


C.W. Bunn, Chemical Crystallography (Ox- 
ford, 1946), p. 271. 

H. Lipson and W. Cochran, The Determination 
of Crystal Structures (1953), pp. 140, 236. 
C.R. Eller, "Le photosommateur harmonique 
et ses possibilités," Bull. Soc. Franc. Min. 
Crist? -703 a TA1955)5 

C.R. Eller, "Sur un nouvel appareil pour le 
developement par voie optique des series de 
Fourier," Compt. Rend. Acad. Sci. , 232, 1122, 
2333 (1931). sus 

L. Bru, M. P. Rodriguez, and M. Cubero, 
"Application of Eller's optical machine to the 
determination of crystal structures by the di- 
rect method," Acta Cryst., 7, 625 (1954). 

L. Bru, M.P. Rodriguez, and A. L. L. Castro, 
"Application of Eller's optical machine to the 
determination of crystal structures by the 
trial-and-error method," Acta Cryst.,7, 625 
(1954). ba 

L. Bru, M.P. Rodriguez, and M.G. Gea, 
"Application of Eller's optical machine to the 
refinement of crystal structure," Acta Cryst., 
7, 625 (1954). 

M.M. Wolfson, "A photoelectric structure- 
factor machine," Acta Cryst. ,4, 3, 250 (1951). 
W. Hughes and C.A. Taylor, "Apparatus used 
in the development of optical-diffraction me- 
thods for the solution of problems in x-ray 
analysis, " J. Sci. Instr. ,30, 4, 105 (1953). 

R. Fiirth and R.W. Pringle, "A new photoelec- 
tric method for synthesis and analysis, " Phil. 
Mag.,35, 643 (1944). 


Electrical Analog Machines 


26. 


Pg be 


28. 


Ze > 


L.R. Maxwell, "An electrical method for com- 
pounding sine functions," Rev. Sci. Instr. ,11, 
47 (1940). — 
T.B. Rymer and C.C. Butler, "An electrical 
circuit for harmonic analysis and other calcu- 
lations," Phil. Mag.,35, 606 (1944). 

G. Hagg and T. Laurent, "Machine for the 
summation of Fourier series," J. Sci. Instr., 
23, 155 (1946). 

C.A. Beevers and D. McEvan, "A machine for 
the rapid summation of Fourier series," J. Sci. 
Instr. ,19, 150 (1942). 

I. W. Ramsay, H. Lipson, and D. Rogers, 
Collection: Computing Methods and the Phase 
Problem in X-Ray Crystal Analysis (X-Ray 
Crystal Analysis Lab., Pennsylvania State 


College, 1952), p.-119. 

31. L.V. Azaroff, "A one-dimensional Fourier 
analog computer," Rev. Sci. Instr.,25, 471 
(1954). me 


JEN E> Genera! Reviews on the Prin- 


ciples of Constructing Pulse Com- 


puters 


32. H.H. Aiken and G.M. Hopper, Automatically 
Controlled Calculating Machine [Russian trans- 
lation from the English by M.A. Bykhovskii, 
Usp. Mat. Nauk, 4, 119 (1948)]. 

33. M.L. Bykhovskii, "Fundamentals of electron- 
ic mathematical discrete-counting machines," 
Usp. Mat. Nauk, 3, 69 (1949). 

34. L.D. Kudryavtsev, "Principles of producing 
arithmetical operations on calculating ma- 
chines," Usp. Mat. Nauk, 3, 104 (1950). 

35. Murray, Theory of Mathematical Machines 
[Russian translation] (Inizdat, 1949). 

36. Hartree, Modern Development of Computing 
Machines [Russian translation] (GNB). 

37. I. Ya. Akushskii, "Brief sketch of computer- 
analyzers," Izv. otd. Tekh. Nauk AN SSSR, 8, 
1081 (1946). - 


V. Application of Computer-Analyzers 


38. P.A. Shaffer, V. Schomaker, and L. Pauling, 
"The use of punched cards in molecular struc- 
ture determinations," J. Chem. Phys. ,14, 648 
(1946). re 

39. E.G. Cox, L.Gross, and G.A. Jeffrey, "Use 
of punched-card tabulating machines for crys- 
tallographic Fourier syntheses," Nature, 159, 
433 (1947). 

40. M.L. Hodgson, C.J.B. Clew, and W. Cochran, 
"A punched-card modification of the Beevers — 
Lipson method of Fourier synthesis," Acta 
Cryst. ,2, 113 (1949). 

41. EB! GoCoxtand’ Gi A. Jeffrey, "The use of 
'Hollerith' computing equipment in crystal- 
structure analysis," Acta Cryst. ,2, 341 (1949). 

42. J. Donohue and V. Schomaker, "The use of 
punched cards in molecular structure deter- 
minations. III. Structure-factor calculations 
of x-ray crystaNography," Acta Cryst. ,2, 344 
(1949). a 

43. M.D. Grems and J.S. Kasper, "An improved 
punched-card method for crystal-structure 
factor calculations," Acta Cryst. ,2, 6,347 
(1949). as 

44, E.G. Cox, L. Gross, and G.A. Jeffrey, "A 
Hollerith technique for computing three-dimen- 


204 


45. 


46. 


47. 


48. 


WOle 


REVIEW 


sional differential Fourier syntheses in x-ray 
crystal-structure analysis," Acta Cryst. ,2, 
351 (1949). = 
E. Cox, Collection: Computing Methods and 
the Phase Problem in X-Ray Crystal Analysis 
(X-Ray Crystal Analysis Lab., Pennsylvania 
State College, 1952), p. 1382. 

E.W. Hughes, Collection: Computing Methods 
and the Phase Problem in X-Ray Crystal- 
Analysis (X-Ray Crystal Analysis Lab. , 
Pennsylvania State College, 1952), p. 141. 

N. Kitz and B. Marchington, "A method of 
Fourier synthesis using a standard Hollerith 
senior rolling total tabulator," Acta Cryst., 
6, o20 (1993), 

C.R. Hudgexs and A.M. Ross, "Computing 
Fourier syntheses in x-ray crystal-structure 
analysis," Anal. Chem. ,25, 734 (1953). 


Application of Universal 


Blectronic Computers 


49, 


ol. 


52. 


J.M. Bennett and J.C. Kendrew, "The com- 
putation of Fourier syntheses with a digital 
electronic calculating machine," Acta Cryst., 
5, 109 (1952). 

F. Ordway, Collection: Computing Methods 
and Phase Problem in X-Ray Crystal Analysis 
(X-Ray Crystal Analysis Lab. , Pennsylvania 
State College, 1952), p. 148. 

S.W. Wayer and K.N. Trueblood, "Three- 
dimensional Fourier summations on a high- 
speed digital computer," Acta Cryst., 6, 427 
(1953). * 

F.R. Ahmed and D. Cruickshank, "Crystallo- 
graphic calculations on the Manchester 
University electronic digital computer," Acta 
Cryst.,'6; 765 (1953). 

Els Thompson and D. Caminer, "The use of a 
high-speed automatic calculator in the refine- 
ment stages of crystal-structure determina- 
tion, " Acta Cryst., 7, 260 (1954). 


A.S. Douglas, "The use of digital computers 
in structure analysis," Acta Cryst., 7, 624 


A.S. Douglas,"Application of the EDSAC to 
structure determination by a direct method," 


M.M. Woolfson and W. Cochran, "Some ex- 
perience of direct methods of Fourier struc- 


D.W.J. Cruickshank,”The use of Manchester 
University electronic digital computer in de- 
tailed structure refinement, " Acta Cryst. 7, 


W. Cochran and A.S. Douglas, "The use of a 
high-speed digital computer for the direct de- 
termination of crystal structures," Proc. Roy. 


54. 
(1954). 
Doe 
Acta Cryst. ,7, 624 (1954). 
56 
ture determination, " Acta Cryst. 7, 624 
(1954). 
sin 
624 (1954). 
58. 
Soc, A227, 1171,,486601955).. 
VII. Pepinsky Machine 
ag) 


60. 


R. Pepinsky, Collection: Computing Methods 
and Phase Problems in X-Ray Crystal Analysis 
(X-Ray Crystal Analysis Lab., Pennsylvania 
State College, 1952), p. 167. 

R. Pepinsky, "R-PAC: Analog recorder-play- 
back computer for crystal analysis," J. Appl. 
Phys. , 24, 663 (1953). 


61. R. Pepinsky, "Letters to the Editor," J. Appl. 
Phys. ,24, 1520 (1953). 

VIII. Magnetic Machines 

62. A.D. Booth, "The application of MAGIC to the 


63. 


64. 


65. 


determination of crystal structures," Acta 
Cryst. , 7, 626 (1954). 

An. Vang, "Static magnetic storage and delay 
line," J. Appl. Phys. ,21, 49 (1950). 

An. Vang, "Magnetic delay-line storage, " 
Proc. Inst. Radio Eng., 38, 626 (1950). 

An. Vang, "Magnetic triggers," Proc. Inst. 
Radio Eng., 39, No. 4 (1951). 


